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CHAPTER  1 


INTRODUCTION 


1-1  Motivation  and  Previous  Research 

A  theoretical  investigation  of  the  nonlinear  effects  in  the  reflection  and 
refraction  of  plane  finite-amplitude  acoustic  waves  at  the  interface  between  two 
lossless  fluids  is  presented  in  this  dissertation.  This  research  may  in  one  sense 
be  considered  an  outgrowth  of  the  author’s  master’s  thesis  work  (Cotaras  1985), 
which  dealt  with  finite-amplitude  effects  in  long  range  underwater  acoustic  prop¬ 
agation.  That  work  is  rigorously  valid  only  up  to  the  point  that  the  acoustic 
waves  interact  with  either  the  surface  or  the  bottom  of  the  ocean  or  pass  through 
a  caustic  (focal  region).  Although  this  work  does  not  contain  any  discussion  of 
what  happens  at  a  caustic,  it  does  have  application  to  the  reflection  and  refraction 
of  finite-amplitude  sound  at  both  the  surface  and  the  bottom  of  the  ocean. 

Evidence  from  both  inside  and  outside  the  field  of  nonlinear  acoustics 
indicates  that  the  laws  governing  the  reflection  and  refraction  of  finite-amplitude 
sound  may  be  different  from  the  well-known  linear  theory  results — Snell’s  law  and 
the  law  of  specular  reflection.  From  outside  the  field  of  acoustics  come  the  unusual 
effects  that  are  observed  during  the  reflection  and  refraction  of  plane  shock  waves, 
effects  such  as  the  stem  shock.1  These  effects  have  been  studied  for  some  time 
and  are  still  of  interest  today  [see,  for  example,  Polachek  and  Seeger  (1951); 
Jahn  (1956,  1957);  Henderson  (1966);  Miura  and  Glass  (1983)].  The  results  of 
the  research  are  not,  however,  of  direct  use  in  acoustics  for  the  following  reasons: 
(1)  The  pressure  change  at  the  shock  is  typically  much  larger  than  that  in  finite- 
amplitude  acoustics,  and  (2)  the  results  indicate  only  what  happens  at  the  shock, 
not  what  happens  to  a  continuous  signal.  From  within  the  field  of  nonlinear 

1  Upon  reflection  from  an  interface,  the  intersection  point  of  the  incident  and  reflected  shocks 
detaches  from  the  interface.  The  intersection  point  is,  however,  still  ‘connected’  to  the  interface 
by  a  third  shock — the  stem  shock.  For  more  details,  see,  for  example,  Courant  and  Friedrichs 
(1976,  |>.  334). 
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acoustics,  evidence  of  a  possible  difference  from  linear  theory  is  suggested  by 
Blackstock  (1959).  Studying  the  normal  reflection  of  finite-amplitude  sound  from 
a  rigid  wall,  Blackstock  noted  that  the  iinear  theory  result  of  pressure  doubling 
is  not  rigorously  correct.  Pressure  doubling  is  actually  the  small-signal  limit  of 
the  proper  law,  which  is  that  excess  sound  speed  (c  —  Cq)  doubles.2 

Despite  the  evidence  of  a  possible  deviation  from  the  linear  theory,  the 
problem  of  reflection  and  refraction  of  finite-amplitude  sound  has  remained  al¬ 
most  untouched.  Van  Buren  and  Breazeale  (1968)  dealt  with  the  reflection  of 
finite-amplitude  ultrasonic  waves  by  decomposing  the  incident  finite-amplitude 
signal  into  its  Fourier  components  and  then  reflecting  each  component  indepen¬ 
dently  using  linear  theory.  In  their  procedure  the  reflected  components  are  then 
recombined  to  reconstitute  the  reflected  wave.  Thus  they  did  not  directly  address 
the  quest  ion  of  nonlinear  effects  in  reflection  and  refraction. 

Two  Chinese  reseachers  have  recently  published  pertinent  work:  Qian 
(1982)  on  reflection  of  finite-amplitude  acoustic  waves  and  Feng  (1983)  on  re¬ 
flection  and  refraction.  Both  use  perturbation  methods  and  Lagrangian  coor¬ 
dinates.  Qian  examined  that  which  is  a  special  case  for  us — the  reflection  of 
finite-amplitude  sound  from  a  rigid  wall  at  oblique  incidence.  The  paper  has  a 
few  algebraic  errors3  but  is  otherwise  correct  for  the  case  it  treats.  On  the  other 
hand,  the  work  by  Feng  appears  to  contain  some  fundamental  errors  in  both  the 
wave  equations  and  the  boundary  conditions.  After  expanding  all  the  acoustic 
variables  in  a  small  parameter  e,  Feng  obtains  the  0(e)  and  0(e2)  wave  equations 
and  boundary  conditions.  However,  Feng  then  assumes  that  the  interaction  terms 
in  the  0(c2)  wave  equations  [his  Eqs.  (13)  and  (15)]  are  zero.  This  is  true  only  for 
the  case  of  simple  wave  motion.4  For  compound  wave  flow,  which  is  what  occurs 
on  the  fluid  I  side  (the  incident  and  reflected  fields  co-exist),  the  terms  do  not 
cancel.  Feng’s  wave  equations  are,  therefore,  not  valid  on  the  fluid  I  side.  Fur¬ 
thermore,  in  his  0(e2)  particle  velocity  boundary  condition,  Feng  fails  to  account 
for  the  variation  of  the  normal  to  the  interface.  It  is  shown  in  this  dissertation 
that  the  variation  vanishes  only  in  the  special  case  of  normal  incidence. 

An  earlier  presentation  of  some  of  the  results  of  this  research  (Cotaras 
and  Blackstock  1987)  apparently  stimulated  some  interest  >n  the  problem.  For 

2For  finite-amplitude  waves,  the  sound  speed  c  differs  from  the  small-signal  sound  speed  c o: 
see  the  next  section. 

3For  example,  in  Eq.  (13a)  an  unneccesary  l/(cg)2,  where  cj,  is  the  small-signal  sound  speed 
in  fluid  I,  is  introduced;  this  error  is  repeated  in  Eq.  (19). 

4 For  the  case  of  plane  waves,  simple  wave  flow  may  be  defined  as  a  wave  field  that  consists 
of  progressive  waves  only,  that  is,  waves  propagating  in  one  direction  only.  If  progressive  and 
regressive  waves  co-exist,  then  the  wave  field  is  referred  to  as  compound. 
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example,  Shu  and  Ginsberg  (1988)  are  examining  the  reflection  and  refraction  of 
finite-amplitude  sound  at  a  fluid-solid  interface.  At  the  time  of  writing  of  this 
dissertation,  however,  their  work  was  not  available  for  examination. 

1-2  The  Propagation  Speed  of  Finite-Amplitude  Acoustic  Waves5 

It  is  well  known  that  the  propagation  of  a  finite-amplitude  wave  can¬ 
not  be  accurately  modeled  by  a  small-signal  acoustical  theory;  see,  for  example, 
Blackstock  (1972).  Small-signal  theory  fails  because  the  propagation  speed  of 
a  finite-amplitude  wave  depends  on  the  local  particle  velocity.  The  dependence 
arises  in  two  ways.  (1)  Convection:  Convection  occurs  when  the  fluid  particles 
themselves  are  set  into  motion  by  the  passing  acoustic  wave  and  contribute  their 
own  velocity  to  the  total  wave  speed.  Thus  the  actual  propagation  speed  dx/dt 
may  be  expressed  as 

dx 

dt=c  +  u  '  (u) 

where  c  is  the  sound  speed  (which  is  to  be  distinguished  from  the  propagation 
speed).  (2)  The  nonlinearity  of  the  pressure-density  relationship:  Because  of  this 
nonlinearity,  the  sound  speed  c  is  not,  in  general,  a  constant.  For  an  outgoing 
plane  wave,  the  sound  speed  is  given  by  the  following  [see,  for  example,  Beyer 
(1974)]: 

B 

c=co  +  — u  ,  (1.2) 

where  Co  is  a  true  constant  referred  to  as  the  small-signal  sound  speed  (the  value 
that  appears  in  tables),  and  A  and  B  are  the  first  and  second  coefficients  of 
the  Taylor  series  expansion  of  the  pressure  density  relation.  When  the  above 
equations  are  combined,  the  resulting  equation  is  usually  written  as 

—  =  Cq  +  0u  ,  (1.3) 

where  (3  =  1  -f  B/2A  is  called  the  coefficient  of  nonlinearity  and  has  the  value  1.2 
for  air  and  3.5  for  water.  Equation  (1.3)  is  sometimes  called  the  wavelet  speed, 
a  wavelet  being  a  given  point  on  the  waveform  of  a  propagating  wave. 

5 For  a  more  thorough  introduction  to  the  physical  origins  of  nonlinear  acoustics  and  its 
applications,  see,  for  example,  Hamilton  (1986). 
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1-3  Scope  of  the  Research 

The  primary  objective  of  this  dissertation  is  to  investigate  nonlinear 
effects  in  reflection  and  refraction  at  a  plane  interface  between  two  lossless  flu¬ 
ids.  However,  because  losses  are  present  in  all  real  fluids  and  because  losses 
must  eventually  be  accounted  for  in  a  complete  treatment  of  reflection  and 
refraction,  the  equation  for  finite-amplitude  wave  motion  in  a  homogeneous, 
thermoviscous  fluid  with  a  single  relaxation  mechanism  is  derived  in  this  disser¬ 
tation.  Although  thermoviscous  effects,  relaxation  effects,  and  nonlinear  effects 
have  been  studied  separately,  the  derivation  (from  first  principles)  of  the  wave 
equation  that  accounts  for  all  three  effects  has  not  been  previously  attempted. 
It  is  therefore  included  in  this  dissertation. 

The  dissertation  is  divided  into  three  parts.  Presented  in  the  first  part 
are  the  basic  equations  for  a  homogeneous,  thermoviscous  fluid  with  a  single 
relaxation  mechanism  and  a  method  for  ranking  the  terms  in  the  equations. 
The  fluid  is  assumed  to  be  initially  quiet,  irrotational,  uniform,  and  in  thermo¬ 
dynamic  equilibrium.  Moreover,  effects  of  rotational  flow  are  ignored.  When 
presenting  the  basic  equations,  we  make  a  fundamental  assumption:  The  mag¬ 
nitude  of  the  deviation  from  thermodynamic  equilibrium  is  assumed  to  be  small, 
and  linear  relationships  between  each  thermodynamic  flux  and  all  the  thermo¬ 
dynamic  forces  are  assumed  to  hold.  The  possibility  of  cross-effects  between 
the  different  thermodynamic  fluxes  and  forces  is  pointed  out,  but  neglected. 
Ranking  the  terms  in  the  basic  equations  according  to  their  relative  impor¬ 
tance  requires  that  additional  fundamental  assumptions  be  made  about  (1)  the 
amplitude  of  the  acoustic  signal,  (2)  the  magnitude  of  the  transport  coefficients, 
and  (3)  the  magnitude  of  the  dispersion  caused  by  the  relaxation.  Use  of  the 
ranking  system  enables  us  to  develop  simplified  forms  of  the  equations.  The 
wave  equation  for  finite- amplitude  signals  in  a  thermoviscous  fluid  with  a  single 
relaxation  mechanism  is  then  developed. 

Analyzed  in  the  second  part  of  this  dissertation  is  the  reflection  and 
refraction  of  finite-amplitude  plane  waves  that  are  obliquely  incident  on  an  ini¬ 
tially  plane  interface  between  two  lossless,  immiscible  fluids.  First,  the  bound¬ 
ary  conditions  at  the  interface  are  examined.  The  interface  is  assumed  to  be 
initially  planar  and  coincident  with  the  z  —  0  plane.  Moreover,  the  effects  of 
surface  tension  and  body  forces  at  the  interface  are  neglected.  Second-order 
perturbation  analysis,  which  is  sometimes  referred  to  as  quasilinear  analysis, 
is  employed  to  analyze  reflection  and  refraction  of  an  abitrary  incident  signal. 
The  notation  closely  follows  the  work  of  Naze  Tjotta  and  Tjotta  (19S7).  The 
0(f)  and  Off2)  systems,  where  e  is  a  small  parameter,  are  then  solved.  The 
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0(e2)  system  accounts  for  not  only  the  nonlinearity  in  the  wave  equation  but 
for  the  finite  displacement  of  the  interface  and  the  variation  of  the  normal  to 
the  interface.  After  the  expressions  for  the  0(e)  and  0(e2)  reflected  and  trans¬ 
mitted  fields  are  obtained,  it  is  noted  that,  to  0(e2),  Snell’s  law  and  the  law  of 
specular  reflection  hold. 

In  the  third  part  of  the  dissertation,  we  develop  two  different  ‘modified 
forms’  of  Snell’s  law — forms  which  appear  to  indicate  that  refraction  has  a 
slight  amplitude  dependence.  The  ‘modified  forms’  are  developed  by  means 
other  than  perturbations  and  appear  to  be  correct  to  second-order.  However, 
in  deriving  the  ‘modified  forms’,  one  of  which  was  reported  by  Cotaras  and 
Blackstock  (1987),  a  large  assumption  is  made:  that  simple  wave  flow  exists 
on  the  incident  side  of  the  interface.  This  assumption  is  quantified  by  using 
the  results  of  the  previous  section  to  develop  the  0(e)  and  0(e2)  conditions 
for  simple  wave  flow.  We  next  expand  each  of  the  ‘modified  forms’  of  Snell’s 
law  and  find  that,  to  0(e2),  the  two  ‘modified  forms’  are  equivalent.  Moreover, 
the  ‘modified  forms’  yield,  at  0(e) ,  ordinary  Snell’s  law  and,  at  0(e2),  one  of 
the  0(e2)  conditions  for  simple  wave  flow.  Since  simple  flow  is  assumed  in  the 
derivation,  the  ‘modified  forms’  are,  to  0(e2),  equivalent  to  ordinary  Snell’s 
law.  Thus  we  metaphorically  deflate  the  balloon  that  carried  us  through  this 
dissertation:  the  hope  of  finding  a  deviation  from  Snell’s  law  at  0(e2). 


CHAPTER  2 


ON  THE  RANKING  OF  TERMS  IN  THE  BASIC  EQUATIONS 
FOR  A  HOMOGENOUS,  THERMOVISCOUS  FLUID  WITH  A 
SINGLE  RELAXATION  MECHANISM 


Discussed  in  this  chapter  are  the  basic  equations  for  a  homogeneous, 
thermoviscous  fluid  with  a  single  relaxation  mechanism  and  a  method  for  rank¬ 
ing  the  terms  in  the  equations.  The  equations  presented  are  the  conservation 
of  mass  (continuity  equation),  the  conservation  of  linear  momentum  (Newton’s 
second  law),  and  the  conservation  of  entropy.  To  the  entropy  equation,  however, 
we  add  the  second  law  of  thermodynamics:  entropy  production  is  positive  for 
irreversible  processes  (loss  mechanisms).  Other  relations  from  thermodynamics 
that  are  presented  are  an  equation  of  state,  Gibb’s  equation,  and  the  first  law 
of  thermodynamics:  the  conservation  of  energy.  The  energy  equation  is  manip¬ 
ulated  into  the  form  of  the  entropy  equation,  and  entropy  production  is  seen  to 
be  the  sum  of  the  products  of  the  various  thermodynamic  fluxes  and  forces.  A 
fundamental  assumption  is  then  made.  The  magnitude  of  the  deviation  from 
thermodynamic  equilibrium  is  assumed  to  be  small,  and  linear  relationships  be¬ 
tween  each  thermodynamic  flux  and  all  the  thermodynamic  forces  are  assumed  to 
hold.  The  possibility  of  cross-effects  between  the  different  thermodynamic  fluxes 
and  forces  is  pointed  out,  but  neglected.  The  fluid  is  assumed  to  be  initially  quiet 
and  in  thermodynamic  equilibrium.  Moreover,  the  signal  is  assumed  to  be  far 
from  all  boundaries  for  all  time.  Additional  fundamental  assumptions  are  then 
made  about  (1)  the  amplitude  of  the  acoustic  signal,  (2)  the  magnitude  of  the 
transport  coefficients,  and  (3)  the  magnitude  of  the  dispersion  caused  by  the  re¬ 
laxation.  The  terms  in  the  equations  are  then  ranked  according  to  their  relative 
importance.  Use  of  the  ranking  system  enables  us  to  readily  obtain  simplified 
forms  of  the  basic  equations  (Lighthill  1956).  Specifically,  wc  develop  simplified 
forms  of  the  equations  appropriate  for  (1)  small-signal  wave  motion  in  a  lossless 
fluid,  (2)  small-signal  wave  motion  in  a  homogeneous,  thermoviscous,  relaxing 
fluid,  and  (3)  finite-amplitude  wave  motion  in  a  homogeneous,  thermoviscous, 
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relaxing  fluid.  As  mentioned  in  Chapter  1,  it  is  our  intention  to  neglect  all  losses 
iii  the  analysis  presented  in  later  chapters.  The  terms  that  represent  the  effects 
of  viscosity,  heat  conduction,  and  relaxation  are,  however,  included  at  th's  stage 
for  completeness. 


2-1  Basic  Equations  for  a  Homogeneous,  Thermoviscous  Fluid  with 
a  Single  Relaxation  Mechanism 


We  start  with  two  classical  equations:  the  conservation  of  mass  and  the 
conservation  of  linear  momentum, 


^  +  V-(pu)  =  0 


(2.1) 


Duk  _  d<7jk 
^  Dt  dxj 


(2.2) 


where  p  is  the  density,  u*  is  the  Hh  component  of  the  particle  velocity  u  (k  = 
1,2,3),  a is  the  stress  tensor,  and  DjDt  is  the  material  derivative, 


Dt  dt  '  dxi 


(2.3) 


Throughout  this  work,  we  follow  the  summation  convention;  that  is,  a  repeated 
index  implies  a  summation  over  the  index.  Equations  (2.1)  and  (2.2)  are  from 
Thompson  [1984,  Eqs.  (1.54)  and  (1.59)],  but  the  effect  of  body  forces  such  as 
gravity  is  neglected  in  the  linear  momentum  equation. 

To  the  pair  of  classical  equations,  we  add  some  thermodynamic  relations: 
a  state  equation,  Gibb’s  equation,  and  the  entropy  equation.  The  state  equation 
gives  the  functional  relationship  between  the  thermodynamic  variables  when  the 
fluid  is  in  thermodynamic  equilibrium.  (Thermodynamic  equilibrium  is  precisely 
defined  below.)  For  a  thermoviscous  fluid  with  a  single  relaxation  mechanism, 
one  form  of  the  state  equation  is1 


P=P(p,s,q)  ,  (2.4) 

'Any  three  thermodynamic  variables  may  be  chosen  as  the  three  independent  variables.  The 
advantage  of  choosing  the  advancement  variable  q  and  the  entropy  s  as  two  of  the  three  is  that, 
if  a  fluid  is  always  in  thermodynamic  equilibrium  (very  common  assumption  in  acoustics),  both 
q  and  s  take  on  their  equilibrium  values,  which  are  constants.  Thus,  the  state  equation  reduces 
to  P  —  I‘(p). 
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where  P  is  the  total  pressure,  s  is  the  entropy  per  unit  mass,  and  q  is  the  degree 
of  advancement  of  the  relaxation  mechanism,  for  example,  a  chemical  reaction. 
It  turns  out  that  the  following  alternate  state  equation  is  required: 

T  =  T(p,  s,  q)  ,  (2.5) 


where  T  is  the  absolute  temperature.  Another  relation  between  the  thermo¬ 
dynamic  variables  is  given  by  Gibb’s  equation,  which  for  a  fluid  with  a  single 
relaxation  mechanism  is 


T ds  =  de - -  dp  —  Adq 


(2.6) 


where  e  is  the  internal  energy  per  unit  mass  and  A  is  the  affinity  of  the  relaxation 
process.  The  entropy  equation  may  be  stated  in  such  a  way  that  it  incorporates 
the  second  law  of  thermodynamics: 


Ds  dJk 
P~Dt  +  dTk 


=  T 


T  >  0 


(2.7) 


where  Jk  is  the  fcth  component  of  the  entropy  flux  vector  J  and  T  is  the  entropy 
source  strength.2 

The  entropy  source  strength  is  an  important  variable  for  two  reasons. 
First,  the  fluid  is  defined  as  being  in  a  state  of  thermodynamic  equilibrium  when 
no  entropy  is  produced,  that  is,  when  T  =  0.  Second,  linear  phenomenological 
relations  may  be  introduced  by  way  of  the  entropy  source  strength.  The  entropy 
source  strength  may  be  written  as  the  sum  of  the  product  of  thermodynamic 
fluxes  and  forces, 

T  =  JkXk  ,  (2.8) 

where  Jk  is  a  generalized  thermodynamic  flux  and  Xk  is  a  generalized  thermo¬ 
dynamic  force.  Since  thermodynamic  equilibrium  is  defined  as  zero  entropy  pro¬ 
duction,  it  is  required  that  all  thermodynamic  fluxes  and  forces  be  zero  at  equi¬ 
librium, 

Jk  =  0  and  Xk  —  0  at  equilibrium  (2.9) 


If  the  departures  from  equilibrium  are  small ,  linear  phenomenological  relations 
between  the  thermodynamic  fluxes  and  forces  may  be  assumed, 


Jk  =  LikXi  ,  (2.10) 

•The  material  in  this  section  is  based  on  the  material  in  Chaps.  Ilf,  IV',  and  V  in  the  hook 
by  I’rigogine  (1%1)  and  Chaps.  Ill,  IV,  and  VI  in  the  hook  by  de  Groot  and  Mazur  (1981). 
References  to  other  works  are  made  when  appropriate. 
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where  L,*  are  called  the  phenomenological  coefficients.  Note  that  although  a  lin¬ 
ear  relationship  between  the  fluxes  and  forces  is  assumed,  the  coefficients  are  not 
assumed  to  be  constants,  and  the  possibility  of  interference  between  the  differ¬ 
ent  fluxes  and  forces  is  permitted.  The  interference  between  the  thermodynamic 
fluxes  and  forces  is,  however,  restricted  by  Curie’s  principle,  which,  briefly  stated, 
is  that  for  an  isotropic  fluid  only  fluxes  and  forces  of  similar  tensorial  nature  may 
interact.  Thus,  a  scalar  flux  may  contain  contributions  from  scalar  forces,  but 
not  from  vector  forces  or  tensor  forces.  Moreover,  the  number  of  coefficients  is 
limited  by  Onsager’s  relation,  which  states  that  L,k  = 


Finding  the  appropriate  form  of  T 


To  find  the  form  of  T  appropriate  for  a  thermoviscous  fluid  with  a  single 
relaxation  mechanism,  we  start  with  the  energy  equation  and  rearrange  it  so  that 
it  is  in  the  form  of  the  entropy  equation.  In  rearranging  the  equation,  we  assume 
local  equilibrium  and  split  the  stress  tensor  into  two  parts,  a  pressure  term  and 
a  viscous  stress  term. 

We  start  with  the  energy  equation;  see,  for  example,  Thompson  [1984, 
Eq.  (1.60)], 


D  (  u2\ 

"Di  e  +  T/  = 


0  .  dQ, 

diyikui)  - 


(2.11) 


where  Q,  is  the  tth  component  of  the  heat  flux  vector  Q.  If  the  kinetic  energy 
relation  (which  is  obtained  by  multiplying  the  linear  momentum  equation  by  u*.) 
is  subtracted  from  the  energy  equation,  the  resulting  equation  is  [see,  for  example, 
Thompson  (1984,  Eq.  1.70)] 


De  _  dui  dQi 

^  Dt  'kdxk  dxi 

The  stress  tensor  o,k  may  be  divided  into  two  parts: 


(2.12) 


&ik  =  ~Pf>ik  ~  v(k  > 


(2.13) 


where  cr'ik  is  the  viscous  stress  tensor,  —P6,k  is  the  pressure  stress,  and  6,k  is  the 
Kroneckcr  delta  function,  which  is  equal  to  unity  if  i  —  k,  but  zero  otherwise.3 

■’The  choice  of  sign  on  the  viscous  stress  tensor  in  Eq.  (2.13),  which  is  arbitrary,  differs 
from  some  work  in  hydrodynamics  [see,  for  example,  Landau  and  Lifshitz  1959,  Eq.  (15.2)], 
but  conforms  to  some  work  in  thermodynamics  [see,  for  example,  de  Groot  and  Mazur  1984, 
Eq.  (Ill  35)]. 
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The  viscous  stress  tensor  is  assumed  to  be  symmetric.4  Using  Eq.  (2.13)  in 
Eq.  (2.12)  and  using  the  continuity  equation,  Eq.  (2.1),  leads  to  the  following 
form  of  the  energy  equation: 

Dt  p2  Dt  ~  pG'kdxk  pdxk  ' 

Equation  (2.14)  may  be  simplified  using  Gibb’s  equation.  However, 
Gibb’s  relation  is,  strictly  speaking,  only  valid  in  equilibrium.  Since  the  fluid 
is  not  always  in  equilibrium,  we  must  assume  that  the  fluid  particles  are  in  a 
state  of  local  equilibrium.  If  we  follow  a  fluid  particle  that  is  in  local  equilibrium, 
Gibb’s  equation  may  be  written  as 


Ds  De  P  Dp  Dq 

Dt  ~  Dt  p2  Dt  Dt  ' 

Use  of  Eq.  (2.15)  in  Eq.  (2.14)  yields 

Ds  1  ,  dui  1  dQi  pA  Dq 

p~Dt  =  ~f<7ikdTk  ~  fdTk  ~  Tm 


(2.15) 


(2.16) 


To  place  Eq.  (2.16)  in  the  form  of  the  entropy  equation,  Eq.  (2.7),  wc 
must  perform  some  manipulations  that  initially  appear  arbitrary.  The  resulting 
expressions  are,  however,  unique;  see  de  Groot  and  Mazur  (1984,  pp.  24-25).  The 
simpler  of  the  manipulations  is 


1  dQi  d  (Qt\  Qi  dT 
T  dx,  dxt\Tj  T2  dx. 


(2.17) 


The  more  complex  manipulation  involves  expanding  a[k  and  dut/dxk  into  trace 
components  and  remainders  that  have  zero  trace.  We  rewrite  the  viscous  stress 
tensor  as  follows: 


a',k  =  °\k  +  ¥*°'n 


(2.18) 


where  a\k  is  a  symmetric  tensor  since  <k  is  assumed  symmetric.  The  gradient  of 
the  velocity  is  rewritten  as 


On  i 
d.vk 


1 

(  0u{ 

.  0uk\ 

2  x  du( 

2 

\dxk 

+  aJ~j 

1  (  dux 

+  2  la^T 


duk  \  1  du( 
dx,J  3  ‘k  dx t 


4 This  is  tantamount  to  assuming  that  the  internal  angular  momentum  is  randomly  dis¬ 
tributed  amongst  the  molecules,  which  is  generally  the  case.  If  this  is  not  assumed,  then  it 
turns  out  that  another  phenomenological  coefficient,  in  this  case,  the  rotational  viscosity,  is  re- 
i|uired.  Moreover,  the  asymmetry  of  the  viscous  stress  tensor  decays;  see  de  Groot  and  Mazur 
(l!)81,  Chap  XII). 
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The  first  term  on  the  right-hand  side,  which  is  symmetric  and  has  zero  trace, 
indicates  shearing  motion  but  without  a  volume  change.  The  second  term,  which 
is  antisymmetric,  is  associated  with  a  rotation  but  again  without  a  volume  change. 
The  third  term  on  the  right-hand  side,  which  consists  solely  of  diagonal  elements, 
indicates  the  volume  change.  Since  the  product  of  symmetric  and  antisymmetric 
tensors  is  zero,  we  may  write 


dui 

'“dTk 


<7., 


'  du;  duk 
dxk  +  dxi 


- 


du, 

dx. 


,  i  / 

+  3  13  dxt 


(2.19) 


Use  of  Eqs.  (2.17)  and  (2.19)  in  Eq.  (2.16)  yields  the  desired  form  of  the  energy 
equation: 


Ds  d 
^  Dt  dx, 


Qi 

r 


QLdT__K!L 
T2  dxi  2  T 

1  (7jj  du, 


dui  duk 
dxk  dx,  t 
pA  Dq 


2  c  du( 

-  '^07, 


3  T  dx,  T  DT 


(2.20) 


Comparing  Eq.  (2.20)  with  the  entropy  equation,  Eq.  (2.7),  we  identify  the  fol¬ 
lowing: 

Ji  =  Y  (2-21) 

and 


T  =  - 


<2,  dT 


aik 

2  T 


'  du ,  duk 
dxk  dx 


k\  2c  du( 

I)  ” 


1  ajj  dut  _  pADq^ 

3  T  dx,  T  Dt  '  ( 


T2  dx, 

Introduction  of  the  linear  phenomenological  relations  and  constants 

We  are  now  in  a  position  to  utilize  Eqs.  (2.8)  and  (2.10),  the  linear 
phenomenological  equations.  Recalling  Curie’s  principle  and  Onsager’s  relations, 
we  write 

Qi  =  .  (2-23) 


T2  dx, 


l2 

( — i  i  dxik' 

\  -S  j. 

du. 

2  T 

\dxk  dxt/ 

)  ~  3°>k 

dx. 

1  t 

L n 

du. 

,  A 

T 

dx. 

•> 

Dq 

T  21 

du. 

,  A 

pm 

T 

dx. 

■  ^22  Y 

1 

(2.24) 

(2.25) 


(2.26) 
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where 


L\2  ~  L 


21 


(2-27) 


The  coefficients  Lx ,  £2,  T11 ,  and  Z-2 2  must  be  positive  in  order  to  satisfy  the  second 
law  of  thermodynamics,  T  >  0.  For  simplicity,  we  shall  neglect  the  possiblity 
of  cross-effects  throughout  the  remainder  of  the  work,  that  is,  we  assume  that 
T12  =  L  21  =  0. 

We  now  identify  some  of  the  phenomenological  coefficients  as  the  trans¬ 
port  coefficients  in  well-known  linear  relations  (for  example,  Fourier’s  law  of  heat 
conduction): 


..  _  Lx 

K  ’ 

(2.28) 

_  L* 

2T  ’ 

(2.29) 

II 

CD 

=3. 

(2.30) 

where  n  is  the  thermal  conductivity,  p  is  the  shear  viscosity,  and  fi#  is  the  bulk 
viscosity.  Use  of  Eqs.  (2.28)-(2.30)  in  Eqs.  (2.23)-(2.25)  yields  the  following  well- 
known  linear  phenomenological  laws  (recall  that  the  cross-effect  in  Eq.  (2.25)  is 
neglected): 

Q  =  -kVT  ,  (2.31) 


/  dui 

2  c  du' 

\dxk 

dXi) 

- 

H;  =  Mb(V-u) 


(2.32) 

(2.33) 


Obtaining  the  fourth  phenomenological  constant  and  relation  takes  a 
little  more  work.  First,  for  future  reference,  we  redefine  the  phenomenological 
coefficient  as 

0  =  ■  (2-3-1) 

If  cross-effects  are  neglected,  this  definition  results  in  the  following  linear  phe¬ 
nomenological  relation: 

wt  =  ■  <235) 


It  turns  out.  however,  that  different  forms  of  the  linear  phenomenological  relation 
and  its  associated  coefficient,  which  are  referred  to  as  the  rate  equation  and  the 
relaxation  lime,  arc  of  great  use.  The  rate  equation  and  the  relaxation  time  are 
obtained  as  follows:  The  affinity  A  =  A(p,T,q)  is  expanded  in  a  Taylor  series. 
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The  expansion  is  then  specialized  for  equilibrium,  and  the  result  is  subtracted 
from  the  original  expansion.  The  results  follow  directly.  In  general,  for  small 
deviations  from  equilibrium,  A  may  be  approximated  as 


A  = 


(P  ~  Po )  + 
o 


(T  -  T0)  + 


(<7  -<? o) 
o 


(2.36) 


where  the  subscript  0  denotes  the  static  value  of  the  variable.  (This  notation 
is  used  throughout  this  work.)  As  noted  earlier,  all  thermodynamic  forces 
and  fluxes  are  zero  at  equilibrium;  thus,  the  equilibrium- — static — value  of  A 
is  zero.  The  equilibrium  value  of  the  advancement  variable  q  is  denoted  by  a 
superscript  *. 


*U-»-  £ 


T.1J 


,  ,  /  dA 

\P-P0)+  I 


(T-T0)  + 


P’l/  0 


dA 

dq 


(?*-9o) 


pJ) 


(2.37) 

Equation  (2.37)  indicates  that  q *  is  a  function  of  only  two  thermodynamic  vari¬ 
ables,  p  and  T.  Note  that  the  static — equilibrium— values  of  q  and  q'  are  equal, 
that  is, 

<?o  =  ?o  • 

Subtracting  Eq.  (2.37)  from  Eq.  (2.36)  yields 


A  = 


(?-?*) 


(2.38) 


This  is  consistent  with  our  earlier  assumption  that  the  thermodynamic  fluxes 
and  forces  be  linearly  related.  Use  of  Eq.  (2.38)  in  the  linear  phenomenological 
relation  (again  neglecting  cross-effects),  Eq.  (2.26),  leads  to 


£q  _  {q  -  <?*) 

Dt  T0 


(2.39) 


where 


t0  =  - 


pT 


^(wL)0 


(2.40) 


The  value  of  r0  is  positive  because  L2 2  must  be  positive  to  satisfy  the  second  law, 
T  >  0,  and  because  ^§~|  is  positive;  see  de  Groot  and  Mazur  (1984,  p.  201) 
for  a  complete  explanation. 
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Introduction  of  the  linear  phenomenological  relations  into  the  momen¬ 
tum  and  entropy  equations 

The  terms  in  the  entropy  equation,  Eq.  (2.7),  that  were  unknown  are  the 
entropy  flux,  J,  and  the  entropy  source  strength,  T.  The  entropy  flux  is  given  in 
Eq.  (2.21).  The  entropy  source  strength  may  be  obtained  if  the  linear  relations 
between  the  thermodynamic  forces  and  fluxes  given  in  Eqs.  (2.31),  (2.32),  (2.33), 
and  (2.35)  are  substituted  into  the  relation  for  T,  Eq.  (2.22), 


k  l  OT 


H  [du,  duk 


T  —  _  1  '  '  »  K  __  2  c 


T2  \dxij  2 T  dxfc  dx. 


/iB  ( du( 


T  \dx( 


~  )  +—A1  ■  (2.41) 


Note  that  all  the  terms  are  positive  and  the  squares  of  variables. 

The  linear  phenomenological  relations  arc  now  introduced  into  the  mo¬ 
mentum  equation.  Recall  that  the  viscous  stress  tensor  was  split  into  two  parts: 
a  trace  component  and  a  symmetric  tensor  with  zero  trace.  Use  of  Eqs.  (2.32) 
and  (2.33)  in  Eq.  (2.18)  yields 


(  du,  duk\  2  dut  1 


(2.42) 


Use  of  Eqs.  (2.13)  and  (2.42)  in  Eq.  (2.2)  yields  the  following  form  of  the  linear 
momentum  equation: 

Du,  dP  ,  du,  ,  (  ,  p\  d  due 


Du i  _  dP  du, 

^  Dt  dii  ^  ^  dxkdxk 


+  (m  + 1)  d7,^r,  •  {2M) 


which  may  also  be  written  in  vector  form  as 

PJX  =  -VP  +  /tV2u  +  (/iB  +  V( V-u)  .  (2,14) 

In  obtaining  Eq.  (2.43),  we  assumed  that  the  transport  coefficients  /t,  /i,  and 
fi-B  arc  constants.5  Equation  (2.44)  is  the  well-known  form  of  the  momentum 
equation;  see  Landau  and  Lifshitz  (1959,  Eq.  (15.6)]. 


’This  assumption  is,  strictly  speaking,  not  necessary;  it  is,  however,  convenient.  It  is  not 
necessary  her  a  use  it  turns  out  that  the  ranking  scheme,  which  is  introduced  later  in  the  chapter, 
indicates  that  the  terms  which  arise  if  the  transport  coefficients  arc  not  assumed  constant  are 
negligible.  As  an  example,  we  now  examine  one  of  the  terms  containing  a  transport  coefficient, 
The  gradient  of  this  term  is 

V(/<|(V-u)  =  /(DV2u  +  (V-u)V/t|(  . 

Hut  since  /m  is  a  function  of  any  three  stale  variables,  say  P,  s,  and  q ,  we  may  expand  /in  in 
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Summary  of  basic  equations  for  a  homogeneous,  thermoviscous  fluid 
with  a  single  relaxation  mechanism 

The  basic  equations  for  a  homogeneous,  thermoviscous  fluid  with  a  single 
relaxation  mechanism  are  the  continuity  equation,  the  linear  momentum  equa¬ 
tion,  the  entropy  equation,  the  equation  of  state,  the  alternate  equation  of  state, 
and  the  rate  equation.  The  equations  are,  respectively, 

^  +  V.(pu)  =  0  ,  (2.1) 

^a.VP  +  ,y»u+(W  +  f)v(V,)  ,  (2.44) 

Ds  _  [VT\  k  (dT\2  n 
P  Dt  K  \  T  )  ~  T2  \&cj  +  2 T 

+  ^ 

T 

P=P(p,s,q)  , 

T  =  T(p,  s,  q)  , 

£l  _  (*?  ~  ?") 

Dt  T0 


OUl  OUk  _  2  OU £ 
dxk  dx,  ^  lk  dx( 

'du'\\P'P  A2 
+  ~tA  ’ 


(2.45) 

(2.4) 

(2.5) 


(2.39) 


2-2  The  Ranking  System  and  Its  Fundamental  Assumptions 

In  this  section  the  thermodynamic  variables  are  defined  to  be  the  sum  of 
a  static  value  and  a  small  fluctuation.  Fundamental  assumptions  are  then  made 
about  the  magnitude  of  the  acoustic  fluctuation,  the  magnitude  of  the  transport 
coefficients,  and  the  amount  of  dispersion  caused  by  the  relaxation  mechanism. 
The  system  for  ranking  the  terms  in  the  basic  equations  is  then  presented. 

a  Taylor  series, 


/'b  =  (/'b)o  + 


(r-P0)  + 


(s  -  So)  + 

0 


(?-?o)+  ••• 

0 


Use  of  the  ranking  scheme  indicates  that  the  only  term  in  that  is  not  negligible  is 

(/in)0V-u.  Thus  the  ranking  system  will  show  us  that  the  variation  in  the  transport  coefficients 
is  a  negligible  effect,,  but  it  is  simpler  for  us  to  treat  them  as  constants  from  this  point  on. 
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We  may  define  each  thermodynamic  variable  to  be  the  sum  of  a  static 
value,  which  is  denoted  by  the  subscript  0,  and  a  small  fluctuation,  which  is 
denoted  by  the  superscript  /, 


p  =  p0  +  p  , 

(2,16) 

P  =  Po  +  P  , 

(2.-17) 

T  =  To  +  r  , 

(2,18) 

5  =  So  +  s'  , 

(2.19) 

u  =  u0  +  u'  , 

(2,50) 

q  =  q0  +  q  , 

(2,51) 

*  _  *  +! 

<!=%  +  <! 

(2,52) 

Note  that  because  the  fluid  was  assumed  to  be  initially  quiet,  the  static  value  of 
the  particle  velocity  is  zero,  that  is,  |u0|  =  0.  We  therefore  drop  the  unnecessary 
prime  from  u'  and  refer  to  the  acoustic  fluctuation  in  the  particle  velocity  sim¬ 
ply  as  u.  Recall  that  the  static  values  of  q  and  q ’  are  the  same  (i.e.,  </0  =  </*,). 
liquations  (2.46)  through  (2.52)  may  be  substituted  directly  in  the  continuity, 
momentum,  entropy,  and  rate  equations,  but  this  is  deferred  until  after  the  rank¬ 
ing  system  is  presented. 

It  is  well  known  that  for  most  fluids  the  linear,  lossless,  wave  equation 
(the  classical  wave  equation)  is  a  good  model  for  describing  the  behaviour  of 
small-signal  acoustic  waves.  The  classical  wave  equation  may  be  obtained  by 
retaining  only  linear,  lossless  terms  in  the  basic  equations.  A  linear,  lossless  term 
is  a  linear  term  with  a  coefficient  that  does  not  involve  a  transport  coefficient  or 
the  derivative  of  the  equilibrium  variable  q*  with  respect  to  another  thermody¬ 
namic  variable.  Examples  are  dp'/di  and  p0(^#u).  The  linear,  lossless  terms  an' 
therefore  ranked  as  first-order  terms. 

On  the  other  hand,  if  one  wants  a  better  approximation  an  approxima¬ 
tion  that  accounts  for  the  weak  effects  of  viscosity,  heat  conduction,  relaxation, 
and  nonlinearity  then  more  terms  must  be  retained.  The  terms  to  be  retained 
are  linear  relaxation  terms,  linear  dissipation  terms,  and  quadratic  nonlinear¬ 
ity  terms.  'These  terms,  which  are  defined  below,  are  ranked  as  second-order 
terms  and  represent  the  most  important  effects  of  viscosity,  heat  conduction, 
relaxation,  and  nonlinearity.  A  linear  relaxation  term  is  a  linear  term  with  a 
coefficient  that  involves  a  derivative  of  (/*  but  not  a  transport,  coefficient.  An  ex 
•'tuple  is  />'  (  )  .  A  linear  dissipation  term  is  a  linear  term  with  a  coefficient 
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that  contains  a  transport  coefficient  but  not  a  derivative  of  g*.  For  example, 
pV2  u  is  a  linear  dissipation  term.  Terms  containing  a  quadratic  nonlinear  term 
with  a  coefficient  that  does  not  contain  either  a  transport  coefficient  or  a  deriva¬ 
tive  of  <7*  are  referred  to  as  quadratic  nonlinearity  terms;  examples  are  p'  (du/dt) 
and  p'  (V-u).  The  names  of  the  various  terms  are  reflective  of  their  physical 
interpretation. 

Since  the  effects  of  viscosity,  heat  conduction,  relaxation,  and  nonlinear¬ 
ity  are  assumed  small,  any  term  representing  the  interaction  of  any  two  second- 
order  effects  would  be  expected  to  be  negligible.  Such  terms  are  encompassed  in 
the  third  major  category,  higher-order  terms.  An  example  of  a  higher-order  term 
is  the  cubic  term,  p'( u»V)u. 

To  quantify  the  ideas  in  the  previous  paragraphs,  we  now  state  three 
fundamental  assumptions  about  the  magnitude  of  the  various  terms:  First,  the 
magnitude  of  the  particle  velocity  |uj  is  assumed  to  be  small  with  respect  to  the 
sound  speed  e0,  that  is, 

|u|  <  Co  . 

It  turns  out  that  this  assumption  also  implies  that  the  magnitude  of  the  pressure 
fluctuation  7/  is  small  with  respect  to  poc$,  and  that  the  magnitude  of  the  density 
fluctuation  p  is  small  with  respect  to  p0  (see  Appendix  B), 

I?' I  <  PoCq  , 

Ip'I  <  Po  ■ 

To  aid  in  the  subsequent  ranking,  we  state  that  the  nondimensional  acoustic 
fluctuations  p'/PoUe  p'/Pin  and  u/cq  are  of  order  e,  where  e  is  a  small  parameter. 
Our  second  fundamental  assumption  is  that  for  the  highest  frequency  of  interest 
the  effects  of  shear  viscosity,  bulk  viscosity,  and  heat  conduction  are  small.  This 
means  that,  when  suitably  nondimensionalized  (see  Appendix  B),  the  transport 
coefficients  p,  pg,  and  k  are  all  small  quantities  of  order  u,  where  v  is  a  small 
parameter.  Our  third  fundamental  assumption  is  that  the  dispersion  caused  by 
the  relaxation  is  of  order  Af,  where  M  is  a  small  nondimensional  parameter 
that  is  formally  defined  later.  At  this  point  it  is  appropriate  to  recall  that  two 
fundamental  assumptions  were  previously  made:  the  deviation  from  equilibrium 
is  assumed  to  be  small  enough  that  the  thermodynamic  fluxes  and  forces  may  be 
linearly  related,  and  the  cross-effects  are  assumed  negligible. 

Use  of  the  ranking  system  makes  it  relatively  easy  to  deal  with  special 
cases.  For  the  case  of  small-signal  wave  motion  in  a  lossless  fluid,  only  terms 
of  order  e  need  to  be  retained.  For  finite-amplitude  wave  motion  in  a  lossless 
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fluid,  terms  of  order  £  and  £2  must  be  retained.  For  finite-amplitude  signals  in  a 
homogeneous,  thermoviscous,  relaxing  fluid,  terms  of  order  e,  e2,  Me,  ve  must 
all  be  retained.  Thus,  the  basic  equations  can  be  simplified  and,  in  most  cases, 
combined  to  form  a  wave  equation  corresponding  to  the  situation  at  hand. 

It  is  important  to  note  that  not  all  second-order  nonlinearity  terms  gen¬ 
erate  similar  effects.  The  effects  of  nonlinearity  may  be  roughly  divided  into  two 
catagories:  cumulative  effects,  that  is,  effects  which  grow  with  distance,  and  non- 
cumulative  or  local  effects.  Since  no  a  priori  method  exists  to  indicate  whether 
a  particular  second-order  nonlinearity  term  will  result  in  effects  that  grow  with 
distance,  all  second-order  terms  must  be  retained.  Under  some  circumstances 
local  effects  can  be  important  because,  although  they  do  not  grow  with  distance, 
they  can  propagate.  Thus,  a  local  effect  generated  at  one  location  may  influence 
(or  render  impossible)  the  measurement  of  a  local  effect  at  a  second  location. 

If  Eqs.  (2.46)  through  (2.52)  are  substituted  into  the  continuity,  mo¬ 
mentum,  entropy,  and  rate  equations  as  well  as  the  Taylor  series  expansions  of 
the  state  equations  and  if  the  terms  that  are  of  higher  order  are  neglected,  the 
following  equations  are  obtained: 

^■  +  p0V-a  =  -p'V-u-u-Vp'  ,  (2.53) 


(>0^  l-Vp'=-/,'^-/>o(u-V)u  +  (,v!u+(/IB  +  £)  V(V-u)  ,  (2.54) 

PoT0—  =  kV2T'  (2.55) 

K  +  n.vq'=q- '-Z<1  ,  (2.56) 

Ot  Tq 


?'  =  P 


+?'(¥l  +W)7 


(2.57) 


I  =  p 


W  £ 


■  UJ\2  °  * 

a  _  "h  2  '•P  )  o  2 


.  (2.58) 


It  turns  out,  that  both  the  momentum  equation  and  the  rate  equation  may  be 
further  simplified,  but  this  is  deferred  until  later. 
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2-3  The  Sound  Speed  in  a  Relaxing  Fluid 


The  defining  relation  for  the  sound  speed  is  usually  written  as 


c 


2 


(2.59) 


Since  three  independent  thermodynamic  variables  are  required  to  completely 
specify  the  state  of  a  fluid  with  a  single  relaxation  mechanism  and  since  Eq.  (2.59) 
only  specifies  how  two  of  them  are  to  be  used,  the  sound  speed  for  a  relaxing  fluid 
may  be  defined  in  a  multitude  of  ways.  It  turns  out  (as  we  see  in  the  next  chap¬ 
ter)  that  the  following  two  defining  relations  play  an  important  role  in  acoustic 
wave  propagation  in  a  relaxing  fluid: 


(c~)2  = 


d_P_ 

dp 


Stq 


=  -^P{p,stq) 


(2.60) 


(c0)2  = 


dP_ 

dp 


M=«* 


C\ 


(2.61) 


Equation  (2.60)  is  the  defining  relation  for  the  frozen  sound  speed  c°°,  the  sound 
speed  for  a  fixed  ( frozen )  value  of  q  that  is  not  necessarily  the  equilibrium  value 
q" .  Equation  (2.61)  is  the  defining  relation  for  the  equilibrium  sound  speed  c°,  the 
sound  speed  for  the  equilibrium  value  of  q  —  qm  that  is  not  necessarily  a  constant. 
It  turns  out  (as  we  see  in  the  next  chapter)  that  the  frozen  sound  speed  is  the 
appropriate  sound  speed  for  a  fluid  with  a  very  long  relaxation  time  and/or  for 
a  very  high  frequency  signal.  In  this  case,  the  signal  disturbs  a  fluid  particle  at 
a  rate  that  is  much  faster  than  the  rate  of  the  relaxation  mechanism,  that  is, 
the  fluid  is  effectively  frozen  in  a  single  state.  Thus,  q  is  constant.  On  the  other 
hand,  the  equilibrium  sound  speed  is  the  appropriate  sound  speed  for  a  fluid  with 
a  very  short  relaxation  time  and/or  for  a  very  low  frequency  signal.  In  this  case, 
the  signal  disturbs  a  fluid  particle  at  a  rate  that  is  much  slower  than  the  rate  of 
the  relaxation  mechanism.  The  advancement  variable  q  is  thus  always  equal  to 
its  equilibrium  value  q' ,  even  though  q *  changes  with  the  passing  acoustic  wave.6 

When  evaluated  at  static  conditions,  the  values  of  the  frozen  and  equi¬ 
librium  sound  speeds  are  denoted  with  a  subscript  0, 


(tfr  = 


(2.62) 


f'  For  a  more  det  ailed  cxpl.anation  of  the  roles  of  the  frozen  and  equilibrium  sound  speeds, 
see,  for  example,  Vincent i  and  Kruger  (1965,  Chap.  VIII,  Secs.  3  and  d). 
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/  dP 

(c°o)2=  £- 

dp 


s,q=q’ 


(2.63) 


Magnitude  of  the  difference  between  the  two  sound  speeds 

When  the  ranking  system  was  presented,  we  assumed  that  the  dispersion 
caused  by  the  relaxation  mechanism  was  a  small  quantity  of  order  M. .  An  expres¬ 
sion  for  the  nondimensional  difference  between  the  equilibrium  and  frozen  sound 
speeds,  which  is  denoted  m  and  was  previously  assumed  small,  may  be  obtained 
in  the  following  manner:  The  total  differential  of  the  pressure  is  first  obtained 
and  is  then  specialized  for  equilibrium.  The  difference  between  the  frozen  and 
equilibrium  sound  speeds  follows  directly.  The  total  differential  of  the  pressure  is 


dp 


,  dP 

,  dP 

dp+17 

a,q 

ds  -f  — 
dq 

p>q  * 

dq 


p,a 


(2.64) 


However,  at  equilibrium  the  change  in  the  entropy  is  zero,  and  the  advancement 
variable  q  takes  on  its  equilibrium  value  <7*, 


dp 


dp  + 


S,q 


dP 

dq 


dq" 


(2.65) 


Rearranging  using  the  defining  relations  for  the  equilibrium  and  frozen  sound 
speeds  and  noting  the  following  identity  from  calculus,7 


dP 

dp 


dp 

dq 

dq 

P,  dP 

=  -1 


we  obtain 

where  m  is  given  by 


(c°)2  =  (c°°)2  (1  —  m)  , 

dpi 


m  — 


dq 


P,3 


df 

dp 


‘  For  a  general  proof  of  the  identity 


dx 

dy 

dz 

0^ 

a  dz 

r  dx 

=  -l  , 


(2.66) 

(2.67) 

(2.68) 


see,  for  example,  Van  Wylen  and  Sonntag  (1976,  p.  366). 
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Rearranging  Eq.  (2.67)  yields  a  more  readily  interpreted  definition  for  m, 


(c°°)2  -  (c0)2 


(coo)2 


=  m 


(2.69) 


The  magnitude  of  m  was  assumed  to  be  of  order  M. ,  where  M.  is  a  small  quantity. 
As  with  the  frozen  and  equilibrium  sound  speeds,  m  has  a  static  value  that  is 
denoted  mo. 

As  the  value  of  m  tends  to  0,  the  effects  of  relaxation  become  small.  At 
the  same  time,  the  values  of  q  and  q‘  tend  to  the  same  constant,  equilibrium 
value,  f/0  =  In  this  limiting  case,  the  equilibrium  and  frozen  sound  speeds  are 
equal,  and,  since  q  is  no  longer  a  variable,  we  recover  the  usual  definition  of  the 
sound  speed  as  given  in  Eq.  (2.59), 


=  (c°°)2  =  (c0)2  .  (2.70) 

We  denote  the  limiting  case  of  m  — ►  0  in  our  notation  by  dropping  the  specifica¬ 
tion  on  q  from  the  partial  differential.  This  notation  is  also  used  to  denote  when 
it  is  inappropriate  to  distinguish  between  the  equilibrium  and  frozen  forms  of  a 
thermodynamic  fluid  property;  see  below.  When  evaluated  at  static  conditions, 
the  sound  speed  for  the  m  =  0  case  is  referred  to  as  the  small-signal  sound  speed 
and  is  denoted  Cq, 


To  show  that  the  difference  between  the  equilibrium  and  frozen  forms 
of  a  first-order  term  is  a  second-order  term,  we  may  rearrange  the  definition  of 
m  and  then  multiply  through  by  a  first-order  term,  say,  the  fluctuation  in  the 
density, 

(tfW  =  (cS)V  +  «*>(<$)V  • 

Since  the  term  m0(c°)2//  is  of  second  order,  it  is  not  appropriate  to  distinguish 
between  the  equilibrium  and  frozen  forms  of  the  sound  speed.  (A  higher-order 
term  would  be  introduced.)  The  above  expression  is  thus  written  correct  to 
second  order  as 

(^)V  =  (^)V  +  m0(co)V  •  (2-72) 

We  may  generalize  the  foregoing  discussion.  Since  the  difference  between 
the  equilibrium  and  frozen  forms  of  a  thermodynamic  fluid  property  is  a  second- 
order  term,  it  is  appropriate  to  distinguish  between  the  equilibrium  and  frozen 


m 


c2  = 


dP_ 

dp 


forms  of  a  thermodynamic  fluid  property  only  when  the  property  is  a  coefficient 
of  a  first-order  term.  This  is  because  the  difference  between  the  equilibrium 
and  frozen  forms  of  a  first-order  term  is  a  second-order  term  (of  order  Aft), 
whereas  the  difference  between  the  equilibrium  and  frozen  forms  of  a  second- 
order  term  is  a  higher-order  term.  In  second-order  terms,  we  therefore  neglect 
the  effects  of  relaxation.  This  is  indicated  notationally  by  not  specifying  what  is 
to  be  done  with  q  in  the  partial  derivatives  that  are  the  defining  relations  for  the 
thermodynamic  fluid  properties.8 

2-4  First-Order  Acoustics:  Small  Signals  in  a  Lossless  Fluid 

In  this  section  the  case  of  small-signal  wave  motion  in  a  lossless  fluid 
is  considered.  The  terms  in  the  basic  equations  that  pertain  to  nonlinearity, 
viscosity,  heat  conduction,  and  relaxation  are  dropped,  and  it  is  then  seen  that 
the  basic  equations  for  this  case  are  composed  only  of  first-order  terms.  It  is  also 
noted  that  neglecting  relaxation  and  heat  conduction  eliminates  the  need  for  s 
and  q  as  variables.  The  motivation  for  developing  some  first-order  relations  at 
this  time  is  to  assist  us  in  simplifying  the  second-order  relations  later.  First-order 
relations  are  of  assistance  because  the  dependent  variables  in  second-order  terms 
may  be  replaced  using  first-order  relations  without  changing  the  overall  level  of 
approximation;  see,  for  example,  Lighthill  (1956). 

If  all  second-order  terms  are  neglected,  the  following  first-order  forms  of 
the  continuity  and  momentum  equations,  Eqs.  (2.53)  and  (2.54),  are  obtained: 


±  +  PoV-u  =  0  ’ 

(2.73) 

du  , 

A>W  +  Vp=°  . 

(2.74) 

Equations  (2.73)  and  (2.74)  are  referred  to  as  the  linear  continuity  and  momentum 
equations. 

We  now  discuss  the  rate  and  entropy  equations  or,  more  accurately,  the 
lack  of  need  for  them.  Since  the  fluid  is  assumed  lossless,  the  fluid  is  assumed 

8If  the  difference  between  the  frozen  and  equilibrium  forms  of  second-order  terms  were  not 
negligible,  that  is  to  say,  if  m  were  not  small,  then  frozen  and  equilibrium  forms  of  both  the  loss 
coefficient  (defined  below)  and  the  parameter  of  nonlinearity  (also  defined  below)  would  emerge. 
However,  in  this  work,  m  is  assumed  to  be  a  small  parameter,  and  the  differences  between  the 
equilibrium  and  frozen  forms  of  terms  containing  the  loss  coefficient  and  the  parameter  of 
nonlinearity  are  of  higher  order. 
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to  be  in  a  state  of  thermodynamic  equilibrium  at  all  times.  The  variation  of  s 
and  q  for  any  fluid  particle  is  thus  zero.  Moreover,  because  the  fluid  is  assumed 
homogeneous,  both  s  and  q  are  constant  throughout  the  fluid.  The  rate  and 
entropy  equations  thus  reduce  to 

q  =  q*  =  constant  , 


s  =  constant  . 

Because  s  and  q  are  constant  throughout  the  fluid,  the  state  equation, 
Eq.  (2.4),  reduces  to 

P  =  P(p,  s  =  constant,  q  =  constant)  .  (2.75) 

Equation  (2.75)  is  sometimes  referred  to  as  the  isentropic  equation  of  state.  The 
linearized  Taylor  series  expansion  of  Eq.  (2.75)  is 


However,  as  noted  above,  it  is  inappropriate  to  distinguish  between  the  equi¬ 
librium  and  frozen  forms  of  the  sound  speed  when  relaxation  is  not  important. 
The  linearized  Taylor  series  expansion  of  the  state  equation  for  a  lossless  fluid  is 
therefore  written  as 


/ 

P 


2  ! 
=  CoP 


(2.76) 


Equations  (2.73),  (2.74),  and  (2.76)  are  the  first-order  forms  of  the  conti¬ 
nuity,  momentum,  and  state  equations  and  are  valid  for  small-signal  wave  motion 
in  a  lossless  fluid.  These  equations  may  be  combined  to  form  wave  equations  in 
p'  and  p\ 


(2.77) 


1  d2p' 
Cq  dt 2 


(2.78) 


2-5  Second-Order  Forms  of  the  Basic  Equations 

In  this  section  a  consistent  second-order  form  of  the  momentum  equation 
is  obtained.  The  order  £2  terms  in  both  the  continuity  and  momentum  equations 
are  then  expressed  in  terms  of  a  new  variable  called  the  Lagrangian  density  C. 
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The  Lagrangian  density  is  related  to  some  of  the  local  finite-amplitude  effects. 
Next,  the  second-order  forms  of  the  entropy,  state,  and  rate  equations  are  ob¬ 
tained.  These  equations  are  then  combined  into  a  single  equation  that  is  referred 
to  as  the  modified  state  equation.  Last,  the  integral  of  the  modified  state  equation 
is  obtained. 

Continuity  and  momentum  equations 

No  terms  in  the  expanded  continuity  equation,  Eq.  (2.53),  are  of  higher 
order.  Thus,  the  second-order  form  of  the  continuity  equation  is  exact.  Note 
that  the  second-order  terms  in  the  continuity  equation  are  quadratic  nonlinearity 
terms  and  need  only  be  retained  when  finite- amplitude  effects  are  of  interest. 

A  consistent  second-order  approximation  of  the  expanded  momentum 
equation,  Eq.  (2.54),  is  now  obtained  for  a  signal  propagating  far  from  boundaries 
in  an  initially  irrotational  fluid.  We  first  introduce  the  vorticity  f?, 

f?  =  Vxu  .  (2.79) 

The  terms  (u-V)u  and  V(V*u)  in  the  momentum  equation  may  be  expanded 
using  vector  identities  [Gradshteyn  and  Ryzhik  1980,  Eq.  (10.31.3)]  and  the  def¬ 
inition  of  the  vorticity  to  give 

(u-V)u  =  |V(u2)  +  ttxu 

and 

V(  V*u)  =  V2u  +  VxJ? 

Inserting  these  relations  into  Eq.  (2.54)  yields 

Po  ^  +  Vp'  =  (pB  +  |p)  V2u  -  |p0Vu2  -  p'  ~ 

— p0(f)Xu)  -f  ^Pb  +  ^  X  O  , 

where  u7  ~  u-u.  In  Appendix  B,  however,  it  is  shown  that  for  a  signal  prop¬ 
agating  far  from  any  boundaries  in  an  initially  irrotational  fluid,  the  vorticity 
is  a  higher-order  term.  The  vorticity  terms  in  the  above  relation  are  therefore 
neglected.9  The  remaining  terms  form  the  second-order  approximation  of  the 
momentum  equation, 

Po  +  Vp'  =  pQuV V2u  -  p'  ^  -  ipoVu2  ,  (2.80) 

9T!ic  interested  reader  is  referred  to  Appendix  B  where  vorticity  and  the  irrotational  flow 
assumption  are  discussed  at  length. 
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where  v  is  the  kinematic  viscosity  and  V  is  is  a  dimensionless  number  called  the 
viscosity  number, 


v 


(2.81) 


V  EE 


4  MB 
3  +  n 


(2.82) 


The  viscosity  number  indicates  the  relative  importance  of  the  bulk  viscosity  to  the 
shear  viscosity.  Estimates  of  the  viscosity  number  and  the  other  nondimensional 
numbers  that  appear  are  given  in  Appendix  B. 

A  quantity  that  aids  in  the  interpretation  of  local  nonlinear  effects  is  the 
Lagrangian  density  £, 


£  = 


(pT 

2poCq 


(2.83) 


The  two  terms  on  the  right-hand  side  are,  respectively,  the  kinetic  and  potential 
energy  densities.  It  turns  out  that  terms  that  contain  the  Lagrangian  density 
correspond  to  local  nonlinear  effects;  see,  for  example,  Aanonsen  et  al.  (1984). 
The  primary  advantage  of  the  term  £  is  that  it  makes  it  easy  to  see  what  terms  are 
to  be  dropped  when  the  wave  motion  is  progressive;  see  Hamilton  and  Blackstock 
(1988).  The  Lagrangian  density  and  local  effects  are  discussed  in  greater  detail 
in  later  chapters. 

We  now  introduce  the  Lagrangian  density  in  the  second-order  forms  of 
the  continuity  and  momentum  equations,  Eqs.  (2.53)  and  (2.80).  Recall  that  first- 
order  relations  may  be  used  to  modify  second-order  terms.  Thus  the  second-order 
terms  in  the  continuity  equation  may  be  rewritten  as 


p'V*  u  =  — 


1  w 

2p0Co  dt 


(2.84) 


and 


u-Vp' 


Po  du2 
2  dt 


(2.85) 


and  the  term  p'  (du/dt)  in  the  momentum  equation  may  be  rewritten  as 


du  _  V(p') 

dt  2p0Co 


Use  of  Eqs.  (2.84),  (2.85),  and  (2.86)  as  well  as  the  definition  of  the  Lagrangian 
density  in  the  continuity  and  momentum  equations  yields  the  desired  equivalent 
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second-order  forms  of  the  continuity  and  momentum  equations, 

dp'  „  1  dip')2  1  dC 

-nT  +  —  - 4 — XT - h  —  ~Z7 

at  p0c$  at  cq  at 


and 


O 

Po^  +  Vp'  =  p0t/VV2u  -  V£ 
at 


(2.87) 

(2. 88) 


Entropy,  state,  and  rate  equations 


The  rank-ordering  process  leads  to  the  linearized  form  of  the  expanded 
entropy  equation,  Eq.  (2.55).  Both  terms  in  Eq.  (2.55)  are  of  second  order.  To 
eliminate  T'  from  the  equation,  we  use  the  first-order  approximation  of  Eq.  (2.58), 


(2.89) 


Substituting  this  expression  into  the  entropy  equation  and  then  using  the  first- 
order  wave  equation  in  terms  of  the  density,  Eq.  (2.78),  leads  to 

d s'  _  K 

dt  p0To4 

This  relation  may  be  integrated  once  noting  that  the  integration  constant  must 
be  zero  to  satisfy  quiet  conditions.  Recalling  our  discussion  on  the  order  of  the 
difference  between  the  equilibrium  and  frozen  forms  of  a  thermodynamic  fluid 
property,  we  see  that  the  following  is  an  equally  valid  second-order  approximation 
of  the  entropy  in  a  relaxing  fluid: 

«  (dT\ 
poT0c%  V  dp 


s'  = 


\ 

Jo* 


(2.90) 


Note  that  the  effects  of  finite  amplitude,  viscosity,  and  relaxation  do  not,  within 
the  second-order  approximation,  contribute  to  the  entropy  production.  Thus,  the 
entropy  equation  need  only  be  considered  when  heat  conduction  is  of  interest. 

The  ranking  process  also  resulted  in  a  simplified  version  of  the  Taylor 
series  expansion  of  the  state  equation,  Eq.  (2.57).  Recalling  our  earlier  discus¬ 
sion  on  the  order  of  the  difference  between  the  equilibrium  and  frozen  forms  of 
thermodynamic  fluid  properties,  we  may  rewrite  Eq.  (2.57)  as  follows: 


p'  =  (ov +*'(?£ 


P,SJ 


+  fa  iL  (P'Y 

Po  2/1  KP  ! 


(2.91) 
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where  B/A  is  referred  to  as  the  parameter  of  nonlinearity, 


B 

A 


(2.92) 


The  defining  relation  for  the  frozen  sound  speed,  Eq.  (2.60),  was  used  in  obtaining 
Eq.  (2.91).  The  second-order  terms  in  Eq.  (2.91)  must  be  included  when  the 
effects  of  heat  conduction,  relaxation,  or  finite-amplitude,  respectively,  are  of 
interest. 

A  consistent  second-order  form  of  the  expanded  rate  equation,  Eq.  (2.56), 
is  now  obtained.  To  do  such,  we  require  some  information  about  the  order  of  q'. 
The  linearized  Taylor  series  expansion  of  q1  is 


q\p\p\s')  =  p‘ 


,  ( 

OP 


p-V  o 


-11 


P,*> 


,  US 


P,p) 


However,  use  of  Eq.  (2.66)  in  the  expansion  for  q'  leads  to 

dq 

dp 


q'{p\p',s')  =  [  7T 


Noting  Eq.  (2.72),  we  see  that  the  first  term  on  the  left-hand  side  of  the  above 
equation  is  of  second  order.  Since  s'  is  also  of  second  order,  the  quantity  q'  must 
be  of  second  order.  Thus  the  convective  term  in  the  rate  equation  must  be  of 
higher  order,  and  a  consistent  second-order  approximation  of  the  rate  equation 
is 


dq'  ^  q«  -  <?' 
dt  r0 


(2.93) 


The  rate  equation  must  be  included  when  relaxation  effects  are  of  interest. 


Combining  the  entropy,  state,  and  rate  equations 

The  second-order  approximations  of  the  entropy,  state,  and  rate  equa¬ 
tions  may  be  combined  into  a  single  equation  that  is  referred  to  as  the  modified 
state  equation.  Our  starting  point  is  the  time  derivative  of  the  rate  equation, 
Eq.  (2.93): 

r4Br)  +  l<’-’">  =  0  -  <2M> 

Our  procedure  is  as  follows:  The  entropy  equation  is  substituted  into  the  state 
equation  with  the  help  of  some  thermodynamic  identities  from  Appendix  A.  The 
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resulting  equation  is  then  rearranged  to  yield  an  expression  for  dq'  jdt.  This 
expression  and  an  expression  for  dq“'/dt  are  then  substituted  into  Eq.  (2.94). 
This  yields  the  modified  state  equation.  An  equivalent  form  of  the  modified  state 
equation  is  then  obtained  by  integration. 

The  entropy  and  state  equations,  Eqs.  (2.90)  and  (2.91),  respectively, 


are  now  combined  via  the  term  ^  ^  and  then  rearranged  to  yield  an 

expression  for  dq1  jdt.  Use  of  the  static  form  of  Eq.  ( A.  1 1 )  from  Appendix  A  and 
the  time  derivative  of  the  entropy  equation  yields 


l  dP 

\  ds1  k  p0  /  dT 

[  ds 

Jo  5F  =  M 

d2p' 

dt2 


However,  noting  the  definition  of  the  thermal  expansion  coefficient,  Eq.  (A. 22), 
and  using  the  static  form  of  a  thermodynamic  identity  from  Appendix  A, 
Eq.  (A. 26),  we  see  that 


( dP_  \  M 

l  *  Jo  at 


K 

Po 


d2p' 
dt 2  ’ 


(2.95) 


where  c„0  and  cPo  are  the  specific  heats  at  constant  volume  and  pressure,  respec¬ 
tively.  Equation  (2.95)  can  be  rearranged  and  placed  in  the  following  form: 


(dP  \  ds'  v{~i  —  1)  d2 p' 

Pr  ’ 

where  7  is  the  ratio  of  specific  heats  and  Pr  is  the  Prandtl  number, 


(2.96) 


7  =  ~  1  (2.97) 

C»>0 

Pr  =  ^  .  (2.98) 

K 

The  Prandtl  number  indicates  the  relative  importance  of  viscosity  to  heat  con¬ 
duction.  Estimates  of  the  Prandtl  number  are  given  in  Appendix  B.  Substituting 
Eq.  (2.96)  into  the  time  derivative  of  the  state  equation,  Eq.  (2.91),  and  rear¬ 
ranging  leads  to  an  expression  for 


(OP  \  dq'  _  dp'  dp'  1/(7-!)  d2p'  Cp  d  ,dp' 

\  dq  dt  dt  ’  dt  Pr  dt2  Po  A  9  dt 


(2.99) 
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An  expression  for  is  obtained  from  the  time  derivative  of  the  lin¬ 
earized  Taylor  series  expansion  of  qm', 


dq *'  /  dq' 


dt 


dp 


\dj_  (d£_ 

Jo  dt  +  \ds 


P/  0 


d£_ 

dt 


The  second-order  term  on  the  right-hand  side  drops  out  because  the  change  in 
the  entropy  at  equilibrium  is  zero.  Thus,  the  above  relation  simplifies  to 


dq*  _  /  dq * 
[dp 


dt 


dp1 


3 /  0 


dt 


(2.100) 


An  expression  for  §-t{q'  —  q*')  may  be  obtained  by  combining  Eqs.  (2.100) 
and  (2.99), 


dP 

dq 


P'aJ  o 


d ,  /  dp’ 

Wt(q  ~q  )  =  Si 


dP 


(C)2+  ITT 


dq 


p,s/  0 


'dr 

.  dp 


3 /  0 


dj_ 

dt 


^(7-1)  d2p'  _  c|  5  ,d£_ 
Pr  dt 2  p0  AP  dt 


However,  using  Eqs.  (2.66),  (2.67),  and  (2.68),  we  see  that  the  coefficient  of  the 
second  term  on  the  right-hand  side  is  simply  the  equilibrium  sound  speed.  Thus, 
we  obtain 


dP_ 

dq 


'p.v  o 


a.,  .,  dP'  0JaP'  Ki-ijav  4B  m 

-gid-i  )  =  W-M  gr — —W~*a'W  '(2'101) 


The  modified  state  equation  is  now  formed.  Substitution  of  Eqs.  (2.99) 
and  (2.101)  into  Eq.  (2.94)  yields  an  expression  that  may  be  integrated  once  with 
respect  to  time.  The  integration  constant  must  be  zero  in  order  to  satisfy  quiet 
conditions.  The  result  is 


d  i 

T°  Ft1 


(#)V 


-  1)  ^  _  £o  B_,  J\2 
Pr  dt  p0  2/i  , 


+  [p  “  (<£)  V  - 


Hit  ~\)dp'  ,,2 

Pr  dt  po  2 A  [P  ’ 


0  .  (2.102) 


Equation  (2.102)  is  the  modified  state  equation,  a  combination  of  the  entropy, 
state,  and  rate  equations.  The  basis  for  our  earlier  comments  about  the  rela¬ 
tionship  between  the  equilibrium  and  frozen  sound  speeds  and  the  relaxation 
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time/signal  frequency  is  now  apparent.  If  the  relaxation  time  is  very  large, 
Eq.  (2.102)  simplifies  into  an  expression  that  depends  only  on  the  frozen  sound 
speed.  Conversely,  for  a  small  relaxation  time,  only  the  equilibrium  sound  speed 
appears. 

Equation  (2.102)  may  be  integrated  by  following  the  procedure  of 
Rudenko  and  Soluyan  (1977,  p.  83).  Use  of  Eq.  (2.69)  leads  to  the  following: 


To 


d  ( „/  /  _o \ 2  j  "(7-1)  dp'  eg  B  \ 

) 

,  (_>  /  ji\2  "(7-1)  dp’  eg  B 

+  (p  -  <s>)  p  -  -pT"  -  -  2a(p  > 


idp' 


-Tn0T0c0~  -  0 


The  above  equation  is  in  the  form  of  a  linear  first-order  differential  equation  that 
can  be  solved  by  means  of  an  integrating  factor  [see,  for  example,  Kreider  et  al. 
(1966,  p.  97)].  The  result  is 


,  =  (c.)V  +  ^li)f  +fJA(pr  +  ma4J^e-^dy  .  (2.103, 


The  integral  in  Eq.  (2.103)  may  be  integrated  by  parts,  and  the  resulting  expres¬ 
sion  is 


+  '  (2I04) 


Equation  (2.104)  is  referred  to  as  the  integral  of  the  modified  state  equation. 
Equation  (2.104)  corresponds  to  Eq.  (IV-1.20)  in  Rudenko  and  Soluyan’s  book 
(1977)  except  that  Eq.  (2.104)  also  accounts  for  heat  conduction.  The  second, 
third,  and  fourth  terms  on  the  right-hand  side  of  the  Eq.  (2.104)  are  associated 
with  heat  conduction,  finite-amplitude  effects,  and  relaxation,  respectively. 

For  the  special  case  of  u>t0  <C  1  it  is  difficult  to  distinguish  between  other 
loss  terms  and  the  leading-order  effects  of  relaxation.  Integrating  Eq.  (2.103)  by 
parts  and  neglecting  terms  of  order  O[e(u;ro)2]  and  higher  yields 


p'  =  (c2)V+(drzil  +  ^,0)|:  +  iA(/)!  .  (2.  ,05) 


In  particular,  the  leading  order  effects  of  relaxation  are  difficult  to  distinguish 
from  bulk  viscosity  because  it  is,  apparently,  difficult  to  devise  a  measurement 
technique  which  isolates  the  bulk  viscosity  from  the  relaxation. 


i 
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Summary  of  second-order  results 

In  this  section,  consistent  second-order  approximations  of  the  continuity, 
momentum,  entropy,  state,  and  rate  equations  for  a  homogeneous,  thermoviscous 
fluid  with  a  single  relaxation  mechanism  were  obtained.  The  equations  are,  re¬ 
spectively,  Eqs.  (2.53),  (2.80),  (2.90),  (2.91),  and  (2.93).  The  Lagrangian  den¬ 
sity  was  introduced  into  the  continuity  and  momentum  equations,  and  equivalent 
forms  of  the  equations,  Eqs.  (2.87)  and  (2.88),  were  obtained.  The  entropy,  state, 
and  rate  equations  were  combined  to  form  a  new  equation,  which  is  called  the 
modified  state  equation.  Two  equivalent  forms  of  the  modified  state  equation, 
Eqs.  (2.102)  and  (2.104),  were  obtained. 


2-6  Small-Signal  Wave  Motion  in  a  Homogeneous,  Thermoviscous, 
Relaxing  Fluid 


In  this  section  the  basic  equations  are  simplified  for  the  case  of  small- 
signal  propagation  in  a  homogeneous,  thermoviscous,  relaxing  fluid.  Rather  than 
re-deriving  the  relations  from  the  basic  equations,  we  merely  drop  the  quadratic 
nonlinearity  terms  from  the  second-order  forms  of  the  basic  equations.10  Since 
no  quadratic  nonlinearity  terms  exist  in  either  the  entropy  or  rate  equations, 
Eqs.  (2.90)  and  (2.93),  respectively,  they  remain  unchanged.  Neglecting  quadratic 
nonlinearity  terms  in  the  continuity  yields  the  same  equation  that  was  obtained 
for  small-signal  propagation  in  a  lossless  fluid,  Eq.  (2.73).  Neglecting  quadratic 
nonlinearity  terms  in  the  momentum  and  state  equations  leads  to  the  following: 


du  j  , 

Po  +  Vp  =  p0t/VV2u 
ot 


(2. 10C) 


and 


p'  =  (OV  +  s' 


dP_ 

ds 


+  ?' 


(2.107) 


’"If  one  wanted  to  derive  approximate  forms  of  the  basic  equations  valid  for  small-signal 
propagation  in  a  homogeneous,  thermoviscous  fluid  with  a  single  relaxation  mechanism  starting 
from  the  basic  equations,  the  procedure  is  identical  to  that  used  for  finite-amplitude  signals 
except  that  nonlinear  terms  may  be  neglected  from  the  outset.  It  turns  out,  however,  the  effort 
saved  in  neglecting  nonlinearity  from  the  outset  is  small. 
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Dropping  the  quadratic  nonlinearity  term  from  the  modified  state  equations, 
Eqs.  (2.102)  and  (2.104),  yields 


_  9  (  ,  /oox2  /  "(7-1)  dp 

T°  ax  p  -  (c° ]  p  H 


,  (  /  /  Os2  /  "(7  -  1  )dp'\ 

-)  +  {”  -{c°)p  —prsi) 


_/  t  _oo\2 _/  ,  "(7-1)  dp' 


dy  . 


=  0 
(2.108) 

(2.109) 


The  forms  of  the  continuity,  momentum,  entropy,  state,  and  rate  equa¬ 
tions  that  are  appropriate  for  small-signal  propagation  in  a  homogeneous,  ther- 
moviscous  fluid  with  a  single  relaxation  mechanism  are,  respectively,  Eqs.  (2.73), 
(2.106),  (2.90),  (2.107),  and  (2.93).  The  appropriate  forms  of  the  modified  state 
equation  are  Eqs.  (2.108)  and  (2.109). 


2-7  Epilogue 

In  this  chapter,  the  basic  equations  for  a  homogeneous  thermoviscous 
fluid  with  a  single  relaxation  mechanism  were  presented.  The  deviation  from 
equilibrium  was  assumed  to  be  small,  and,  accordingly,  linear  relations  between 
the  thermodynamic  fluxes  and  forces  were  used.  The  possibility  of  a  cross-effect 
between  the  bulk  viscosity  and  relaxation  was  neglected.  Assumptions  were  then 
made  about  the  magnitude  of  the  acoustic  signal,  the  magnitude  of  the  transport 
coefficients,  and  the  amount  of  dispersion.  A  ranking  system  was  introduced 
and  used  to  simplify  the  basic  equations  into  specialized  forms  for  (1)  small 
signals  in  lossless  fluids,  (2)  small  signals  in  homogeneous,  thermoviscous,  relaxing 
fluids,  and  (3)  finite-amplitude  signals  in  homogeneous,  thermoviscous,  relaxing 
fluids.  Moreover,  the  Lagrangian  density  was  introduced  into  the  continuity  and 
momentum  equations,  and  the  entropy,  state,  and  rate  equations  were  combined 
to  form  a  new  equation,  the  modified  state  equation. 


CHAPTER  3 


WAVE  EQUATIONS  FOR  FINITE-AMPLITUDE  SIGNALS  IN  A 
HOMOGENEOUS,  THERMOVISCOUS  FLUID  WITH  A  SINGLE 
RELAXATION  MECHANISM 


3-1  Introduction 

In  this  chapter  second-order  forms  of  the  wave  equation  for  finite- ampli¬ 
tude  signals  in  a  thermoviscous,  relaxing  fluid  are  developed.  The  wave  equations 
are  expressed  in  terms  of  the  acoustic  pressure  p\  the  velocity  potential  <j>,  and 
two  new  variables — the  modified  acoustic  pressure  P'  and  the  modified  velocity 
potential  <I>.  The  modified  variables  are  defined  in  such  a  way  that  the  Lagrangian 
density  C  does  not  appear  in  a  wave  equation  that  is  expressed  in  terms  of  P'  or 
Use  of  the  modified  pressure  and  modified  velocity  potential  was  developed 
by  Naze  Tjptta  and  Tjptta;  see,  for  example,  Aanonsen  et  al.  (1984). 

3-2  Wave  Equation  in  Terms  of  the  Acoustic  Pressure  p' 

The  wave  equation  in  terms  of  the  acoustic  pressure  p'  is  derived  in  this 
section.  The  procedure  is  as  follows:  The  continuity  and  momentum  equations 
containing  the  Lagrangian  are  combined  to  form  one  equation.  This  equation 
and  the  linear  state  equation  are  then  used  to  eliminate  the  density  fluctuation 
in  the  modified  state  equation.  The  result  is  the  wave  equation  in  terms  of  the 
acoustic  pressure  p' .  After  a  modicum  of  rearrangement,  the  wave  equation  may 
be  integrated  and  placed  in  an  equally  valid  form  that  is  referred  to  as  the  integral 
of  the  wave  equation  in  terms  of  the  acoustic  pressure  p' . 

Combining  the  continuity  and  momentum  equations,  Eqs.  (2.87)  and 
(2.88),  is  a  straightforward  procedure:  The  time  derivative  of  the  continuity  equa¬ 
tion  is  subtracted  from  the  divergence  of  the  momentum  equation.  The  viscosity 
term  may,  however,  be  rearranged  so  that  it  may  later  be  readily  combined  with 
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the  heat  conduction  term.  Use  of  first-order  forms  of  the  continuity  equation, 
the  state  equation,  and  the  wave  equation  in  p',  Eqs.  (2.73),  (2.76),  and  (2.77), 
respectively,  leads  to  the  following  first-order  relation: 


V2(V-u) 


i  ay 

po4  dt 3 


Use  of  the  above  yields  the  desired  form  of  the  combined  continuity  and  momen¬ 
tum  equation, 


ay 
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(3.1) 


The  wave  equation  in  terms  of  the  acoustic  pressure  is  obtained  by  elim¬ 
inating  the  density  fluctuation  p’  in  first  and  second-order  terms  (using,  respec¬ 
tively,  Eqs.  (3.1)  and  the  linear  state  equation,  Eq.  (2.76))  in  the  modified  state 
equation,  Eq.  (2.102).  The  result  is 
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where  b/ p0  is  the  diffusivity  of  sound  (Lighthill  1956,  Eq.  (19)]  and  0  is  the 
coefficient  of  nonlinearity, 


6  =  ^(V  +  I_i) 


4 

=  PB  +  gP  +  K 


(3.3) 


0  = 


(3.4) 


Note  that  if  the  relaxation  time  r0  is  very  large,  Eq.  (3.2)  simplifies  to  a  wave 
equation  that  involves  only  the  frozen  sound  speed.  On  the  other  hand,  if  the 
relaxation  time  is  very  short,  the  converse  occurs;  Eq.  (3.2)  simplifies  to  a  wave 
equation  that  involves  only  the  equilibrium  sound  speed.  The  third,  fourth,  and 
fifth  terms  in  both  major  components  of  Eq.  (3.2)  represent,  respectively,  the 
effects  of  viscosity  and  heat  conduction,  local  nonlinearity,  and  growing  non¬ 
linearity.  Considered  together,  the  first  and  second  terms  in  both  the  major 
components  represent  linear,  lossless  wave  motion  at  the  frozen  and  equilibrium 
sound  speeds,  respectively. 
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Equation  (3.2)  may  be  integrated  in  the  same  way  that  the  modified  state 
equation,  Eq.  (2.102),  was  integrated.  First,  however,  it  must  be  rearranged  using 
the  following: 
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Substitution  of  Eq.  (3.5)  into  Eq.  (3.2)  leads  to 
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where  O2  is  called  the  d’Alembertian  operator  and  the  subscript  0  indicates  that 
this  d’Alembertian  operator  uses  the  equilibrium  sound  speed, 

2  _2  Id2 
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°  ( c°0)2dt 2  ‘ 

As  was  the  case  with  the  modified  state  equation,  Eq.  (3.6)  is  in  the  form  of  a 
linear  first-order  differential  equation  that  may  be  solved  by  way  of  an  integrating 
factor  (see,  for  example,  Kreider  et  al.  (1966,  p.  97)).  The  result  is 
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which  is  here  called  the  integral  of  the  wave  equation  in  terms  of  the  acoustic  pres¬ 
sure.  The  equation  is  a  second-order  approximation  that  consistently  accounts 
for  the  effects  of  viscosity,  heat  conduction,  relaxation,  and  nonlinearity.  Each  of 
the  five  terms  in  Eq.  (3.7)  represents  a  specific  acoustic  phenomenon.  The  three 
terms  on  the  left-hand  side  represent,  respectively,  (1)  small-signal  wave  motion 
in  a  lossless  fluid,  (2)  the  dissipation  of  the  signal  that  is  caused  by  viscosity 
and  heat  conduction,  and  (3)  the  effects  of  relaxation.  The  two  terms  on  the 
right-hand  side  of  the  equation  represent,  respectively,  (1)  the  local  nonlinear 
effects  and  (2)  the  cumulative  nonlinear  effects.  Clearly,  if  either  thermoviscous 
or  relaxation  effects  are  negligible,  the  appropriate  form  of  the  wave  equation 
is  obtained  by  setting  b  or  m,  respectively,  to  zero.  Similarly,  if  nonlinearity  is 
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not  of  interest,  then  the  appropriate  form  of  the  wave  equation  is  obtained  by 
neglecting  the  two  terms  on  the  right-hand  side. 

An  equation  that  is  equivalent  to  Fq  (3.7)  may  be  obtained  by  combining 
Eq.  (3.1)  with  the  integral  of  the  modified  state  equation,  Eq.  (2.104).  This 
eliminates  the  need  for  the  foregoing  integration,  but  Eq.  (3.2)  would  not  be 
obtained. 

As  we  noted  in  the  previous  chapter  when  examining  the  modified  state 
equation,  it  is  difficult  to  distinguish  between  other  loss  terms  and  the  leading- 
order  effects  of  relaxation  for  the  special  case  of  ojt0  1.  Integrating  the  re¬ 
laxation  term  in  Eq.  (3.7)  by  parts  and  neglecting  terms  of  order  O[e(u>T0)2}  and 
higher  yields 
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Note  that  the  definition  of  the  diffusivity  b  may  be  generalized  to  include  the 
effects  of  relaxation  in  this  special  case  of  u>t0  1. 


3-3  The  Wave  Equation  in  Terms  of  the  Velocity  Potential  </> 

The  particle  velocity  may  be  expressed  as  the  gradient  of  a  scalar  po¬ 
tential  because  the  fluid  is,  within  the  second-order  approximation,  irrotational 
(see  Appendix  B), 

u  =  .  (3.9) 

The  gradient  of  the  velocity  potential  may  therefore  be  used  to  replace  the  particle 
velocity  in  all  equations.  Moreover,  first  and  second-order  expressions  for  the 
pressure  and  the  wave  equation  in  terms  of  <f>  may  be  obtained.  An  expression 
for  the  Lagrangian  £  in  terms  of  p  may  also  be  obtained. 

Before  obtaining  the  second-order  form  of  the  wave  equation  in  <f>,  we 
develop  the  other  aforementioned  relations.  Integration  of  the  first-order  wave 
equation  in  terms  of  the  particle  velocity  and  setting  the  arbitrary  static  value  of 
the  potential  <f>o  to  zero  leads  to 

=  °  •  (3.10) 

The  second-order  form  of  the  momentum  equation,  Eq.  (2.88),  becomes  (after 
integration) 

p' =  ~Po^  +  Port?2*  ~  C  , 


(3.11) 
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the  first-order  approximation  of  which  is 

,  d<j> 
p  = 


(3.12) 


The  Lagrangian  £,  which  is  defined  in  Eq.  (2.83),  may  also  be  expressed  in  terms 
of  <f>, 
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(3.13) 


Note,  however,  that  the  term  d2<fr  expands  identically  as 
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where  D2  is  the  d’Alembertian  operator  using  the  small-signal  sound  speed, 
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Noting  Eq.  (3.10),  we  see  that,  to  second  order,  the  Lagrangian  may  be  expressed 
in  terms  of  the  velocity  potential  as  follows: 


•  (3.15) 

With  the  above  relations,  we  may  now  readily  express  the  second-order 
wave  equation,  Eq.  (3.7),  in  terms  of  <j>.  Substitution  into  first  and  second-order 
terms  involving  the  pressure  using  Eqs.  (3.11)  and  (3.12),  respectively,  yields 


d  n2l  ,  b  d4^ 
dt Uo<t>  +  p04  dt 4 


+ 


12°  /  e-(t-y)/ro 

c2  J  dy4 

—  OO 


dy 


2-<Lr  +  fUL(?£\ 

Pocldt 2  c^dt2  j 


(3.16) 


A  higher-order  term,  (i/V/cjOn2^,  was  neglected  in  obtaining  the  above  relation. 
In  order  to  integrate  the  above  equation,  we  rearrange  the  relaxation  term  using 
the  following: 
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Equation  (3.17)  can  be  verified  by  integrating  the  left-hand  side  by  parts  and  by 
applying  Leibnitz’s  theorem  to  the  right-hand  side.  Substituting  Eq.  (3.17)  into 
Eq.  (3.16)  and  then  integrating  with  respect  to  time  leads  to 
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By  expanding  the  nonlinearity  terms  using  the  definition  of  f3  and  £,  Eqs.  (3.4) 
and  (3.13),  we  see  that  Eq.  (3.18)  may  be  written  as 
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Equations  (3.18)  and  (3.19)  are  equivalent  forms  of  the  wave  equation  in  terms 
of  the  velocity  potential  </>.  The  terms  in  Eq.  (3.18)  have  the  same  meaning  as 
the  corresponding  terms  in  Eq.  (3.7). 


3-4  Wave  Equations  in  Terms  of  the  Modified  Pressure  P'  and  the 
Modified  Velocity  Potential  $ 


In  this  section,  the  wave  equation  is  expressed  in  terms  of  two  new  vari¬ 
ables,  the  modified  acoustic  pressure  P'  and  the  modified  velocity  potential  4>, 
which  are  defined  below. 

To  motivate  the  definition  of  P'  and  to  find  the  wave  equation  in  terms 
of  it,  we  start  with  the  wave  equation  in  terms  of  the  acoustic  pressure,  Eq.  (3.7). 
Equation  (3.7)  may  be  rearranged  with  the  aid  of  Eq.  (3.15)  to  obtain 
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The  modified  acoustic  pressure  P'  is  now  defined  as 

■  (3-2,) 

To  first  order,  the  modified  and  the  ordinary  acoustic  pressures  are  equal.  The 
difference,  a  second-order  term,  accounts  for  local  nonlinear  effects.  The  wave 
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equation  in  terms  of  modified  acoustic  pressure  is,  accordingly, 
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In  Eq.  (3.22),  no  single  term  explicitly  accounts  for  the  local  effects.  The  effects 
are  instead  accounted  for  in  the  definition  of  the  modified  acoustic  pressure. 

An  expression  for  the  modified  pressure  P'  in  terms  of  the  velocity  poten¬ 
tial  <{>  motivates  the  definition  of  the  modified  velocity  potential  4>.  If  Eq.  (3.11) 
is  substituted  into  the  definition  of  P1,  Eq.  (3.21),  the  following  expression  is 
obtained: 

The  modified  velocity  potential  is  now  defined  as 

1  d 


$  = 

2  eg  dt 


(3.23) 


Thus  the  modified  acoustic  pressure  is  related  to  the  modified  velocity  potential 

by 
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P1  =  -p0-~  +  pQi/W2$ 
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(3.24) 


Inserting  Eq.  (3.24)  into  Eq.  (3.22)  and  integrating  noting  Eq.  (3.17) 
leads  to  the  wave  equation  in  terms  of  the  modified  velocity  potential, 
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As  with  the  wave  equation  in  terms  of  the  modified  pressure,  Eq.  (3.22),  no  single 
term  in  Eq.  (3.25)  accounts  for  local  nonlinear  effects.  These  effects  are  instead 
accounted  for  in  the  modified  velocity  potential  $. 

Equation  (3.25)  may  be  obtained  directly  from  the  wave  equation  in 
terms  of  the  velocity  potential,  Eq.  (3.19).  Such  an  approach  yields  the  same 
motivation  for  the  definition  of  the  modified  velocity  potential. 


3-5  Summary 

In  this  chapter,  wave  equations  in  terms  of  the  acoustic  pressure  //, 
the  velocity  potential  <j>,  and  two  new  variables — the  modified  acoustic  pressure 
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P'  and  the  modified  velocity  potential  $ — were  derived.  They  are,  respectively, 
Eqs.  (3.7),  (3.19),  (3.22),  and  (3.25).  The  wave  equations  in  terms  of  the  modified 
variables  contain  one  less  term  than  their  ordinary  counterparts  because  the 
modified  variables  themselves  account  for  the  local  nonlinear  effects. 


CHAPTER  4 


ON  THE  BOUNDARY  CONDITIONS  AT  THE  INTERFACE 
BETWEEN  TWO  LOSSLESS,  IMMISCIBLE  FLUIDS 


4-1  Introduction 

In  this  chapter,  the  boundary  conditions  at  the  interface  between  two 
inviscid,  immiscible  fluids  are  examined.  The  fluids  are  assumed  to  be  initially 
quiet,  and  the  interface  is  assumed  to  be  initially  planar  and  coincident  with 
the  2  =  0  plane.  Moreover,  we  neglect  the  effects  of  surface  tension  and  body 
forces  at  the  interface.  Expressions  that  are  correct  to  second  order  are  obtained 
for  each  of  the  two  boundary  conditions — the  kinematic  condition  and  the  force 
balance  condition  (Newton’s  second  law).  Although  our  analysis  in  later  chapters 
is  restricted  to  two  dimensions  (plane  waves  obliquely  incident  on  an  initially 
plane  interface),  the  approximate  forms  of  the  boundary  conditions  developed  in 
this  chapter  are  for  three  dimensions. 

The  superscripts  I  and  II  are  used  to  indicate  the  different  fluids  through¬ 
out  the  remainder  of  this  dissertation. 

4-2  Kinematic  Condition 

The  author  can  think  of  no  better  way  to  express  the  kinematic  boundary 
condition  than  to  quote  from  Prandtl  and  Tietjens  (1934), 

The  kinematic  boundary  conditions  at  the  surface  of  contact  between 
a  liquid  and  a  solid  body,  and  also  between  two  immiscible  fluids 
(water  and  oil,  water  and  air,  etc.),  must  clearly  be  such  that  neither 
vacuum  nor  interpenetration  can  occur.  The  necessary  consequence 
of  this  is  that  the  normal  components  of  the  velocities  of  the  two 
media  are  equal  on  each  side  of  the  surface  of  contact,  ...  . 
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Obtaining  the  normal  component  of  the  particle  velocities  at  the  interface  requires 
knowledge  of  the  normal  to  the  interface.  Consider  Figure  4.1,  which  depicts  the 
interface  between  fluids  I  and  II  at  time  t,  where  t  >  0.  Because  the  interface 
initially  lies  in  the  z  —  0  plane,  the  equation  for  the  interface  may  be  written  as 


z-f(x,y,t)  =  0  ,  (4.1) 

or,  alternatively,  as 

F{x,y,z,t)  =  0  .  (4.2) 

A  vector  normal  to  the  interface  N(x,y,  z,t)  is  [see,  for  example,  Thomas  and 


Figure  4.1  Interface  between  fluids  I  and  II  at  time  t  >  0 
Finney  (1979)] 

N{x,y,z,t)  =  -j^  +  k  ’  (4-3) 

or,  in  terms  of  the  function  F, 

N  =  VF  .  (4.4) 


Note  that  for  later  convenience  N  is  used  rather  than  the  unit  normal  V.F/jV.F|. 
If  the  interface  has  a  velocity  of  Uf,  then  the  normal  velocity  of  the  interface  and 
the  normal  velocity  of  a  particle  at  the  interface  on,  say,  the  fluid  I  side  must 
match, 

Up*  N  =  u1*  N  ,  (4.5) 

where  u1  is  the  velocity  of  fluid  I.  The  situation  is  the  same  on  the  fluid  II  side, 

up-  N  =  u11-  N  ,  (4.6) 

where  u11  is  the  velocity  of  fluid  II.  The  kinematic  boundary  condition  between 
the  two  inviscid,  immiscible  fluids  may  thus  be  written  as  follows:1 

u1- N  =  un- N  on  F(x,y,  z,t)  =  0  .  (4.7) 

For  future  reference,  we  note  that  u1  and  u11  may  be  written  in  terms  of  their  x, 
y ,  and  2  components  as  follows: 

u1  =  iu1  +  jv1,  +ku>1  ,  (4.8) 

u11  =  itt11  +  jt>n  +  ku;11  .  (4.9) 

It  turns  out  that  it  is  important  to  know  what  happens  to  a  particle 
that  is  initially  on  the  interface.  Physically,  we  may  reason  that  if  the  normal 
velocity  of  a  fluid  particle  at  the  interface  is  equal  to  the  normal  velocity  of  the 
interface  itself,  then  a  particle  that  is  initially  on  the  interface  must  remain  on 
the  interface.  This  may  be  shown  mathematically  by  following  the  developments 
of  Myers  (1980)  and  Lamb  (1932).  The  velocity  of  a  point  on  the  interface, 
r  =  ii-fj?/-f-k2,  is  given  by 

dF  dr  _  _ 

~m+Tt'VF  =  0  on  F(x^y^z^)- 0  ■ 

Noting  that  dr/dt  =  Up  and  substituting  using  Eq.  (4.5)  yields 
dF  DF 

—  +  u'-VF  =  —  =  0  on  F(x,y,z,t)  =  0  . 

Following  the  proof  given  by  Lord  Kelvin,  Lamb  (1932,  p.  7)  shows  that,  for  an 
inviscid  fluid,  an  interface  whose  equation  satisfies  DF/ Dt  =  0  always  consists  of 


'Note  that  if  the  fluid  were  viscous,  a  second  kinematic  condition  would  emerge.  This 
kinematic  condition  would  govern  the  particle  velocity  tangential  to  the  interface. 
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the  same  fluid  particles,  that  is,  that  a  particle  that  is  initially  on  the  interface 
remains  on  the  interface. 

Note  that  the  kinematic  boundary  condition  as  expressed  in  Eq.  (4.7)  is 
to  be  applied  on  the  interface,  the  location  of  which  is  time  varying.  Consequently, 
the  normal  to  the  interface  is  also  time  varying.  However,  because  the  fluid 
is  assumed  initially  quiet,  any  motion  of  the  interface  must  be  caused  by  the 
impinging  acoustic  signal.  It  would  appear,  therefore,  to  be  possible  to  describe 
the  motion  of  the  normal  to  the  interface  in  terms  of  the  impinging  acoustic 
signal.  Since  the  kinematic  boundary  condition  already  depends  directly  on  the 
particle  velocity,  an  expression  for  the  normal  in  terms  of  the  particle  velocity  is 
sought. 

It  turns  out,  however,  to  be  convenient  to  find  first  an  expression  for 
the  normal  in  terms  of  the  particle  displacement.  To  that  end,  the  expression 
for  the  normal  to  interface  given  in  Eq.  (4.3)  is  used,  and  expressions  for  df /dx 
and  df  /dy  are  obtained  in  terms  of  particle  displacement.  The  displacement  of 
a  fluid  particle  that  is  at  location  x,y,z  at  time  t  is  denoted  £, 


=  iA(x,y,z,t)  +  jB(x,y,z,t)  +  kC{x,y,z,t)  ,  (4.10) 

where  A,  B ,  and  C  represent,  respectively,  the  components  of  displacement  from 
the  initial  (rest)  position  of  the  particle.  The  particle  displacement  and  the 
particle  velocity  are  related  in  the  following  manner: 

u  =  §3¥  +  u-v<  '  <4'u> 

Since  a  particle  that  is  initially  on  the  interface  stays  on  the  interface,  the  dis¬ 
placement  of  the  interface  given  by  Eq.  (4.1)  and  the  z  component  of  Eq.  (4.10) 
for  a  particle  on  the  interface  must  be  equal,  that  is, 


z  =  /(x,  y,  t)  =  C  {x,y,z,t)\zsf(x  v  t) 


(4.12) 


The  following  relationships  between  the  derivatives  of  /  and  C  may  be  readily 
derived: 


dl 

dx 


9C 

dx 


i  _  bd  ’ 
1  dz 


(4.13) 


i-f? 


(4.14) 
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Inserting  Eqs.  (4.13)  and  (4.14)  into  Eq.  (4.3)  and  rearranging  yields  the  desired 
form  of  the  normal  in  terms  of  the  particle  displacement, 


.dC  .dC  ,  (,  dC 

-JaF  +  k(1-a7 


(4.15) 


*=/(*.»  .0 


Since  no  approximations  have  been  made,  Eq.  (4.15)  is  an  exact  expression  for 
the  normal  to  the  moving  interface. 

Next,  an  expression  for  the  normal  that  is  correct  to  first  order  (linear 
in  the  acoustic  variables)  is  obtained.  Performing  a  binomial  expansion  of  the 
denominator  of  Eq.  (4.15)  and  retaining  only  linear  terms  leads  to 


.dC  .dC  , 

N  -  -i  -z - J  -z — h  k 

dx  °y  lz=/(r,y,0 

Noting  that  the  first-order  approximation  of  Eq.  (4.11)  is 

dt 


u  = 


dt 


(4-16) 


(4.17) 


and  that  the  fluid  is  assumed  to  be  initially  at  rest,  we  see  that  the  first-order 
approximation  of  C(x,y,  z,t)  is 


C(x,y,z,t)  =  J  w(x,y,z,t) 


dt 


*=/(*.».«) 


(4.18) 


Use  of  Eq.  (4.18)  leads  to  the  first-order  approximation  of  the  normal  expressed 
in  terms  of  the  particle  velocity, 


lJ-y 


dt  - f  k 


z=/(x,y,t) 


(4.19) 


An  expression  of  the  kinematic  boundary  condition  applied  on  the  mov¬ 
ing  interface  that  is  correct  to  second  order  is  now  obtained.  Substitution  of 
Eq.  (4.19)  into  the  kinematic  boundary  condition,  Eq.  (4.7),  yields  the  following 
expression,  which  is  correct  to  second  order: 


,  [dw1  ,  ,  fdw1  ,  j 

— ul  /  — —  dt  —  vl  — —  dt  +  tv1  = 
J  Ox  J  ay 


'dw1 

By 

'  chc11 


_  U”  / dt  ~  J~1T~  dt  +  14,11  on  z  =  ■  (4-20) 
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Since  the  first-order  approximation  of  Eq.  (4.20)  is 

wl  =  w11  on  z  —  f(x,y,t )  , 


we  see  that  Eq.  (4.20)  may  be  written 

(u11  -u1)  J  ^—dt  +  (v11  -v1)  J  ^—dt  +  wl  =  iu11  on  z  =  f(x,y,t)  .(4.21) 


Equation  (4.21)  is  the  kinematic  boundary  condition  applied  on  the  moving  in¬ 
terface  correct  to  second  order. 

An  expression,  which  is  correct  to  second  order,  for  the  kinematic  bound¬ 
ary  condition  to  be  applied  at  z  =  0  rather  than  on  the  moving  interface  2  = 
f(x,  y,  t )  is  now  obtained.  This  step  is  required  because  the  wave  equation  is 
most  conveniently  solved  subject  to  boundary  conditions  on  boundaries  that  do 
not  move.  Application  of  a  Taylor  series  expansion  to  Eq.  (4.21)  and  retaining 
terms  up  to  second  order  yields 


dw11 

dy 


dt 


+  wl  +  f(x,y,t)  —  -  wn  +  f(x,y,t) -Q—  on  z  ~  0  . 

Substitution  of  the  leading-order  term  in  the  Taylor  series  expansion  of  Eq.  (4.12) 
into  Eq.  (4.18)  leads  to 


/(x,y,t)  =  J w(x,  y,  z,t)  < 
Combining  the  two  foregoing  relations  yields 

(u11  -  ul)  J  dt  +  {v11  ~  vl)  J  dt  +  ( 


z=0 


(4.22) 


dw1  dwu ' 


w11  dt 


dz  dz 
=  u,11  —  ui1  on  z  =  0  . 


(4.23) 


Equation  (4.23)  is  the  kinematic  boundary  condition  that  applies  at  z  =  0  correct 
to  second  order.  Dropping  the  second-order  terms  yields 

w1  =  w11  on  2  =  0  ,  (4-24) 


which  is  the  well  known  small-signal  form  of  the  kinematic  boundary  condition: 
the  normal  particle  velocities  balance  at  the  interface.  Note  that  in  the  small- 
signal  approximation  the  motion  of  the  interface  is  assumed  negligible. 
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4-3  Force  Balance  Condition 

The  second  boundary  condition — the  force  balance  condition — is  ob¬ 
tained  by  applying  Newton’s  second  law  to  the  interface.  Consider  Fig.  4.2,  which 
shows  a  material  region  (abbreviated  MR)  that  spans  the  interface.  Newton’s  sec- 


Figure  4.2  Material  region  that  spans  the  interface 

ond  law  may  be  expressed  in  a  general  form  as  follows  [see,  for  example,  Panton 
(1984,  p.  107)]: 

—  /  pudV  =  f  pF  dV  +  [  RdS  on  z  =  f(x,y,t)  ,  (4.25) 

dt  J MR  J MR  J MR 

where  F  is  the  body  force  per  unit  mass  vector,  and  R  is  the  surface  force  per 
unit  area  vector.  Since  body  forces  (gravity,  for  example)  are  neglected,  the 
first  term  on  the  right-hand  side  drops  out.  Moreover,  the  surface  force  vector 
reduces  to  the  pressure  applied  normal  to  the  material  region  because  viscosity 
and  surface  tension  are  assumed  negligible.  If  the  material  region  is  collapsed 
about  the  massless  interface,  the  left-hand  side  of  the  Eq.  (4.25)  goes  to  zero, 
which  leaves 


pl  =  pn  on  z  =  f(x,  y,  t) 


(4.26) 
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Equation  (4.26)  is  the  exact  boundary  condition  that  is  to  be  applied  on  the 
moving  interface. 

An  expression  for  the  pressure  balance  boundary  condition  on  z  —  0  that 
is  correct  to  second  order  may  be  obtained  by  applying  a  Taylor  series  expansion 
to  Eq.  (4.26).  Application  of  a  Taylor  series  expansion  and  retention  of  terms  up 
to  second  order  yields 


p 1  +  /(*>  o  =  p11  +  /(*»  0  -Jj- 


on  z  =  0 


Use  of  Eq.  (4.22)  in  the  above  relation  leads  to  the  desired  form  of  the  pressure 
balance  boundary  condition  on  z  =  0  that  is  correct  to  second  order, 


P 


ir 


on  z  —  0 


(4.27) 


In  obtaining  Eq.  (4.27),  the  first-order  form  of  the  particle  velocity  balance 
has  been  used  in  the  second-order  term.  Dropping  the  second-order  terms  in 
Eq.  (4.27)  leads  to 

p1  =  pn  on  z  =  0  .  (4.28) 

This  is  the  well-known  small-signal  form  of  the  pressure  boundary  condition:  the 
pressures  balance  at  the  interface  between  two  fluids. 


CHAPTER  5 


AN  ANALYSIS  OF  THE  REFLECTION  AND  REFRACTION  OF 
FINITE-AMPLITUDE  PLANE  WAVES  AT  A  PLANE 
FLUID-FLUID  INTERFACE  USING  SECOND-ORDER 
PERTURBATION  THEORY 


5-1  Introduction 

In  this  chapter  we  analyze  the  reflection  and  refraction  of  finite-am¬ 
plitude  plane  waves  that  are  obliquely  incident  on  an  initially  plane  fluid-fluid 
interface.  A  second-order  perturbation  analysis  method,  sometimes  referred  to 
as  quasilinear  analysis,  is  used.  The  notation  and  method  closely  follow  the 
work  of  Naze  Tj0tta  and  Tj0tta  (1987).  The  boundary  condition  at  the  source 
is  arbitrary.  Special  attention  is  given  to  the  0(e2)  source  boundary  condition 
and  the  motion  of  the  interface.  The  chapter  is  divided  as  follows:  In  the  first 
section,  each  acoustic  variable  is  expanded  in  a  power  series  in  e,  where  e  is  a  small 
parameter  (the  peak  particle  velocity  normalized  by  Co).  The  (9(e)  and  (9(e2) 
interface  boundary  conditions,  wave  equations  for  fluids  I  and  II,  and  relations 
between  the  various  acoustic  field  variables  are  obtained.  In  the  second  section, 
the  0(e)  system  is  solved,  and  expressions  for  the  0(e)  reflected  and  transmitted 
acoustic  fields  are  found.  In  the  third  section,  we  solve  the  0(e2)  system  and 
find  expressions  for  the  0(e2)  reflected  and  transmitted  acoustic  fields.  Although 
some  interpretation  of  the  results  is  given  as  they  are  obtained,  the  majority  of 
interpretation  of  results  is  deferred  to  the  next  chapter. 

It  is  useful  to  note  at  this  time  that  oblique  incidence  reflection  and  re¬ 
fraction  of  plane  waves  is  inherently  a  two-dimensional  problem.  Without  loss  of 
generality,  we  therefore  take  the  incident  ray  to  lie  in  the  x,z  plane.  The  subse¬ 
quent  motion — reflected  wave,  transmitted  wave,  and  motion  of  the  interface — is 
independent  of  the  coordinate  y.  Thus,  it  is  appropriate  for  our  work  to  define 
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the  propagation  direction  of  the  incoming  signal  as 


ninc  =  i  sin  $>DC  +  k  cos  0mc 


(5.1) 


where,  as  above,  0  is  measured  from  the  +z-axis;  see  Fig.  5.1.  (Recall  that  the 
interface  is  assumed  to  initially  lie  in  the  z  =  0  plane.)  This  orientation  of  the 
incident  signal  to  the  coordinate  system  is  assumed  through  the  remainder  of  this 
work.  Moreover,  all  angles  are  measured  from  the  -f  z-axis.  The  superscript  inc 
is  used  throughout  this  work  to  indicate  an  acoustic  variable  associated  with  the 
incident  signal;  similarly,  the  superscripts  refl  and  trans  are  used  to  indicate  the 
reflected  and  transmitted  signals. 
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Figure  5.1  Diagram  showing  coordinates  that  are  required  for  our  analysis 


5-2  Basic  Equations,  Relations,  and  Boundary  Conditions 

Obtained  in  this  section  are  0(e)  and  0(t2)  forms  of  the  following:  (1) 
the  wave  equations  for  fluids  I  and  II,  (2)  the  relations  between  the  modified 
acoustic  potential  $  and  some  of  the  other  field  variables,  namely  P\  <f>,  p\  and 
u,  and  (3)  the  interface  and  source  boundary  conditions. 
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Wave  equations  of  0(e)  and  0(e2) 

The  lossless  version  of  the  wave  equation,  Eq.  (3.25),  is 


(5.2) 


To  obtain  0(e)  and  0(e2)  forms  of  Eq.  (5.2),  we  expand  the  modified  velocity 
potential  $  in  a  power  series  in  terms  of  a  small  parameter  e  (which  is  defined 
later  from  the  source  boundary  condition), 


$  =  e$(1)  +  e2${2)  +  •  -  •  .  (5.3) 

If  the  expansion  is  inserted  into  Eq.  (5.2)  and  if  coefficients  of  like  powers  of  e  are 
equated,  then  the  following  0(c)  and  0(e2)  forms  of  the  wave  equation  for  fluid 
I  are  obtained: 


□2$m  =  0 


'(i) 


°%)  = 


d 

(4 )4  dt 


(5.4) 

(5.5) 


Recall  that  the  superscript  I  is  used  to  indicate  an  acoustic  variable  or  a  fluid 
constant  that  is  associated  with  fluid  I.  The  corresponding  pair  of  equations  that 
govern  the  0(e)  and  0(e2)  acoustic  fields  in  fluid  II  are 


□  2*") 


=  0 


P11  d  /a»;y 

(41)4  dt  \  dt  t 


(5.6) 

(5.7) 


0(e)  and  0(e2)  relations  between  $  and  P',  <f>,  p',  and  u 

Since  the  wave  equation  is  expressed  in  terms  of  the  modified  acoustic 
potential  $,  the  0(e)  and  0(e2)  forms  of  the  relations  between  $  and  some  of  the 
other  field  variables,  namely,  P',  (ft,  p',  and  u,  are  required.  We  start  by  expanding 
the  modified  acoustic  pressure,  the  velocity  potential,  and  the  acoustic  pressure 
in  power  series  in  e: 


P'  —  t^(i)  +  t2 P(i)  +  ’  ’  ’ 

4>  =  e^(i)  -f  t2 <f)(?)  +  ••• 
P'  =  eP(l)  +  f2p(  2)  +  •  •  • 


(5.8) 

(5.9) 

(5.10) 
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Note  that  for  convenience  the  primes  in  the  expanded  forms  of  P'  and  p'  have 
been  dropped.  The  particle  velocity  u  may  be  expanded  in  a  power  series  as 
follows: 

u  =  eu(!)  +  c2U(2)  +  •  •  •  •  (5.11) 

The  0(e)  and  0(e2)  forms  of  u  may  be  written  in  terms  of  their  x  and  ^  compo¬ 
nents, 


U(i)  =  iu(i)  +  kio(i)  ,  (5-12) 

U{2)  =  iu(2)  +  ku;(2)  .  (5.13) 

Alternatively,  the  components  of  u  may  be  directly  expanded  as  follows: 

u  =  eu(1)  +  e2«(2)  -f  •  •  •  ,  (5.14) 

w  =  et«(i)  +  e2tu(2)  +  •  •  •  .  (5.15) 

For  the  case  of  propagation  in  a  lossless  fluid,  the  expression  for  the 
modified  acoustic  pressure  P'  given  in  Eq.  (3.24)  simplifies  to 


P'  =  -Po 


dt 


(5.16) 


Inserting  the  expansions  for  P'  and  $  into  Eq.  (5.16)  and  then  equating  coeffi¬ 
cients  of  like  powers  of  e  yields  the  following  0(e)  and  0(e2)  relations  for  P'  in 
fluids  I  and  II: 


pi  .... 

Mi)  “ 

Po  dt  ’ 

(5.17) 

pi  _ 

M2)  “ 

p°  dt  ’ 

(5.18) 

pH  _ 

Mi)  “ 

n  a$|i) 

h  gt  . 

(5.19) 

pH  ... 

Mm  - 

P°  dt  ■ 

(5.20) 

The  modified  acoustic  potential  $  is  defined  in  Eq.  (3.23)  in  terms  of 
the  acoustic  potential  <f>.  Using  the  expansions  for  $  and  <f>  leads  to  the  0(e)  and 
0(e2)  forms  of  the  relation  between  $  and  for  fluids  I  and  II, 

*0)  =  ^(i)  ’ 


(5.21) 
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*(2)  ~  ^<2>  2(4)2  dt  ’ 

(5.22) 

(5.23) 

^(2)  ~  ^  2(c“)2  dt  • 

(5.24) 

Similarly,  the  0(e)  and  0(e 2)  forms  of  the  relation  between  the  modified 
velocity  potential  $  and  the  acoustic  pressure  p'  for  both  fluids  I  and  II  may  be 
found.  Using  the  expansions  for  p',  P\  and  <f>  in  the  defining  relation  for  P', 
Eq.  (3.21),  and  then  rearranging  and  substituting  using  Eqs.  (5.1 7)— (5.20),  we 
obtain 


P(2)  = 


p!!) 


I  d  !  Po 
-Po  g-(  *(2)  -  J 

II  ^  ^II 

~Podi*w  ’ 


(5.25) 

(5.26) 

(5.27) 

(5.28) 


Use  of  Eq.  (5.11)  in  Eq.  (3.9)  and  substitution  of  Eqs.  (5.21 )— (5.24)  into 
the  results  yields  the  0(e)  and  0(e2)  forms  of  the  relation  between  the  modified 
velocity  potential  $  and  the  particle  velocity  u  for  both  fluids  I  and  II: 


“(i)  =  V*!,,  ,  (5.29) 

=  +  -  (5-30) 
u(!)  =  ,  (5.31) 

u8i  =  Jt  •  <5'32> 


0(e)  and  0(e2)  forms  of  the  interface  boundary  conditions 

The  0(e)  and  0(e2)  forms  of  the  two  dimensional  versions  of  the  interface 
boundary  conditions,  Eqs.  (4.23)  and  (4.27),  may  be  obtained  by  substituting  in 
the  expansions  for  p',  u,  and  w.  Equating  coefficients  of  like  powers  of  e  yields  the 
0(e)  arid  0(e2)  forms  of  the  expression  for  the  normal  particle  velocity  balance 
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at  the  interface, 


4)  ,=0  =  w(')\,=0  ’ 


(5-33) 


w'£Lo-whLo=  (u")-“(i  >)/ 


dw11 

dx 


dt 


+ 


i) 


2=0 


dz 


dz 


IV. 


W)dt 


and  the  0(e)  and  0(e2)  forms  of  the  pressure  balance, 

p|')Lo  =  p"»L  ■ 
*l,=o  -  pi»Lo  =  “  lh)  / 


z=u 

(5.34) 


(5.35) 

(5.36) 


2= 0 


0(e)  and  0(e2)  forms  of  the  source  boundary  condition 

The  boundary  condition  at  the  source  is  arbitrary,  and  our  results  may 
therefore  be  used  to  analyze  transients  or  harmonic  motion.  Although,  strictly 
speaking,  our  source  must  be  infinite  in  extent  because  we  have  assumed  infinite 
plane  waves,  we  assume  that  the  reflected  signal  does  not  interact  with  the  active 
source.  This  is  reasonable,  however,  because  practical  experiments  that  approx¬ 
imate  plane  waves  can  be  devised.  Moreover,  the  theory  developed  here  may  be 
used  as  a  single  component  of  the  spatial  Fourier  decomposition  of  a  directive 
source,  in  which  case  the  geometry  may  be  selected  such  that  the  reflected  signal 
does  not  return  to  the  active  source. 

The  boundary  condition  at  the  source  may  be  written  as  follows: 

n  d 

u  =  — r  ■37‘s’(*  -  To)  at  r  =  r0  ,  (5.37) 

Cq  ot 

where  r  =  r0  defines  a  plane,  5  is  any  given  function  of  t,  and  r0  is  defined  as 
n  *r0/co.  The  coefficient,  the  time  derivative,  and  the  phase  shift  are  introduced 
for  later  convenience,  specifically,  to  make  it  easier  to  match  the  solution  in  terms 
of  *he  modified  velocity  potential.  Using  tfi*»  expansion  of  the  particle  velocity 
given  in  Eq.  (5.11)  and  expanding  5  in  a  power  series  as 

S  =  cS(i)  +  e2^{2)  +  *  •  •  (5.3S) 

leads  to  the  0(e)  and  0(t2)  forms  of  the  arbitrary  particle  velocity  source  condi¬ 
tion, 

U(1)  =  "5i5(1)(<~r°)  at  r  =  r°  ’ 


(5.3!)) 
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n  d 

u(2)  =  -^r  -  To)  at 


(5.40) 


5-3  Solution  of  the  0(e)  System 

In  this  section,  the  0(e)  system  is  solved.  The  0(e)  system  consists  of 
the  0(e)  wave  equations  for  fluids  I  and  II,  Eqs.  (5.4)  and  (5.6),  and  the  0(e) 
boundary  conditions.  The  0(e)  boundary  conditions  consists  of  the  0(e)  interface 
boundary  conditions,  Eqs.  (5.33)  and  (5.35),  and  the  0(e)  boundary  condition  at 
the  source,  Eq.  (5.39).  To  these  we  add  the  Sommerfeld  radiation  conditions  for 
fluids  I  and  II,  that  is,  that  the  reflected  and  transmitted  waves  propagate  away 
from  the  interface  in  fluids  I  and  II,  respectively.  First,  the  general  solutions  of 
the  wave  equations  in  fluids  I  and  II  are  found.  Next,  application  of  the  interface 
boundary  condition  leads  first  to  the  law  of  specular  reflection  and  Snell’s  law 
as  conditions  for  the  validity  of  the  solutions  and  second  to  the  reflection  and 
transmission  coefficients.  Finally,  the  incident  signal  is  matched  to  the  source 
boundary  condition. 

To  solve  the  0(c)  wave  equations,  we  first  note  that  because  they  are 
homogenous,  no  particular  solutions  are  required.  If  the  reflected  and  transmitted 
waves  are  assumed  to  be  planar,  the  general  solution  of  the  0(c)  wave  equations 
for  fluids  I  and  II  may  be  written  as  follows: 


(5.41) 

*")(■•.  o = *sr(r(‘5“)  • 

(5.42) 

The  definitions  of  rmc,  r^n ,  and  r{™ns  are 

inc  —  .  n'nC  *r 

4  ’ 

(5.43) 

refl  _  t  n(i!  *r 

r(i)  — 1  gi  ’ 

(5.44) 

,  nftr*  -r 

_trans  _  t  U) 

T(l)  ~  1  _I1 

(5.45) 

where  and  n^8  are  unit  vectors  in  the  direction  of  propagation  of  the 
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reflected  and  transmitted  waves,  respectively  (see  Fig.  5.2),1 

=  i  sin^'®  +  k  cos  0"^  where  cos  0^  <  0  ,  (5.16) 

n\\ ^  =  i  sin  0\\?s  +  k  cos  where  cos  O^f8  >  0  .  (5,1 7) 

The  assumption  that  the  reflected  and  transmitted  waves  are  planar  is  tested 


Figure  5.2  Diagram  showing  the  angle  of  incidence  and  the  0(e)  angles  of  re¬ 
flection  AND  TRANSMISSION 

by  checking  whether  the  assumed  solution  satisfies  the  interface  boundary  condi¬ 
tions. 

Use  of  Eqs.  (5.41)  and  (5.42)  in  the  0(e)  boundary  conditions  at  the 
interface,  Eqs.  (5.33)  and  (5.35),  respectively,  leads  to 


cos  0"J 

cos  0^™ 

cos  0inc  3$^ 

<0  dl  ,=0 

4'  dt  z=0 

Co  dt 

'This  nonstandard  definition  of  (lie  reflection  angle  0refl  is  chosen  to  simplify  our  later  work. 
It  turns  out  to  he  easier  to  manage  the  problem  if  all  angles  (more  are  introduced  in  the  O(i’) 
analysis)  are  defined  relative  to  the  same  reference;  we  chose  the  +;-axis. 
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and 


0*S) 

Pa  a*sr 

dt 

,=0  Po  dt 

CO 

1 

o 

II 

N 

where  we  have  used  the  following  chain  rule  result: 


(5.49) 


= 


n  9$(i) 
cq  dt 


(5.50) 


This  relation  is  used  repeatedly  in  later  analysis.  Equations  (5.48)  and  (5.49) 
are  written  with  the  0(e)  incident  signal  on  the  right-hand  side  because  it  is  a 
‘known’  function.  The  form  of  the  0(e)  incident  signal  is  obtained  by  matching 
the  0(e)  source  condition,  which  is  deferred  until  later. 

Snell’s  law  and  the  law  of  specular  reflection  are  now  obtained  as  condi¬ 
tions  for  the  validity  of  the  solutions.  Our  derivation  of  Snell’s  law  and  the  law 
of  specular  reflection  is  a  little  longer  than  some  other  more  common  approaches. 
This  approach  is,  however,  the  same  one  that  we  use  to  solve  the  0(e2)  problem. 
It  is  therefore  used  here  so  that  the  reader  may  become  familiar  with  it.  Equa¬ 
tions  (5.48)  and  (5.49)  represent  a  system  of  two  equations  and  two  unknowns 
with  two  conditions  to  be  met,  namely,  that  the  solutions  be  valid  for  all  x  and 
all  t.  The  conditions  may  be  obtained  by  taking  the  x  and  t  derivatives  of  the 
pair  of  equations,  solving  each  set  independently,  then  forcing  the  solutions  to  be 
the  same.  But  if  the  solutions  are  to  be  the  same,  then  the  equations  generated 
by  taking  the  x  and  t  derivatives  of  the  boundary  conditions  must  be  linearly 
dependent.  Accordingly,  we  arrive  at  the  following  conditions  by  forcing  the  x 
and  t  derivatives  of  Eqs.  (5.48)  and  (5.49)  to  be  linearly  dependent: 


sin  0inc  =  sinfl^J  , 
sin  0'nc  =  nsin^f8 


where  n  is  the  ratio  of  the  sound  speeds, 


n  = 


_ 
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(5.51) 

(5.52) 


(5.53) 


Equation  (5.52)  is  Snell’s  law,  and,  as  we  soon  see,  Eq.  (5.51)  leads  to  the  law  of 
specular  reflection. 

We  now  solve  for  the  angles  and  To  satisfy  Eqs.  (5.51)  and 

(5.46)  simultaneously,  0",c  and  must  be  related  as  follows: 


6 


refl 

(1) 


180°  -  0'mc 


(5.54) 
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To  satisfy  Eqs.  (5.52)  and  (5.47)  simultaneously,  and  6,nc  must  be  related 
by 

=  sin-1(n-1  sin  0mc)  .  (5.55) 

If  n  >  1,  is  real  for  all  0,nc\  if  n  <  1,  is  real  provided  that  |  sin  #,nc|  <  n. 
The  angle  of  incidence  at  which  |sin0mc|  =  n  is  referred  to  as  the  critical  angle. 
In  what  follows  we  assume  for  simplicity  that  is  real.  The  case  of  0|™'s 

being  complex  (incidence  above  critical  angle)  is  briefly  discussed  in  Sec.  5-5. 

Use  of  Eqs.  (5.54)  and  (5.55)  in  Eqs.  (5.43),  (5.44),  and  (5.45)  shows 
that  the  values  of  r,nc,  r^®,  and  are  equal  when  evaluated  at  the  interface. 

For  convenience,  a  new  term  that  is  equal  to  rinc,  r^f®,  and  when  they  are 

evaluated  at  the  interface  is  now  introduced, 


Tff  =  t  - 


xsin  0>nc 


(5.56) 


The  reflection  and  transmission  coefficients,  denoted  R  and  T,  respec¬ 
tively,  are  now  introduced,  and  the  two  interface  boundary  conditions  are  used  to 
yield  expressions  for  the  coefficients.  The  expressions  depend  only  on  the  physical 
properties  of  the  fluid  and  the  angle  of  incidence.  The  coefficients  R  and  T  relate 
the  reflected  and  transmitted  pressures2  at  the  interface  to  the  incident  pressure 
at  the  interface  and  are  defined  as  follows: 


refl 

P(l) 

5= 0 

TlinC 

P(i) 

5=0 

psr 

’I 

U=o 

pi") 

5  =  0 

(5.57) 


(5.58) 


where  p'^j,  p(j®,  and  P(^"s  are  the  pressures  corresponding  to  and 

respectively.  At  this  time,  we  also  define  the  acoustic  impedance  of  fluids  I  and 
II,  respectively, 


Zl 

Z 11 


_  Poco 


cos  9mc 

PqCq 
cos  0\\™ 


(5.59) 

(5.60) 


2 Some  authors  use  coefficients  based  on  the  acoustic  power  of  the  incident,  reflected,  and 
transmitted  waves.  For  clarity,  our  coefficients  are  sometimes  referred  to  as  amplitude  reflection 
and  transmission  coefficients. 


59 


Note  that  the  ratio  of  impedances  is,  accordingly, 

_  n  cos 
Z11  m  cos  0,nc 


(5.61) 


where  m  is  the  ratio  of  the  densities  (not  the  relaxation  variable,  which  is  no 
longer  required), 


-  Po 

m  =  — 7- 

Po 


(5.62) 


Note  that  m  is  defined  as  (fluid  II  density/fluid  I  density)  whereas  n  is  reciprocally 
defined  (fluid  I  sound  speed/fluid  II  sound  speed).  The  reason  for  the  difference 
is  merely  convenience.  Use  of  Eqs.  (5.57)  and  (5.58)  in  Eqs.  (5.48)  and  (5.49) 
yields 


and 

where  we  have  noted  that 


o  n  cos 

1  -R= - &—T 

m  cos  trnc 

l+R=T  , 
cos  0 (J®  =  —  cos  8mc 


Manipulations  of  Eqs.  (5.63)  and  (5.64)  yield  the  traditional  results 

Z11  -  Zl 


R  = 


and 


or,  in  terms  of  n  and  m, 


T  = 


Z 11  +  Z1 

2  Zn 

z1  +  zu 


R  = 

T  = 


m  cos  0mc  —  n  cos 
m  cos  0,nc  -(-  n  cos 
2m  cos  9mc 
m  cos  6mc  +  n  cos 


(5.63) 

(5.64) 

(5.65) 

(5.66) 

(5.67) 

(5.68) 

(5.69) 


We  now  point  out  some  interesting  special  cases.  Note  that  for  intromis¬ 
sion  (no  0(e)  reflection)  to  occur,  the  acoustic  impedances  of  the  two  fluids  must 
be  the  same,  that  is,  Z1  =  Z11,  or,  by  Eq.  (5.61), 


intromission  : 


m 


cosfftr*  =  —  cos  0" 
'  '  n 


(5.70) 
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By  squaring  both  sides  of  Eq.  (5.70),  eliminating  the  cosines  in  favor  of  sines, 
and  using  Snell’s  law,  we  see  that 


sin2  0intro  = 


m2  —  n2 


m 


2  -  1 


(5.71) 


where  0inlro  =  0inc 

(intromission 

following  condition  must  be  met: 


Note  that  for  a  intromission  angle  to  exist,  the 


m2  —  n2 


m2  —  1 


>  0 


Another  special  case  occurs  when  the  impedance  of  fluid  II  is  very  small  in  com¬ 
parison  to  that  of  fluid  I.  The  reflection  coefficient  tends  to  —1,  and  the  interface 
is  referred  to  as  a  pressure  release  interface.  If,  on  the  other  hand,  the  impedance 
of  fluid  II  is  very  large  in  comparison  to  that  of  fluid  I,  the  reflection  coefficient 
tends  to  unity.  In  this  case,  the  interface  is  referred  to  as  rigid.  Note  that  al¬ 
though  the  transmission  coefficient  tends  to  2  in  this  case,  it  can  be  shown  that 
the  acoustic  power  transmitted  tends  to  zero. 

We  now  obtain  the  reflected  and  transmitted  fields  in  terms  of  the  inci¬ 
dent  field.  Note  that  the  definitions  of  R  and  T  may  be  rearranged  to  yield 

and 

I* =  • 

"The  solutions  away  from  the  interface  are  obtained  by  replacing  the  independent 
variable  r9  with  rff  for  the  reflected  field  and  with  rffi"3  for  the  transmitted  field. 
The  resulting  field  relations  may  be  integrated  once  with  respect  to  time  noting 
that  the  integration  constant  must  be  zero  in  order  to  satisfy  the  condition  that 
the  reflected  and  transmitted  fields  be  zero  when  the  incident  field  is  zero.  Thus, 
away  from  the  interface,  the  reflected  and  transmitted  fields  are 

=  «f,!"5(r"!1)  .  (5.72) 

T 

< rans / _ Irans \  _  ,f,inc/_trans \  ir  - ■> \ 

We  now  obtain  an  expression  for  that  matches  the  source  boundary 
condition.  Because  the  reflected  wave  is  assumed  not  to  return  to  the  active 
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source,  the  reflected  wave  is  not  involved  in  matching  the  source  boundary  con¬ 
dition.  We  may  therefore  write 


ninc  dtifi  =  n  dS{l) 
4  dt  r=ro  4  & 


(5.74) 


Accordingly,  to  match  the  source  condition,  the  following  must  be  true: 


nmc  =  n  ,  (5.75) 

o)  •  (5.76) 

r=r0  ' 

Since,  when  evaluated  at  r  =  r0,  the  independent  variables  of  both  functions  are 
equal,  the  solution  away  from  the  interface  may  be  obtained  by  replacing  the 
independent  variable  with  rinc, 

To  obtain  an  expression  for  rather  than  its  time  derivative,  we  must  inte¬ 
grate  the  foregoing  relation  with  respect  to  time.  This  would,  in  general,  mean 
that  an  arbitrary  function  that  is  linearly  dependent  on  r,  say  F(r),  must  be 
determined.  We  know,  however,  that  the  boundary  condition  at  the  source  re¬ 
quires  that  VjF|r=ro  be  zero.  The  function  .F(r)  must,  therefore,  be  a  constant. 
Since  <!>(!)  is  not  a  measurable  quantity  and  all  measurable  quantities  depend  on 
derivatives  of  $(i),  we  are  free  to  choose  this  constant  such  that  $(])  at  the  source 
at  time  t  =  0  is  zero.  We  may,  therefore,  write  the  incident  field  as 


=  S(l|(r~)  -  S„,(To)  .  (5.77) 

Accordingly,  the  reflected  and  transmitted  fields  are 

*0)  =  K  -  Sm(T ■„))  ,  (5.78) 

*;;r  =  ~  (s<.>wn  -  s<.iM  •  (5.-9) 

To  summarize,  we  have  derived  the  following  major  results:  (1)  Snells 
law  and  the  law  of  specular  reflection,  Eqs.  (5.52)  and  (5.54),  respectively,  (2) 
the  expressions  for  the  reflection  and  transmission  coefficients,  Eqs.  (5.6G)  and 
(5.67),  respectively,  and  (3)  solutions  representing  the  0(e)  incident,  reflected, 
and  transmitted  fields,  Eqs.  (5.77),  (5.78),  and  (5.79),  respectively. 
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5-4  Solution  of  the  0(e2)  System 

In  this  section  the  0(e2)  system  is  solved.  The  0(e2)  system  consists  of 
the  0(e2)  wave  equations  for  fluids  I  and  II,  Eqs.  (5.5)  and  (5.7),  respectively,  the 
0(e2)  boundary  conditions  at  the  interface,  Eqs.  (5.34)  and  (5.36),  respectively, 
the  Sommerfeld  radiation  c  mditions  for  fluids  I  and  II,  and  the  0(e2)  boundary 
condition  at  the  source,  Eq.  (5.40).  The  procedure  is  as  follows:  First,  t he 
general  solutions  for  fluids  I  and  II  are  separated  into  two  parts — the  particular 
solution,  denoted  by  the  subscript  p,  and  the  homogeneous  solution,  denoted  by 
the  subscript  h: 


*[i)  =  *|i)p  +  *!■„;.  .  (5.80) 

*ii>  =  *(».  + *P)»  '  (5  81) 

The  particular  solutions  are  chosen  to  satisfy  the  0{ t2)  wave  equations.  The 
general  solutions  are  then  substituted  into  the  0(e2)  interface  boundary  condi¬ 
tions,  and  the  0(e2)  homogeneous  solutions  are  chosen  to  satisfy  these  boundary 
conditions.  Terms  accounting  for  the  finite  displacement  of  the  interface  and 
the  variation  of  the  normal  to  the  interface  are  identified.  The  results  are  then 
summarized. 


The  0(e2)  particular  solution 


Since  the  0(c)  solution  is  known,  it  may  be  substituted  into  the  right- 
hand  side  of  the  0(e2)  wave  equations.  However,  since  the  final  form  of  the  0(e) 
solution  depends  on  the  boundary  condition  at  the  source,  the  functions  4)",’^. 

and  ${jy's  are  used.  The  0(e2)  wave  equations,  Eqs.  (5.5)  and  (5.7),  then 
become 


□a(*i»,+*!*>  -  jjj  | 


.  dt 


21 


+  2 

J*(i) 

dt 

dt 

0U 

d  (d<t>\\™ 

(4'Y 

dt  \ 

dt 

+ 


V  dt  ) 


,  (5.82) 

(5.83) 


These  equations  are  significantly  more  complicated  than  their  0(e)  counterparts. 
The  nonlinear  terms  on  the  right-hand  side  of  Eq.  (5.82)  represent,  respectively, 
the  self-action  of  the  incident  signal,  the  interaction  between  the  incident  and 
reflected  signals,  and  the  self-action  of  the  reflected  signal.  The  nonlinear  term 
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on  the  right-hand  side  of  Eq.  (5.83)  represents  the  self-action  of  the  transmitted 
signal. 

Because  the  wave  equation  for  fluid  II  is  simpler,  we  solve  for  its  partic¬ 
ular  solution  first.  Note  that  the  right-hand  side  of  Eq.  (5.83)  may  be  expressed 
as  a  d’Alembertian, 


0n  d  (d$\\]nt 


(41)4  dt  V  dt 


cr 


P" 


2  (41) 


/  „  trans 

II \3  ("(I)  ^ 


(i) 


dt 


(5.84) 


Equation  (5.83)  may  thus  be  written 


a4 


+ *&>» + j (fp  ("ir  *>  ( 


d<I>^ns\2 

“ ^T“  j 


=  o 


(5.85) 


The  particular  solution  for  fluid  II  is  now  obvious, 

'2(c» ) 


V(2)P 


r  i  trans  [  I1) 

9t/.lh3  ( n(l)  *r)  (  Qt 


(5.86) 


The  remaining  equation, 

□  2^(2)/i  —  0  ,  (5.87) 

is  to  be  solved  subject  to  the  Sommerfeld  radiation  condition  for  fluid  II  and  the 
0(e2)  interface  boundary  conditions.  We  defer  this  to  later.  Note  that  the  partic¬ 
ular  solution  fur  fluid  II  exhibits  amplitude  growth  in  the  direction  of  propagation 
via  the  coefficient  (n^p8  -r).  Moreover,  note  that  the  direction  of  propagation  of 
the  particular  solution  for  fluid  II  is  given  by  Snell’s  law. 

To  obtain  the  particular  solution  of  the  wave  equation  for  fluid  I,  we  fol¬ 
low  the  method  of  Naze  Tjptta  and  Tjptta  (1987).  As  noted  earlier,  the  particular 
solution  may  be  thought  of  as  being  composed  of  three  parts:  the  self-action  of  the 
incident  wave,  the  self-action  of  the  reflected  wave,  and  the  interaction  between 
the  incident  and  reflected  waves.  The  fluid  I  wave  equation  is  first  rearranged  to 
find  the  particular  solution  that  pertains  to  the  interaction  between  the  incident 
and  reflected  waves.  Noting  that  □2(^|"p2  =  0  and  that  IH2($^))2  =  0  leads  to 
the  following  form  of  □2(<I>|1))2: 


□2(4>i„)2  =  2d2(>i>;";  <!>;;?)  = 


2 

( c& )2 


•n 


refl 

(I) 


-1) 


dt  dt 


(5.88) 
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Use  of  the  definitions  of  nmc  and  n(*p  the  law  of  specular  reflection,  and  some 
algebra  simplifies  Eq.  (5.88)  to 


d$(i) 
dt  dt 


(4)2 

4  cos2  &'nc 


(5.S9) 


Note  that  Eq.  (5.89)  may  be  expressed  in  terms  of  the  Lagrangian  £,  which  was 
defined  in  Eq.  (2.83).  The  0(e)  values  of  the  Lagrangian  in  both  fluids  I  and  II 
are  zero.  The  0(e2)  value  of  the  Lagrangian  in  fluid  II  is  also  zero.  However,  the 
0(e2)  value  of  the  Lagrangian  in  fluid  I  is 


7°s 


(5.90) 


With  Eq.  (5.89),  the  particular  solution  that  pertains  to  the  interaction  be¬ 
tween  the  incident  and  reflected  waves  is  at  hand.  The  inhomogeneous  terms  in 
Eq.  (5.82)  that  correspond  to  the  self-action  of  the  incident  and  reflected  waves 
are  of  the  same  form  as  the  inhomogeneous  term  in  the  fluid  II  wave  equation. 
Thus,  they  too  may  be  expressed  as  d’Alembertians  that  are  similar  in  form  to 
Eq.  (5.84).  Using  this  idea  and  Eq.  (5.89)  in  Eq.  (5.82)  leads  to  the  following: 


(5.91) 


The  0(c2)  particular  solution  for  fluid  I  is  therefore 


<f> 


(2  )p 


2(4)3 


(ninc  *r) 


+ 


dt 


+  ("<?>  -r) 


reR 


cos2  0tnc  dt 


dt 


(5.92) 


The  remaining  equation, 

□  %)a  =  0  ,  (5.93) 

must  be  solved  subject  to  the  0{t2)  interface  boundary  conditions  and  the  0(<2) 
source  condition.  Note  that  the  two  terms  in  the  particular  solution  for  fluid 


I  that  correspond  to  the  self-action  of  the  incident  and  reflected  waves  exhibit 
amplitude  growth  in  the  direction  of  propagation.  The  amplitude  growth  is  given 
by  the  coefficients  (nmc  *r)  and  (n"®  *r).  Moreover,  note  that  the  direction  of 
propagation  of  the  term  that  corresponds  to  the  self-action  of  the  reflected  wave 
is  given  the  law  of  specular  reflection. 

The  0(e2)  homogeneous  solution 

Now  that  the  particular  solutions  to  the  <9(e2)  wave  equations  are  known, 
we  may  find  the  0(e2)  homogeneous  solutions.  We  start  by  examining  the  0(e2) 
interface  boundary  conditions  and  introducing  some  new  notation.  The  general 
solutions  for  fluids  I  and  II  are  then  substituted  into  the  0(e2)  interface  bound¬ 
ary  conditions,  and  two  relations  between  the  0(e2)  homogenous  solutions  are 
obtained.  Next,  the  relations  are  solved  for  the  0(e2)  homogeneous  solutions  for 
both  fluids  I  and  II. 

A  noteworthy  difference  between  the  0(e2)  interface  boundary  conditions 
and  their  0(e)  counterparts  is  that,  as  a  result  of  the  particular  solution  being 
evaluated  at  the  interface,  the  0(e2)  interface  boundary  conditions  contain  terms 
that  grow  with  x.  This  makes  sense  physically  because  the  amplitude  of  the 
particular  solution  grows  with  the  distance  propagated.  Since  the  particular 
solution  strikes  the  interface  at  oblique  incidence,  the  distance  propagated  varies 
with  position  along  the  interface;  consequently,  its  amplitude  varies  along  the 
interface.  It  thus  turns  out  that  the  0(e2)  homogeneous  solutions  is  composed  of 
two  parts:  The  part  that  accounts  for  the  aforementioned  particular  solution  has 
amplitudes  which  vary  along  its  wavefront.  The  part  that  accounts  for  everything 
else  has  uniform  amplitude  along  its  wavefront.  The  propagation  direction  of 
both  waves  is  assumed  constant.  The  assumption  of  constant  angle  is  difficult  to 
maintain  for  all  x  and  all  t  if  the  wave  with  the  amplitude  variation  along  the 
wavefront  is  overlooked.3 

The  introduction  of  some  new  notation  makes  matching  the  0(e2)  in¬ 
terface  boundary  conditions  more  manageable.  Recall  that  the  0(e2)  interface 
boundary  conditions,  given  in  Eqs.  (5.34)  and  (5.36),  contain  terms  that  account 
for  the  displacement  of  the  interface,  the  variation  of  the  normal  to  the  interface, 
and  the  0(e2)  pressure  and  particle  velocity  differences.  Singling  out  the  0(t2) 
pressure  and  particle  velocity  differences  for  a  moment,  we  see  that  they  may  be 
expanded  into  their  homogeneous  and  particular  parts  and  written  in  terms  of 

3A  wave  with  amplitude  variation  along  the  wavefront  is  a  solution  of  the  wave  equation.  For 
example,  p  =  A(y)f(x  -  co<)  is  a  solution  of  the  one-dimensional  wave  equation,  plt  —  c2pTT  =  0, 
where  the  subscripts  indicate  differentiation. 
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the  modified  velocity  potential  $, 


(5.94) 


«*S)p 


1  d  d  *  2 

+  2 


1  d  9  (frH  )2 

2 [$ydtdz{  (1)j 


(.5.95) 


Use  of  Eqs.  (5.94)  and  (5.95)  in  Eqs.  (5.34)  and  (5.36)  leads  to  some  lengthy 
expressions.  We  therefore  introduce  the  following  shorthand  notation  for  not 
only  the  terms  in  Eqs.  (5.94)  and  (5.95),  but  also  for  the  terms  in  Eqs.  (5.34) 
and  (5.36)  that  represent  the  variation  of  the  normal  and  the  displacement  of  the 
interface: 


A  accounts  for  the  particular  solution  of  the  particle  velocity  and  is  subdivided 
into  two  components,  Ax  and  Arcf,  where  Ax  represents  those  terms  that, 
when  evaluated  at  the  interface,  have  x  as  a  coefficient,  and  /lrcg  represents 
the  ‘regular  terms’  that  do  not  have  x  as  a  coefficient, 

B  accounts  for  the  nonlinear  relation  between  the  particle  velocity  and  the 
modified  velocity  potential, 

C  accounts  for  the  particular  solution  of  the  pressure  and  is  also  subdivided 
into  two  components,  Cx  and  Creg,  for  the  same  reasons  as  for  Ax  and  .4reg, 

I)  accounts  for  the  nonlinear  relation  between  the  pressure  and  the  modified 
velocity  potential, 

E  accounts  for  the  displacement  of  the  interface  in  the  particle  velocity  bound¬ 
ary  condition, 

E  accounts  for  the  displacement  of  the  interface  in  the  pressure  boundary 
eondit  ion, 
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G  accounts  for  variation  of  the  normal  in  the  particle  velocity  boundary  con¬ 
dition. 


Mathematically,  the  shorthand  notation  is  defined  as  follows: 


A^A~*+xA^-irL--ir,=0  ' 

R  =  1  d  12 _ L_A  L  1 

-  2(4)2  dt  dz{  2(c”)2&  dz{  <l)) 

d*L„  d$JL  „ 

C  =  Creg  +  xCx  =  -  m 

dt  dt  __n 


1  d  d 
0  ~  2(c”)2  dt  dz 


(5.96) 


(5.97) 


(5.98) 


D  s  I  (v‘ +  (3)4)  ~t(v2  +  (w£)  (*">>2  „„  ’ 


^  dw\l)  gw}j)\ 

^  $2  92  ) 

1  fdpl{l)  dp Ji} 


Po  \  dz  dz 


to”)  dt 


to”)  dt 


c  =  (“”) -“<.»)  / 


(5.99) 

(5.100) 


(5.101) 


(5.102) 


Use  of  the  foregoing  definitions  in  the  interface  boundary  conditions  yields  the 
following  more  compact  form  of  the  interface  boundary  conditions: 


=  ~(C„.  +  xC,  +  D)  +  F  .  (5.103) 


—  — (/4reg  +  xAx  -f  B  +  E  +  G)  .  (5.104) 


Explicitly  shown  in  Eqs.  (5.103)  and  (5.104)  is  the  interrelation  between  the  0(c2) 
homogeneous  solutions  for  fluids  I  and  II. 

An  examination  of  the  definitions  of  A  through  G  reveals  that  they  are 
quadratic  forms  of  O(c)  terms.  Moreover,  A  through  G  are  defined  at  2  =  0,  and 
relations  between  the  O(e)  terms  at  2  =  0  are  also  known.  It  is  thus  possible  to 
specify  A  through  G  in  terms  of  the  O(e)  incident  signal,  the  properties  of  fluid  I, 
and  the  nondimensional  ratios  of  the  fluid  properties.  The  general  procedure 
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is  as  follows.  First,  we  use  Eqs.  (5.72)  and  (5.73)  to  express  the  functions 
and  <Flj^ns  in  terms  of  Next,  Snell’s  law  and  the  law  of  specular  reflection, 
Eqs.  (5.52)  and  (5.54),  are  used  to  express  the  coefficients  sin  Orff  and  sin  in 
terms  of  sin0,nc.  Using  Eq.  (5.65),  we  rewrite  the  coefficient  costfj^  in  terms  of 
cos</mc;  the  coefficient  cos#^'8,  which  usually  appears  with  the  coefficient  T,  is 
eliminated  in  favor  of  cos0mc  by  way  of  Eq.  (5.63).  In  general,  we  eliminate  1  in 
favor  of  R  using  Eqs.  (5.63)  and  (5.64),  or  a  combination  thereof.  Also  note  that 
the  time  and  space  derivatives  of  $(q  are  related  by  Eq.  (5.50).  The  resulting 
expressions  for  A  through  G  (with  the  ‘x-dependent’  parts  of  A  and  C  separated 
from  the  ‘regular’  parts)  are 


/l|  eg(  T~0  ) 

•rAI(re) 

B(t») 

f-'reg(  Tq  ) 

*Cx(to) 


~(co)2  4 

cos  0 


m )  \  dt  ) 


(5.105) 


2  =  0 


f,  /?"nq  xsinfl'"'  0 
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_ 1  cosO**  ,  ( 

2M)2  cl  (  B)( 

_ f. 
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m) 


-’=0 


Z  r.O 


2(4) 
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x sin 0“,c  d  /^joV 


dt  \  dt  ) 


D{t9)  = 


(4) 


I  \2 


(5.106) 

(5.107) 

(5.108) 


?=o 

(5.109) 


2  =  0 

(5. 110) 
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2  =  0 


I'{to)  = 


2  0"'c  „ 
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(4)2 


-(1  -  RY{m  -  1)<& 


(,)  dt2 


Ui) 


,  (5.111) 


(5.112) 
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G(t6)  = 


cos  Qmc  sin2  0mc  2  (  m  —  1  \  /  d4>(")  \ 

{~8T) 


(5.113) 


Equations  (5.105)— (5.113)  are  a  significant  result  because,  as  we  soon 
find  out,  the  terms  A  through  G  appear  throughout  the  homogeneous  solutions. 
Obtaining  simplified  forms  of  the  homogenous  solutions  for  special  cases  such  as 
reflection  from  a  rigid  wall  (R  =  1)  or  reflection  from  a  pressure  release  surface 
(/?  =  —  1)  is  now  relatively  straightforward.  (In  the  next  chapter,  we  analyze  the 
special  case  of  intromission,  R  =  0.)  Moreover,  the  dependence  of  the  terms  on 
the  changes  in  physical  properties  of  the  fluids  is  now  apparent,  specifically, 


/ 31  m 


m 

m  —  1 
m 


0"  r  po<j 

pg(4V  l  0 

'  =  ’ 


(5.114) 

(5.115) 

(5.116) 


where  the  symbol  |  ]  indicates  the  jump  of  the  enclosed  quantity  at  the  inter¬ 
face,  that  is,  (property11  —  property1).  Equation  (5.114)  represents  the  jump  in 
a  new  property,  the  ratio  of  the  compressibility  to  the  coefficient  of  nonlinearity, 
whereas  Eqs  (5.115)  and  (5.116)  represent  the  jump  in  the  compressibility  and 
the  density,  respectively.  Note  that  the  terms  that  represent  the  motion  of  the 
interface  vanish  if  both  the  densities  and  small-signal  sound  speeds  are  matched. 
Moreover,  it  turns  out  that  the  general  condition  for  no  reflections  (perfect  trans¬ 
mission)  is  that  the  densities,  the  small-signal  sound  speeds,  and  the  coefficients 
of  nonlinearity  be  matched. 

We  now  solve  the  0(e2)  homogeneous  wave  equations,  Eqs.  (5.87)  and 
(5.93),  subject  to  the  0(e2)  interface  boundary  conditions  as  given  in  Eqs.  (5.103) 
and  (5.104).  We  earlier  noted  that  the  0(e2)  interface  boundary  conditions  con¬ 
tain  terms  that  grow  with  x.  Futhermore,  we  noted  that  these  terms  are  ac¬ 
counted  for  by  a  separate  function — one  that  exhibits  an  x-dependence  when 
evaluated  at  the  interface.  This  function,  which  is  denoted  by  a  subscript  a,  and 
a  constant  amplitude  function,  which  is  denoted  by  the  subscript  5,  satisfy  the 
interface  boundary  conditions.  The  reflected  and  transmitted  signals  are  sought 
in  the  following  form: 


*Un(M)  =  ♦&<’•“>  +  (n&i.  "•)*$,(>•$)  +  ♦SSUrSj) 


(5117) 


70 


*(V(r.  0  =  K™  •r)*ffi,(rSD  +  *ST(TSr)  , 

where  the  amplitude  variation  along  the  interface  is  given  by 

(n(2)oi  ‘r)  a,,d  where 
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n(2)?S  *r 

T(2)6  “ 

41 

n[2)0  =  i  sin  0T$a  +  k  cos  0‘(2)a  , 

n$6  =  isin0$6  +  kcos0$)  , 
ngr  =  isin^a  +  kcos^s  , 
ngjj*  =  i  sin  0g£*  +  k  cos  , 

n(2)aX  =  >  sin  0(2 )oX  +  k  COS  0(2 Jal  , 

nSS.  =  isin^ax  +  ksin0^aX  . 


(5.118) 

(nr2)Ux -r) a,ul 

(5.119) 

(5.120) 

(5.121) 

(5.122) 

(5.123) 

(5.124) 

(5.125) 
(5.120) 

(5.127) 

(5.128) 


That  (n^'jax  *r)^(2)a(T(2)1a)  Is  a  solution  of  the  homogeneous  wave  equation  may 
be  verified  by  direct  substitution  if  it  is  noted  that  the  vectors  n(2)aX  and  11  Uhl 
are  perpendicular  to  the  direction  of  propagation,  that  is, 


n(2)o  *  n(2)oX  =  0  ,  (5-129) 

n(™'n(™l  =  °  •  (5-i;Jo) 


However,  two  possibilities  exist  for  both  n^x  and  nj!^'-  ^  rnay  have  a  com¬ 
ponent  in  either  the  +i  direction  or  the  -i  direction.  We  choose  both 
and  n(2f)ax  to  have  components  in  the  -fi  direction;  see  Fig.  5.3.  Direct  substi¬ 
tution  of  Eqs.  (5.117)  and  (5.118)  into  the  0(e2)  interface  boundary  conditions, 
F,qs.  (5.103)  and  (5.104),  yields  what  appears  to  be  a  system  of  two  equations  and 
four  unknowns  with  four  angle's  to  be  determined  subject  to  the  two  conditions 
that,  the  solution  be  valid  for  all  x  and  all  t.  Fortunately,  we  may  separate  the 
terms  that  exhibit  growth  with  x  and  demand  that  they  equate  separately.  We 


Figure  5.3  Diagram  showing  direction  of 


and  n 


Ho 
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are  thus  left  with  two  systems  of  two  equations  and  two  unknowns  with  two  an¬ 
gles  to  be  determined  subject  to  the  two  conditions  that  the  solutions  be  valid  for 
all  x  and  all  t — one  system  for  the  ‘x-dependent  terms’  and  one  for  the  ‘regular 
terms.’  The  system  for  the  ‘x-dependent  terms’  is 
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l2=0 


—  (Creg  +  D)  +  F  . 


(5.133) 

(5.131) 


Since  the  unknowns  $[2^,  and  aPPear  In  the  ‘regular  term’  system, 
it  is  required  that  ‘x-dependent  term’  system  be  solved  first.  To  solve  the  ‘x- 
rlependent  term’  system,  we  first  find  the  angles  by  applying  the  conditions  that 
the  solutions  be  valid  for  all  x  and  all  t.  This  is  accomplished  in  the  same  way  as 
it  was  in  the  0(e)  case — the  x  and  t  derivatives  of  Eqs.  (5.131)  and  (5.132)  are 
taken.  The  resulting  sets  of  equations  are  then  forced  to  be  linearly  dependent 
thereby  making  their  solutions  identical.  To  make  the  x  and  t  derivatives  of 
Eqs.  (5.131)  and  (5.132)  linearly  dependent,  the  following  two  conditions  must 


be  true: 


sing(l)q 

C0 


sin  9" 
co 


(5.135) 


Accordingly,  we  see  that 


sin  1 


(2)« 


sin  6" 
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(5.13(5) 


0$,  =  0$  =  180°  -  0'w 


(5.137) 
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and 

=  0(T  •  (5-138) 

Thus,  the  signals  that  match  the  particular  solution  at  the  interface  propagate  in 
the  same  direction  as  the  0(e)  signal.  That  is  to  say,  these  0(e2)  homogeneous 
solutions  obey  Snell’s  law  and  the  law  of  specular  reflection.  Since  the  angles 
0(7)a  and  0(2)a9  are  known,  the  angles  and  $(2)al  are  a^so  known, 

0$a±  =  ~  90°  =  90°  ~  6'nC  >  (5.139) 

figs.  =  fig-  +  900  •  (5-140) 

Equations  (5.131)  and  (5.132)  may  now  be  solved  for  and  ^^o*- 
Use  of  Eqs.  (5.137)— (5.140)  in  rearranged  forms  of  Eqs.  (5.131)  and  (5.132)  yields 
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(5.142) 

where  we  have  made  use  of  the  following: 
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(5.144) 

(5.145) 

(5.146) 

(5.147) 


Equations  (5.141)  and  (5.142)  may  be  integrated  once  with  respect  to  time  noting 
that  the  integration  constant  must  be  zero  in  order  to  satisfy  quiet  conditions, 
that  is,  in  the  absence  of  an  incident  signal,  a  reflected  or  transmitted  signal  does 
not  exist.  Since  both  Ax  and  Cx  are  given  in  terms  of  the  time  derivative  of 
the  incident  signal  [see  Eqs.  (5.106)  and  (5.109)],  integration  is  possible  in  closed 
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form.  That  this  integration  is  possible  is  important  because  the  expressions 
and  show  up  explicitly  in  the  0(e2)  form  of  the  particle  velocity.  To  obtain 

expressions  for  and  awaY  from  the  interface,  we  replace  the  current 

independent  variable  rg  with  t”®  for  the  0(e2)  reflected  field  and  with  for 
the  0(e2)  transmitted  field.  Thus,  the  expressions  for  and  are 

**«'!?)  -  -2=?=  ((1  - «>  /  ‘M’S?)  *  +  ^ 1 < ' 1  - *>  /  dl) 

(5.148) 

and 
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x  (/ C,(rHr)dt  -  /  AArUndtJ  .(5.1-19) 

The  system  composed  of  ‘regular  terms’  is  now  solved  in  an  identical 
fashion:  First,  the  x  and  t  derivatives  of  the  pair  of  equations  are  taken,  and  the 
two  conditions  that  force  the  x  and  t  derivatives  of  the  equations  to  be  linearly 
dependent  are  obtained.  The  two  conditions  are 


sin  05$, 


(2)6  _ 


sin  0^  _  sin  0inc 


(5.150) 


-v—r-  (s-18l) 

Thus,  the  angles  at  which  the  0(e2)  homogeneous  solutions  leave  the  interface 
are  the  same  as  their  0(e)  counterparts, 

05$,  =  05$  =  180°  -  0inc  ,  (5.152) 


0(2)6  =  =  180°  -  0inc  >  (5.i52) 

^r  =  0(T  •  (5. 1.53) 

That  is  to  say,  these  0(c2)  homogenous  solutions  also  obey  Snell’s  law  and  the 
law  of  specular  reflection.  Using  Eqs.  (5.152)  and  (5.153)  in  the  ‘regular  term’ 
system  [Eqs.  (5.133)  and  (5.134)]  and  then  rearranging  leads  to  the  following 
expression  for  and  at  z  =  0: 
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where  we  have  again  made  use  of  Eqs.  (5.143)— (5.147).  Integration  of  Eqs.  (5.154) 
and  (5.155)  is,  in  principle,  possible,  but  not  in  closed  form.  This  is  not  a  problem, 
however,  because  and  $(£)?*,  unlike  and  do  not  appear  explictly 

in  any  measurable  quantity,  only  their  derivatives.  To  obtain  expressions  fol¬ 
ly  ^(2)6  anc*  §i  ^(2^6 S  away  from  the  interface,  we  replace  the  independent  variable 
Tg  with  for  the  reflected  field  and  with  r^8  for  the  transmitted  field. 

In  summary,  although  we  have  not  obtained  closed  form  solutions  for 
$[2^  and  $[2)/,,  we  have  obtained  closed  form  solutions  for  the  derivat  ives.  This 
is  significant  because  only  the  derivatives  of  appear  in  measurable  quantities 
such  as  the  pressure  or  particle  velocity.  The  derivatives  of  and  $’!,)/,  are 
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where  n^±  and  n|j“8  are  given  by 

n$x  =  i  cos  einc  +  k  sin  einc  ,  (5.  ICO) 

n(i“x  =  » cos  «8T  ~  k  sin  °(\?S  •  (5. id ) 

All  the  components  of  Eqs.  (5.156)— (5.159)  are  known:  and  are  g>von 

by  Eqs.  (5.148)  and  (5.149),  and  their  derivatives,  <t>^a  and  §-t  ^(2)T>  are  read¬ 
ily  obtained.  The  functions  Jj  and  J7  $(£)£*  are  given  by  Eqs.  (5.154)  and 
(5.155).  It  is  also  noted  that  the  0(e2)  homogeneous  solutions  propagate  in  di¬ 
rections  predicted  by  Snell’s  law  and  the  law  of  specular  reflection. 


Matching  the  0(e2)  source  condition 

The  0(c2)  boundary  condition  at  the  source  is  now  matched  to  the  0( f2) 
solution  evaluated  at  the  source,  r  =  r0.  As  we  noted  when  matching  the  0(c) 
boundary  condition  at  the  source,  the  reflected  signal  is  not  involved  in  matching 
the  boundary  condition  at  the  source.  The  0(e2)  boundary  condition  at  the 
source  is  given  in  Eq.  (5.40).  The  relation  between  U(2)  and  the  modified  velocity 
potential  4>(2)  is  given  in  Eq.  (5.30).  Use  of  Eq.  (5.40)  and  the  general  solution 
for  Eq.  (5.80),  in  Eq.  (5.30)  leads  to 
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(5.102) 

The  first  term  on  the  right-hand  side  of  Eq.  (5.162)  forces  the  0(t2)  particular 
solution  to  vanish  at  the  source.  The  physical  interpretation  is  that  no  self- 
action  of  the  incident  wave  has  occurred  when  the  wave  is  still  at  the  source. 
The  second  term  on  the  right-hand  side  accounts  for  the  nonlinear  relationship 
between  the  modified  velocity  potential  and  the  particle  velocity.  As  we  see  below, 
this  term  disappears  if  the  equation  is  re-expressed  in  terms  of  a  measurable 
quantity  such  as  the  pressure  or  particle  velocity.  The  third  term  on  tin'  right- 
hand  side  accounts  for  the  0(e2)  signal  at  the  source.  Physically,  an  0(c2)  signal 
at  the  source  represents  a  local  effect  at  the  source  such  as  the  finite  displacement 
of  the  source.  This  local  efFect  must  not  be  overlooked  because  it  is  of  the  same 
order  as  the  local  effects  at  the  interface — the  finite  displacement  of  the  interface' 
and  the  variation  of  the  normal  to  the  interface.  Evaluating  V4>J2);i  (neglecting 
reflections)  and  rearranging  leads  to  the  following: 
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A  closed  form  integral  of  Eq.  (5.163)  was  not  found.  If  the  source  condi¬ 
tion  is  harmonic,  however,  the  integral  may  be  obtained.  The  form  of  Eq.  (5.163) 
is  nevertheless  acceptable  because  all  measurable  quantities  depend  on  deriva¬ 
tives  of  $  and  not  on  $  itself.  For  example,  we  may  use  Eq.  (5.26)  to  rearrange 
Eq.  (5.163): 
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lAe  have  thus  found  the  expression  for  the  0(e2)  incident  pressure  in  terms  of  the 
0(c)  incident  pressure  [related  to  d&'fij/dt  by  Eq.  (5.25)]  and  the  0(e2)  pressure 
at  the  source  (related  to  dS^/dt).  Other  such  manipulations  are  possible. 


Summary  of  0(c2)  Results 

Because  closed  form  solutions  have  not  been  obtained  for  ^[^6  and  $(2)6S> 
we  have  no  closed  form  solutions  for  and  $|2).  Solutions  for  any  0(c2)  measur¬ 
able  quantity  may,  however,  be  obtained  because  measurable  quantities  depend 
only  on  the  derivatives  of  $j2)  and  The  derivatives  of  the  functions  $j2)p 

and  may  be  obtained  by  differentiating  Eqs.  (5.92)  and  (5.86),  respectively. 
The  derivatives  of  $j2JA  and  are  given  in  Eqs.  (5. 1 56)-(5 . 1 59) . 

The  major  results  of  this  section  are  as  follows:  (1)  Closed  form  ex¬ 
pressions  have  been  obtained  for  the  0(e2)  particular  solutions  for  fluids  I  and  11. 
Eqs.  (5.92)  and  (5.86),  respectively.  (2)  The  propagation  directions  of  the  growth 
terms  in  the  0(e2)  particular  solutions  for  fluids  I  and  II  are  given  by  Snell’s  law 
and  the  law  of  specular  reflection.  (3)  Closed  form  expressions  have  been  ob¬ 
tained  for  the  derivatives  of  the  0(e2)  homogeneous  solutions  for  fluids  I  and  11, 
Eqs.  (5.156)  -(5.159),  respectively.  (4)  The  propagation  direction  of  the  0(c2) 
homogeneous  solutions  are  given  by  Snell’s  law  and  the  law  of  specular  reflection. 
(5)  Expressions  have  been  obtained  for  the  quantities  A  through  G,  which  repre¬ 
sent.  the  motion  of  the  interface,  variation  of  the  normal  to  the  interface,  and  the 
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nonlinearity  of  the  relations  between  the  acoustic  variables,  Eqs.  (5.105) — (5.1 13), 
respectively.  All  of  the  foregoing  have  been  obtained  assuming  an  arbitrary  source 
function. 

5-5  Complex  Transmission  Angle 

In  this  section  we  briefly  point  out  what  steps  should  be  taken  if  the 
incident  angle  is  greater  than  the  critical  angle.  Above  critical  incidence,  the 
transmission  angle  and  the  retarded  time  Ttran*  are  complex.  Moreover,  the  defi¬ 
nitions  of  the  reflection  and  transmission  coefficients,  Eqs.  (5.57)  and  (5.58),  are, 
strictly  speaking,  no  longer  valid  because  the  reflected  and  transmitted  signals 
are  no  longer  in  phase  with  the  incident  signal.  The  problem  may  be  resolved  by 
using  complex  notation.  We  start  by  assuming  an  incident  harmonic  wave,  say 
=  Re[e_,w‘$)"j(x,  z)],  where  i  =  yf—\  and  $)"j(x,2)  is  a  complex  amplitude. 
All  the  time  derivatives  in  linear  terms  are  then  replaced  by  —  iu>.  To  compute 
the  quadratic  terms  which  contribute  to  the  second  harmonic,  we  must  do  only 
a  little  more  work.  For  example,  to  compute  the  time  derivative  of  (3^>(i )jdt)i 
which  appears  in  the  wave  equation  Eq.  (5.2),  we  note  first  that 


where  4*^  is  the  complex  conjugate  of  $(i).  The  time  derivative  of  (<94>(i)/<9<)2 
is,  accordingly,  given  by 
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Thus,  everything  we  have  done  may  be  applied  to  incident  angles  above  the 
critical  angle  provided  that  (1)  the  time  derivatives  c  .r  terms  are  replaced 
by  —iuj,  (2)  the  time  derivatives  of  quadratic  tern.^  -  replaced  by  —2 iu,  and 
(3)  the  amplitude  of  the  quadratic  terms  is  multiplied  by 

The  more  general  case  of  an  incident  pulse  may  also  be  treated  using 
complex  notation.  The  incident  pulse  is  considered  as  a  superposition  of  harmonic 
waves.  Care  must  exercised  in  this  case,  however,  because  each  harmonic  interacts 
not  only  with  itself  (as  above),  but  also  with  every  other  harmonic.  Thus,  it  turns 
out  that  the  0( c2)  solution  contains  convolution  integrals  of  the  Fourier  transform 
of  the  pulse  with  itself. 
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5-6  Summary 

In  this  chapter,  we  have  analyzed  the  reflection  and  refraction  of  finite- 
amplitude  plane  waves  that  are  obliquely  incident  on  an  initially  plane  fluid-fluid 
interface.  The  procedure  used  is  second-order  perturbation  analysis;  the  source 
boundary  condition  was  arbitrary.  Specific  terms  that  account  for  the  displace¬ 
ment  of  the  interface,  the  variation  of  the  normal  to  the  interface,  and  the  0(e2) 
source  condition  are  identified.  It  is  found  that  the  0(e2)  reflected  and  transmit¬ 
ted  waves  propagate  in  the  same  direction  as  their  0(e)  counterparts.  In  other 
words,  to  0(e2),  no  deviation  from  Snell’s  law  or  the  law  of  specular  reflection  is 
seen.  Expressions  for  the  amplitude  of  the  0(e2)  reflected  and  transmitted  waves 
have  been  obtained. 


CHAPTER  6 


MODIFIED  FORMS  OF  SNELL’S  LAW  BASED  ON  SIMPLE 

WAVE  FLOW 


6—1  Introduction 

Developed  in  this  chapter  are  two  different  ‘modified  forms’  of  Snell’s 
law — forms  which  appear  to  indicate  that  refraction  has  a  slight  amplitude  depen¬ 
dence.  The  ‘modified  forms’  are  developed  by  means  other  than  perturbations. 
Specifically,  the  first  ‘modified  form’  is  obtained  by  matching  the  trace  velocities 
of  the  incident  and  the  transmitted  signals  at  a  fixed  interface.  This  ‘modified 
form’  was  reported  by  Cotaras  and  Blackstock  (1987).  The  second  ‘modified 
form’  is  obtained  by  examining  the  variation  of  the  pressure  along  the  moving 
interface  and  was  suggested  by  Naze  Tj0tta  and  Tjptta  (1988).  While  different 
from  each  other,  both  forms  appear  to  be  correct  to  second-order.  However,  in 
deriving  the  ‘modified  forms,’  it  is  assumed  that  simple  wave  flow  exists  in  both 
fluids.  We  are  now,  however,  in  a  position  to  quantify  the  simple  wave  flow 
approximation.  Using  the  results  of  the  previous  chapter,  we  develop  the  0(e) 
and  0(e2)  conditions  for  simple  wave  flow  to  exist  in  both  fluids  simultaneously. 
These  conditions  are  then  imposed  on  ‘modified  forms.’ 

This  chapter  is  divided  into  two  sections  as  follows.  In  the  first  section, 
conditions  for  simple  wave  flow  in  fluids  I  and  II  are  developed.  The  0(e)  and 
0(e2)  conditions  for  simple  wave  flow  in  fluid  I  are  obtained  by  determining  the 
conditions  for  no  0(e)  and  0(e2)  reflections.  We  then  examine  the  conditions 
required  for  simple  wave  flow  to  exist  simultaneously  in  fluid  II.  Developed  in 
the  econd  section  of  this  chapter  are  the  two  different  ‘modified  forms’  of  Snell's 
law.  It  turns  out  that,  to  0(e2),  the  two  methods  result  in  equivalent  expressions: 
ordinary  Snell’s  law  and  one  of  the  conditions  for  simple  wave  flow,  which  is 
assumed  in  their  derivation.  The  ‘modified  forms’  are,  therefore,  equivalent  to 
ordinary  Snell’s  law. 
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6-2  Conditions  for  0(e)  and  0(e2)  Simple  Wave  Flow 


Simple  Wave  Flow  in  Fluid  I 


In  this  section,  the  results  of  the  previous  chapter  are  used  to  determine 
the  0(e)  and  0(e2)  conditions  for  simple  wave  flow  in  fluid  I.  This  is  done  by 
determining  the  0(e)  and  0(e2)  conditions  for  no  reflections.  The  condition  for 
no  0(c)  reflection  was  noted  in  the  previous  chapter  [see  Eq.  (5.70)]  and  was 
given  the  name  intromission.  Thus,  all  that  remains  to  be  done  is  to  determine 
the  0(e2)  condition  for  no  reflection.  Our  procedure  is  straightforward:  First, 
the  0(e)  condition  is  imposed  on  the  0(e2)  general  solution  for  fluid  I  from  the 
previous  chapter.  The  0(e2)  condition  for  simple  wave  flow  in  fluid  I  is  then 
determined  by  forcing  the  0(e2)  reflection  to  vanish. 

To  impose  the  0(e)  condition,  first  recall  that  the  0(e2)  general  solution 
for  fluid  I  is  divided  into  two  parts,  a  particular  solution  and  a  homogeneous 
solution.  At  intromission  the  0(e)  reflection  is  by  definition  zero,  and,  accord¬ 
ingly,  the  terms  in  the  0(e2)  particular  solution  for  fluid  I  that  correspond  to 
the  interaction  of  the  incident  and  reflected  waves  and  to  the  self-action  of  the 
reflected  wave  vanish.  The  0(e2)  condition  for  simple  wave  flow  in  fluid  I  is  thus 
determined  solely  by  forcing  the  0(e2)  homogeneous  solution  for  fluid  I  to  vanish. 

Imposing  the  0(e)  condition  (R  —  0)  on  the  terms  that  make  up  the 
0(e2)  homogeneous  solution  for  fluid  I  [the  terms  A  through  G  which  are  given 
in  Eqs.  (5.105)— (5.1 13)],  we  obtain  the  following: 
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Note  that  G,  the  term  that  accounts  for  the  variation  of  the  normal,  is  zero  for 
normal  incidence.  (The  condition  for  intromission  at  normal  incidence  is  m  =  n.) 
For  all  other  angles,  however,  m  must  be  unity  for  G  to  be  zero.  However,  if  m 
is  unity,  then  the  0(e)  condition  [see  Eq.  (5.70)]  in  combination  with  Snell’s  law 
indicates  that  n  is  unity  also,  that  is,  n  =  m  —  1.  If  both  m  and  n  are  unity,  then 
E,  F,  and  G,  which  together  account  for  all  interface  motion,  are  zero.  Note  that 
even  with  n  and  m  equal  to  unity,  the  term  /4reg  is  not  zero  unless  /?’  =  0n. 

We  now  determine  the  form  of  the  0(t2)  homogeneous  solution  given 
that  R  —  0.  Using  Eqs.  (6.2)  and  (6.5)  in  Eqs.  (5.148)  and  (5.149),  we  see  that 
for  the  special  case  of  R  —  0 
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Note  that  the  part  of  the  homogeneous  solution  for  fluid  II  that  has  am¬ 

plitude  variation  along  its  wavefront,  is  zero  in  the  case  of  normal  incidence.  For 
oblique  incidence,  $1™“  is  zero  only  if  /freg  is  zero.  Use  of  Eqs.  (6.6),  (6.10),  and 
(6.11)  in  Eqs.  (5.154)  and  (5.155)  leads  to 
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Equations  (6.10)  and  (6.12)  are  the  two  terms  that  make  up  the  0(e2)  homoge¬ 
neous  solution  for  fluid  I.  Since  Eq.  (6.10)  shows  that  is  already  zero,  the 
0(e2)  conditions  for  simple  wave  flow  in  fluid  I  may  be  obtained  by  forcing  <J>^6 
to  be  zero. 

An  examination  of  Eqs.  (6.1),  (6.7),  (6.8),  and  (6.9)  indicates  that  for 
<$(2?6  to  be  zero,  the  coefficients  of  (d$mc/dt)2  and  <923>jjj/<9i2  must  be  zero. 
Setting  the  coefficient  of  d^Q'j^/dt2  equal  to  zero  yields  the  0(e)  condition 
for  intromission,  which  was  previously  assumed.  No  new  conditions  are  therefore 
required  to  force  the  coefficient  of  <92$[")/<9t2  to  be  zero.  On  the  other  hand, 
forcing  the  coefficient  of  (d<$>,nc  /  dt)2  to  be  zero  leads  to  the  following  condition: 
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Note  that  at  normal  incidence,  Eq.  (6.14)  reduces  to  (with  m  =  n) 


(6.15) 


Equation  (6.15)  indicates  that  for  no  0(e2)  reflection  at  normal  incidence,  the 
jump  in  the  quantity  (3/cq  must  be  zero.  For  oblique  incidence  the  condition  is, 
however,  more  complicated.  The  angle  dependence  in  Eq  (6.14)  may  be  removed 
using  the  0(e)  condition  [see  Eq.  (5.71)],  and  we  thereby  obtain  the  following 
form  of  the  condition  for  no  0(e2)  reflection: 
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Thus,  for  oblique  incidence,  simple  wave  flow  exists  in  fluid  1  if  the 
properties  of  the  fluid  pair  are  such  that  both  the  0(e)  and  0(e 2)  conditions  for 
no  reflections  [Eqs.  (5.71)  and  (6.16),  respectively]  are  met  simultaneously.  For 
normal  incidence,  the  number  of  conditions  for  simple  wave  flow  is  the  same — 
two,  but  the  conditions  themselves  are  simplier  in  form.  The  0(()  condition  for 
no  reflection  at  normal  incidence  is  obtained  from  Eq.  (5.70), 


SfoCoJ  =  0 


(6.17) 


The  0( f2)  condition  for  no  reflection  at  normal  incidence  is  Eq.  (6.15). 
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Simple  Wave  Flow  in  Fluids  I  and  II 

We  now  examine  the  conditions  for  simple  wave  flow  to  exist  simulta¬ 
neously  in  fluids  I  and  II.  We  start  by  considering  what  restriction  the  simple 
wave  flow  assumption  places  on  the  relationship  between  the  temporal  and  spa¬ 
tial  derivatives  of  an  acoustic  field  variable.  The  restriction  is  central  to  many 
relations  that  are  derived  assuming  simple  wave  flow. 

As  an  example,  we  develop  the  nonlinear  impedance  relation  for  the 
special  case  of  outgoing  simple  waves.  Our  starting  point  is  the  second-order 
form  of  the  momentum  equation,  Eq.  (2.88).  If  simple  wave  flow  in  a  lossless 
fluid  is  assumed,  Eq.  (2.88)  simplifies  to 

P0¥  +  Vp/  =  °  '  (6-18) 

If  the  acoust  ic  field  consists  of  simple  outgoing  waves,  the  time  and  space  deriva¬ 
tives  of  the  pressure  are  related  by 
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where  n  is  the  direction  of  propagation.  [Recall  that  the  propagation  speed  of  a 
finite-amplitude  wave  is  Co  +  /3(n  *u).]  Substituting  the  foregoing  in  Eq.  (6.18), 
integrating,  and  noting  that  the  integration  constant  must  be  zero  in  order  to 
satisfy  quiet  conditions  leads  to 
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Equation  (6.20)  is  the  second-order  form  of  the  impedance  relation  and  is  valid 
for  simple  outgoing  waves  of  finite  amplitude. 

Equation  (6.19)  indicates  that,  under  the  simple  wave  flow  assumption, 
the  temporal  and  spatial  derivatives  are  simply  related;  moreover,  the  gradient 
of  the  pressure  is  parallel  to  the  direction  of  propagation.  The  particle  motion  is, 
accordingly,  restricted  to  being  parallel  to  the  direction  of  propagation. 

To  obtain  a  condition  for  simple  wave  flow  in  fluid  II,  we  examine  the 
gradient  of  the  0(c)  and  0(c 2)  pressures  in  fluid  II  and  demand  that  they  be 
parallel  to  the  direction  of  propagation.  The  gradient  of  the  0(e)  pressure  in 
fluid  II  is  obtained  by  taking  the  gradient  of  Eq.  (5.27).  Noting  Eq.  (5.50),  we 
find  that,  to  0(f),  the  gradient  of  the  pressure  is  parallel  to  the  direction  of 
propagation.  Tims  no  additional  conditions  are  required  to  assume  simple  wave 
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flow  exists  in  fluid  II  to  0(e).  The  situation  is  more  complicated  at  0(e2).  The 
0(e2)  pressure  in  fluid  II  is  given  by  Eq.  (5.28).  The  second  term  on  the  right- 
hand  side  of  Eq.  (5.28)  is  a  quadratic  form  of  the  0(e)  solution,  and  the  gradient 
of  this  term  is,  accordingly,  parallel  to  the  direction  of  propagation.  The  first  term 
on  the  right-hand  side  of  Eq.  (5.28)  may  be  divided  into  two  parts:  a  particular 
solution  and  a  homogeneous  solution.  From  Eq.  (5.86)  we  see  that  the  fluid  II 
particular  solution  is  also  a  quadratic  form  of  the  0(e)  solution,  the  gradient 
of  which  is  parallel  to  the  direction  of  propagation.  The  fluid  II  homogeneous 
solution  is  given  by  Eq.  (5.118).  By  inspection  we  see  that  for  the  gradient  of 
^(2)h  be  Paradel  to  the  direction  of  propagation,  must  be  zero.  We  noted 

earlier  that,  for  the  special  case  of  R  =  0,  *s  zero  f°r  normal  incidence,  but 

for  oblique  incidence,  is  zero  only  if  y4reg  is  zero.  Furthermore,  y4reg  is  zero 

only  if 


V  P1  m 


(6.21) 


Since  Eq.  (6.21)  is  the  same  as  the  normal  incidence  condition  for  0(e 2) 
simple  wave  flow  in  fluid  I  [Eq.  (6.15)],  no  additional  conditions  are  required 
for  0(e2)  simple  wave  flow  to  exist  simultaneously  in  fluids  I  and  II  at  normal 
incidence.  For  oblique  incidence,  however,  Eq.  (6.21)  represents  an  additional 
condition.  Thus,  for  0(e2)  simple  wave  flow  to  exist  simultaneously  in  both  fluids 
I  and  II  at  oblique  incidence,  the  conditions  expressed  in  Eqs.  (5.71),  (6.1  1).  and 
(6.21)  must  be  met  simultaneously.  The  only  way  all  three  conditions  may  be 
satisfied  simultaneously  is  if  the  static  densities,  the  small-signal  sound  speeds, 
and  the  coefficients  of  nonlinearity  are  matched, 
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,6.23) 

(6.24) 


6-3  Modified  Forms  of  Snell’s  Law 


In  this  section,  two  ‘modified  forms’  of  Snell’s  law  are  developed.  The 
first  is  obtained  by  matching  the  trace  velocities  of  the  incident  and  transmitted 
signals  at  the  interface,  while  assuming  that  the  interface  is  stationary.  The 
second  ‘modified  form’  of  Snell’s  law  is  obtained  by  examining  the  variation  of 
the  pressure  along  the  moving  interface.  The  0(e)  and  0(e2)  results  of  the  two 
methods  are  found  to  be  equivalent.  Moreover,  when  the  conditions  for  simple 
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wave  flow,  which  is  assumed  in  their  derivation,  are  imposed  on  the  ‘modified 
forms’,  they  reduce  to  ordinary  Snell’s  law. 


Trace  velocity  matching1 

This  section  progresses  as  follows:  First,  ordinary  Snell’s  law  is  derived 
by  matching  the  trace  velocities  of  small-signal  waves.  Then,  a  ‘modified  form’  of 
Snell’s  law  is  developed  by  matching  the  trace  velocities  of  finite  amplitude  waves 
at  the  interface.  Simple  wave  flow  is  assumed  on  the  fluid  I  side,  and  the  interface 
is  assumed  stationary.  The  0(c)  and  0(e2)  approximations  of  the  ‘modified  form’ 
are  then  developed. 

The  method  of  matching  the  trace  velocities  of  the  incident  and  trans¬ 
mitted  signals  is  first  demonstrated  by  deriving  Snell’s  law  for  small-signal  waves. 
Shown  in  Fig.  6.1  is  a  plane  wave  obliquely  incident  on  a  plane  interface.  First, 
consider  fluid  I.  We  follow  the  progress  of  a  single  typical  wavelet  of  the  incident 
wave,  in  this  case,  a  peak.  At  time  t  =  0,  the  wavefront  for  this  wavelet  is  AA! .  At 
a  time  At  later,  the  wavefront  has  moved  a  distance  c^At  to  position  B  B' .  The 
increment  the  wavelet  has  moved  along  the  interface,  that  is,  the  trace  distance 
Axab  is 


^ XAB  ~ 


•  (6-25) 

sin  0,nc  v  ’ 

Note  that  Cy/  sin  0mc  is  the  trace  velocity  of  the  wavelet  along  the  interface.  Next 
consider  fluid  11.  The  wavelet  of  the  transmitted  wave  at  time  t  =  0  is  represented 
by  the  wavefront  AA”  and,  at  a  time  At  later,  by  BB".  The  distance  along  the 
ray  is  c]0[At  and  the  trace  distance  is  given  by 


Ax  -  4‘A < 

4B  5;n0trans 


(6.2C) 


Thus,  the  trace  velocity  is  Cy/sin#“““s.  Combining  the  two  equations  leads  to 
Snell’s  law. 


J 


cn 

-o 


(6.27) 


sin  0,nc  sin  0lrans 

It  appears  that,  Snell's  law  is  simply  a  kinematic  condition  t  hat  requires  the  trace 
velocities  to  be  the  same  for  both  the  incident  and  transmitted  waves.  (Note 


1  The  trace  velocity  matching  technique  and  the  resulting  ‘modified  form’  of  Snell’s  law  wore 
presented  at  the  lMtli  meeting  of  the  Acoustical  Society  of  America  (Cotaras  and  Hlackst.ock 
HIST). 
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Figure  G.l  Small-signal  plane  wave  obliquely  incident  on  a  plane  fluid-fluid 
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that  the  trace  velocity  matching  method  may  also  be  used  to  obtain  the  law  of 
specular  reflection.) 

The  trace  velocity  matching  technique  is  now  used  for  finite-amplitude 
waves.  Recall  that  the  propagation  velocity  of  a  simple  outgoing  wave  of  finite 
amplitude  is  +  [3l(n  *u) — not  just  the  small-signal  sound  speed  Cq.  Consider 
Fig.  6.2.  By  invoking  the  same  geometrical  argument  used  for  small-signal  waves, 
we  obtain  the  following  ‘modified  form’  of  Snell’s  law  for  finite  amplitude  waves, 


sin  n'mc 


(C.28) 


Note,  however,  that  simple  wave  motion  was  assumed  in  fluid  I  and  that  the 
motion  of  the  interface  was  neglected. 

The  0(e)  and  0(e2)  forms  of  the  ‘modified  form’  of  Snell’s  law  are  now 
developed.  We  first  expand  Eq.  (6.28)  using  the  techniques  developed  in  Chapter 
5,  specifically,  Eqs.  (5.11 )— (5.1 3).  Using  Eqs.  (5.12)  and  (5.13)  in  Eq.  (6.28),  then 
binomially  expanding  the  denominators  and  retaining  only  the  leading  and  next 
higher-order  terms  leads  to 


sin  flmc  ( 

4  V 


c0 


trans  trans \ 
*U(1)  ) 


(6.29) 


Equating  the  leading-order  terms  in  Eq.  (6.29)  yields  the  following  0(e)  form  of 
the  ‘modified  form’  of  Snell’s  law: 

sin  0>nc  sin0trans 

This  is  merely  the  ordinary  form  of  Snell’s  law.  It  was  noted  in  the  previous 
chapter  that  Snell’s  law  is  one  of  the  conditions  required  for  the  0(e)  solutions 
for  the  wave  equation  to  be  valid  for  all  x  and  all  t.  Equating  the  next  higher- 
order  terms  of  Eq.  (6.29)  and  using  Snell’s  law  leads  to  the  0(e2)  form  of  the 
'modified  form’  of  Snell’s  law, 


,  Pl  ml 


(6.30) 


In  obtaining  Eq.  (6.30),  wc  have  made  use  of  the  following  0(e)  result  ,  which 
was  in  turn  obtained  using  the  0(e)  pressure  balance.  Eq.  (5.35),  and  the  O(e) 
impedance  relation,  p  =  /)0c0(n  *u): 

,,inc  ,.i»c  ^  —  trails  ..trans  t  c  ■  >  i  \ 

11  ’U(l)  —  11  *U(1)  (f»->  1  > 
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Figure  6.2  Finite-amplitude  plane  wave  obliquely  incident  on  a  plane  fluid-fluid 
INTERFACE 
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It  was  noted  in  the  previous  section,  however,  that  Eq.  (6.30)  is  a  condition  for 
simple  wave  flow  to  exist  in  fluid  II  to  0(c2);  see  Eq.  (6.21).  Thus,  to  0(c2),  the 
‘modified  form’  of  Snell’s  law  and  the  ordinary  form  are  equivalent.  That  is  to 
say,  the  ‘modified  form’  of  Snell’s  law  is  a  trivial  non-extension  of  Snell’s  law. 
The  problem  with  the  derivation  of  the  ‘modified  form’  is  that  the  conditions 
imposed  by  the  simple  wave  flow  assumption  were  overlooked. 

In  summary,  it  has  been  shown  that,  to  0(c2),  the  ‘modified  form’  of 
Snell’s  law  that  is  obtained  by  matching  the  trace  velocities  reduces  to  ordinary 
Snell’s  law  when  the  conditions  of  the  simple  wave  flow  are  imposed  on  it. 

Variation  of  pressure  along  the  interface2 

In  this  section,  another  ‘modified  form’  of  Snell’s  law  is  derived — this  one 
by  examining  the  variation  of  the  pressure  along  the  interface.  A  Taylor  series 
expansion  is  then  used  to  obtain  0(c)  and  0(e2)  approximation  of  this  ‘modified 
form’  along  z  =  0.  Simple  wave  flow  is,  however,  again  assumed  on  the  fluid  I 
side,  and  it  turns  out  that  this  ‘modified  form’  of  Snell’s  law  yields  the  same  0(c) 
and  0(c2)  expressions  that  were  obtained  using  the  trace  velocity  method. 

We  first  show  that  the  jump  in  the  variation  of  the  pressure  along  the 
interface  is  zero.  A  simple  way  to  show  this  is  to  describe  the  interface  in  a 
parametric  form  using  only  a  length  coordinate  i  and  the  time  t.  The  pressure 
jump  at  a  point  £  is  then  written 

p\£,  t)  -  pn(£,  t)  =  0  •  (6.32) 

Similarly,  at  a  point  i  -f  Af,  we  may  write 

p\£  +  AM)  -pn{£  +  Af,f)  =  0  .  (6.33) 


Combining  the  two  previous  expressions  and  dividing  by  A£  yields  the  desired 
result, 


dp1  dp 11 

Hi  =  ~df 


(6.34) 


The  variation  of  the  pressure  along  the  interface  may  also  be  expressed  in  terms 
of  the  tangent  to  the  interface  T  and  the  gradient  of  the  pressure, 


X -Vp1  =  T -Vp11 


(6.35) 


The  second  ‘modified  form’  of  Snell’s  law  may  be  obtained  by  applying 
Eq.  (6.35)  along  the  moving  interface.  If  simple  wave  flow  is  assumed  on  the 

2 This  method  was  suggested  6y  Naze  Tjotta  and  Tjfltta  (1988). 
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fluid  I  side,  Eq.  (6.19)  may  be  used  to  obtain  an  expression  for  the  gradient. 
Thus,  along  the  interface,  Eq.  (6.35)  may  be  written  as 

sin  $inc  dpinc  sin0tr““  dp*™  /e 

4  +  fll(rimc  *u'nc)  dt  ~  41  +  /9II(ntr“ls  •utrans)  dt  ’  (  •  ) 


where  the  angles  0,nc  and  0trans  are  defined  relative  to  the  instantaneous  normal. 
Except  for  the  dp/di  terms,  Eq.  (6.36)  is  the  same  as  our  previous  ‘modified 
form’  of  Snell’s  law,  Eq.  (6.28),  which  was  obtained  neglecting  interface  motion 
altogether. 

If  this  ‘modified  form’  of  Snell’s  law  is  to  be  compared  with  the  previ¬ 
ous  one,  the  angles  must  be  referred  to  the  +z-axis  and  the  displacement  of  the 
interface  must  be  accounted  for.  Moreover,  the  entire  expression  must  be  ap¬ 
proximated  to  0(t2).  This  may  be  accomplished  as  follows:  First,  an  expression 
for  the  tangent  is  obtained  using  the  previously  derived  0(e)  expression  for  the 
normal  to  the  interface,  Eq.  (4.19).  Next,  an  0(e2)  expression  for  the  gradient 
of  the  pressure  is  obtained  by  expanding  Eq.  (6.19)  and  using  a  Taylor  series 
expansion  about  z  =  0.  This,  of  course,  assumes  simple  wave  flow  on  the  fluid  I 
side.  The  two  expressions  are  then  combined  to  yield  the  0(e)  and  0(e2)  forms 
of  Eq.  (6.35)  about  z  —  0. 

An  expression  for  the  tangent  may  be  obtained  from  our  expression  for 
the  normal  to  the  interface  N,  Eq.  (4.19).  Since  T  must  be  perpendicular  to  N, 


T  is  given  by 


T  =  ±i  ±  ke  j 


dw'ff* 


(6.37) 


The  +  sign  is  chosen. 

If  simple  wave  flow  is  assumed  on  the  fluid  I  side,  Eq.  (6.19)  may  be  used 
to  obtain  an  expression  for  the  gradient.  Expanding  the  denominator  binominally 
and  accounting  for  the  motion  of  the  interface  using  a  Taylor  series  expansion 
results  in  the  following  expression  for  the  gradient  of  the  pressure  along  z  =  0: 


The  expression  for  the  gradient  of  the  pressure  in  fluid  II  is  similar. 


92 


Using  Eqs.  (6.37)  and  (6.38)  in  Eq.  (6.35)  and  equating  equal  powers  of 
e  yields  the  following  0(e)  and  0(e2)  results: 

sin0inc  <V("j  sin  0trani  8p\\™ 

4  ,=0”  41  at  ,=0  ’  (6‘39) 


sin  9'nc  dp\ sin  flinc  d2p{^ 


P(2)  _  sin  f 

dt  4  dtdz  J 


w  ir8  dt 


cos0incdp‘{")  fOW(\j-  ji  sir 

4  dt  J  ~dT~dt  +  ~ 

sin  sin  fltraiu  ^2P(i) 


du>  '  sin  0inc  /31  dpm  ■ 

dx  4  4  ar 


<92  <9z 


...trans  j/ 

10(1)  dt 


cos  fltrans  ^P(if 8  [(MT,  sin  fltri 
cl,1  dt  J  dx  *  Cq1 


rans  oil  dp}™™ 

^  rl*l  /..trans  _.trans\ 

"?“ar(n  •“'»  > 


(6.40) 

Fortunately,  the  0(e)  and  0(e2)  results  given  in  Eqs.  (6.39)  and  (6.40) 
may  be  simplified.  Use  of  the  time  derivative  of  the  0(e)  pressure  balance  given 
in  Eq.  (5.35)  reduces  Eq.  (6.40)  to  Snell’s  law, 

sin0inc  sin  0trana 

4  ~  41  ' 

Equation  (6.40)  may  be  simplified  by  noting  that,  for  simple  wave  flow,  the  x- 
derivative  of  the  0(e2)  pressure  balance,  Eq.  (5.36),  is 

_  5“ eJl  dJk  _  tw*~  M 

4  at  |t=0  4  dtdz  J  (1)  i=0 

co se^~dp\T)  [<HT  ,| 

ck  dt  J  dx  _ 


...  trails  /j 
Z0(i)  ^ 


sin  0tr 


Antram 

°P(2) 


sinfltranS<92P(5r  r  tra 

4  at  a*  JW{1] 

f  dw\T  J 


(6.41) 


Use  of  Eqs.  (6.41),  (6.31),  and  Snell’s  law  reduces  Eq.  (6.40)  to 
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This  is  the  same  expression  obtained  using  the  trace  velocity  approach,  Eq.  (6.30). 

Thus,  the  ‘modified  form’  of  Snell’s  law  that  is  obtained  by  examining 
the  variation  of  the  pressure  along  the  interface  is,  to  0(e2),  equivalent  to  that 
obtained  by  matching  the  trace  velocities  if  the  results  are  referred  to  z  —  0. 
Moreover,  both  ‘modified  forms’  are,  to  0(t2),  equivalent  to  ordinary  Snell’s  law. 
This  implies  that  neglecting  the  motion  of  the  boundary  in  the  derivation  of 
the  first  ‘modified  form’  is  not  a  problem.  The  problem  is  the  simple  wave  flow 
assumption  that  is  made  when  deriving  both  of  the  modified  forms. 

6-4  Summary 

In  this  chapter,  we  examined  the  conditions  for  0(e)  and  0(e2)  simple 
wave  flow  in  fluid  I  and  fluid  II.  For  oblique  incidence,  the  conditions  are  that  the 
static  densities,  the  small-signal  sound  speeds,  and  the  coefficients  of  nonlinearity 
be  matched.  For  normal  incidence  it  is  only  required  that  p0co  and  poco//^  be 
matched.  Also  examined  in  this  chapter  are  two  ‘modified  forms’  of  Snell’s  law. 
It  was  found  that  the  two  ‘modified  forms’  were,  to  0(e2),  equivalent  and  that 
they  reduce  to  ordinary  Snell’s  law  and  one  of  the  conditions  for  simple  wave 
flow,  which  was  assumed  in  their  derivation.  Thus,  to  0(e2),  they  are  equivalent 
to  ordinary  Snell’s  law. 


CHAPTER  7 


SUMMARY  OF  RESULTS 
AND  PROPOSALS  FOR  FUTURE  WORK 


7-1  Summary 

This  dissertation  is  divided  into  three  parts:  In  the  first  part  the  basic 
equations  for  a  homogeneous,  thermoviscous  fluid  with  a  single  relaxation  mecha¬ 
nism  were  examined.  The  examination  led  to  the  following  results:  (1)  Simplified 
(correct  to  second  order)  forms  of  the  equations  for  a  homogeneous,  thermovis¬ 
cous  fluid  with  a  single  relaxation  mechanism  were  developed.  (2)  The  wave 
equation  for  finite- amplitude  wave  motion  in  such  a  fluid  was  also  developed. 
The  derivation  of  this  wave  equation  from  first  principles  has  not  been  presented 
previously.  The  following  assumptions  were  made  during  the  derivation:  (1)  The 
deviation  from  equilibrium  was  assumed  to  be  small,  and,  accordingly,  linear  re¬ 
lations  between  the  thermodynamic  fluxes  and  forces  were  used.  (2)  The  possibil¬ 
ity  of  a  cross-effect  between  bulk  viscosity  and  relaxation  was  neglected.  (3)  The 
magnitude  of  the  acoustic  signal  was  assumed  to  be  small,  but  large  enough  to 
make  quadratic  nonlinearity  terms  significant.  (4)  The  magnitude  of  transport 
coefficients  was  assumed  to  be  small,  and  (5)  the  amount  of  dispersion  was  also 
assumed  to  be  small. 

In  the  second  part  of  the  dissertation,  the  reflection  and  refraction  of  the 
finite-amplitude  sound  at  an  interface  between  two  lossless  fluids  was  analyzed. 
The  analysis  yielded  the  following  results:  (1)  Expressions  that  are  correct  to 
second  order  were  derived  for  the  kinematic  and  force  balance  boundary  condi¬ 
tions.  (2)  Expressions  for  the  0(c)  and  0(e2)  reflected  and  transmitted  fields 
were  obtained,  and  specific  terms  that  account  for  the  displacement  of  the  inter¬ 
face,  the  variation  of  the  normal  to  the  interface,  and  the  0(e2)  source  condition 
were  identified.  (3)  It  was  noted  that,  to  0(e2),  no  deviation  from  Snell's  law 
or  the  law  of  specular  reflection  is  predicted.  None  of  these  expressions  have 
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been  previously  obtained,  and  all  expressions  were  obtained  for  an  arbitrary  (not 
necessarily  harmonic)  source  condition. 

The  significant  results  of  the  third  part  of  this  dissertation  are  as  follows: 

(1)  The  requirements  for  simple  wave  flow  to  exist  simultaneously  in  fluids  I  and 
II  are  that  the  small-signal  sound  speeds,  the  static  densities,  and  the  coefficients 
of  nonlinearity  be  matched.  (2)  The  two  ‘modified  forms’  of  Snell’s  law  are,  to 
<9(f2),  equivalent.  (3)  When  the  conditions  of  simple  wave  flow  in  fluids  I  and 
II  are  imposed  on  the  ‘modified  forms’  of  Snell’s  law,  they  reduce  to  ordinary 
Snell’s  law. 

7-2  Future  Work 

The  work  presented  in  this  dissertation  is  much  less  extensive  than  the 
work  that  remains  to  be  done.  Two  ideas  for  future  work  on  the  fundamental 
assumptions  in  the  derivation  of  the  basic  equations  are  the  following:  The  first 
is  to  examine  the  possibility  of  nonlinear  thermodynamic  effects.  Such  effects 
occur  because  the  deviation  from  equilibrium  is  not  small  enough  to  be  correctly 
predicted  by  a  linear  model.  For  the  case  of  a  chemical  reaction  as  the  relaxat  ion 
mechanism,  Prigogine  (1961)  points  out  that  the  deviation  from  equilibrium  is 
frequently  large  enough  to  require  the  use  of  nonlinear  relations  between  the 
thermodynamic  fluxes  and  forces.  The  second  idea  for  future  work  is  to  include 
the  effects  of  diffusion.  Diffusion  may  be  important  near  the  interface  of  two 
fluids.  Moreover,  the  effects  of  diffusion  are  known  to  interact  with  the  effects  of 
heat  conduction;  see  Prigogine  (1961). 

Future  work  on  the  problem  of  reflection  and  refraction  of  finite-ampli¬ 
tude  sound  should  be  oriented  towards  reducing  the  idealizations  of  this  work  by 
performing  the  two  following  extensions.  (1)  The  effects  of  losses  in  the  fluids 
should  be  included  in  both  the  wave  equations  and  the  boundary  conditions. 

(2)  The  effects  of  finite  source  size  should  be  included  by  using  results  of  this 
dissertation  as  one  of  the  components  in  the  spatial  Fourier  decomposition  of 
the  source.  The  results  of  this  work  may,  however,  be  readily  used  to  analyze 
two  special  cases  that  are  easily  realized  experimentally,  namely,  reflections  from 
a  rigid  wall  [R  =  1)  and  reflections  from  a  pressure  release  surface  ( R  =  — 1). 
In  the  long  term,  a  numerical  implementation  of  the  results  would  be  useful  for 
computing  the  transmitted  and  reflected  fields  of  a  given  source,  both  near  and 
far  from  the  interface.  In  this  way  many  different  source  receiver  geometries  and 
fluid  pairs  could  be  examined. 


APPENDIX  A 


THERMODYNAMICS  OF  RELAXING  FLUIDS 


A-l  Introduction 

It  is  useful  to  know  the  origins  of  the  fundamental  thermodynamic  rela¬ 
tions  that  apply  to  acoustics.  The  majority  of  the  expressions  that  are  developed 
in  this  appendix  are  well  known  for  nonrelaxing  fluids  (sec,  for  example,  Van- 
Wylen  and  Sonntag  (1976)),  but  not  so  well  known  for  relaxing  fluids.  Good 
references  on  the  thermodynamics  of  relaxing  fluids  and  its  application  to  acous¬ 
tics  do,  however,  exist;  see,  for  example,  de  Groot  and  Mazur  (1962)  and  Bhatia 
( 1967).  The  interested  reader  is  referred  to  these  references  for  more  details.  The 
purpose  of  deriving  fundamental  relations  in  this  appendix  is  merely  complete¬ 
ness. 

The  notation  used  in  this  appendix  is  as  follows:  Frozen  quantities  are 
denoted  by  a  superscript  °°,  equilibrium  values  by  a  superscript  °.  The  reason 
for  this,  which  is  elaborated  in  the  main  text,  is  that  the  frozen  value  is  the 
correct  value  for  infinite  frequency;  the  equilibrium  values  are  correct  for  zero 
frequency.  'File  notation  used  to  denote  the  static  value  of  any  thermodynamic 
variable  is  the  subscript  o-  If  a  fluid  does  not  contain  a  relaxation  mechanism, 
no  difference  between  the  frozen  and  equilibrium  states  exists.  The  notation 
used  to  denote  variables  in  this  situation  is  the  same,  except  the  superscripts 
are  dropped.  For  example,  the  static  values  of  the  specific  heats  at  constant 
pressure  and  constant  volume  for  a  nonrelaxing  fluid  are  denoted  cVq  and  c„0. 
respectively.  All  variables  in  this  appendix  are  dimensional,  and,  consequently, 
no  special  notation  is  required  to  indicate  dimensionality. 

A -2  Relations  from  Gibbs  Equation 

Considered  in  this  section  are  the  relations  that  may  be  obtained  from 
Gibbs  equation.  Use  of  Gibbs  equation  leads  to  differential  equations  for  the 
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pressure,  density,  and  temperature.  Relations  between  derivatives  of  the  temper¬ 
ature  and  pressure  may  then  be  obtained.  We  start  by  stating  one  form  of  Gibbs 
equation, 


de  =  Tds  - ) — -  dp  —  Adq 


(A.l) 


where,  as  in  the  main  text,  e  is  the  specific  internal  energy  (energy  per  unit 
mass),  T  is  the  absolute  temperature,  s  is  the  entropy  (per  unit  mass),  P  is  the 
total  pressure,  p  is  the  density,  A  is  the  affinity  of  the  relaxation  process,  and 
q  is  the  degree  of  advancement  of  the  relaxation  mechanism.  As  was  pointed 
out  in  the  main  text,  three  independent  state  variables  are  required  for  relaxing 
fluids.  In  general,  p,  s,  and  q  are  used  as  the  independent  variables  in  this  work. 
Accordingly,  the  differential  of  the  internal  energy  e  may  be  expressed  as  follows: 


ds 


P.9 


de 

ds  +  — 
dp 


».9 


de 

ipJr  rq 


Comparison  with  Eq.  (A.l)  reveals  that 

de 


T  = 


ds 


A? 


and 


A~~  — 
A~  dq 


$,q 


p,3 


dq  . 


P,3 


Cross-differentiation  between  Eqs.  (A. 3)  and  (A. 4)  reveals  that 


(A. 2) 


( A .  3 ) 


(A. 4) 


(A. 5) 


dT_ 

dp 


3,q 


l_  dP_ 

p 2  ds 


(A. 6) 


Equations  (A. 3),  (A. 4),  and  (A. 6)  are  valid  for  constant,  or  frozen,  q. 

Equations  (A.l)  and  (A. 2)  are  now  examined  at  equilibrium.  At  equi¬ 
librium,  q  takes  on  its  equilibrium  value  q *,  and  A  is  equal  to  zero.  Furthermore, 
only  two  independent  thermodynamic  variables  are  required.  Thus,  at  equilib¬ 
rium,  Eqs.  (A.l)  and  (A. 2)  simplify  to 


de  ~  Tds  - 1 — -  dp 


(A. 7) 
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and 


.  Oc 

<fe  =  — 

ds 


J  de 
ds  4-  — 
dp 


dp  . 


I  A  *) 


I P,q=q*  r  'a. q=q' 

Consequently,  the  equilibrium  definitions  of  the  temperature  and  pressure  are  the 
same  as  the  frozen  definitions,  Eqs.  (A. 3)  and  (A.l),  except  the  derivatives  are 
evaluated  at  q  —  qm , 

de  I 


T  ds 


P.fl  =  9* 


'•-'•5 


(A.  Id) 


a,V=9* 


( 'loss-differentiation  between  the  two  foregoing  relations  leads  to  the  equilibrium 
form  of  Eq.  (A. 6), 


dT 

dp 


l_  d_P_ 
_  P1  ds 


a,q=q 


p,i=q 


Proceeding  as  above,  more  relations  may  be  obtained  from  the  C!i 
'quation  expressed  in  terms  of  the  enthalpy,  h  =  c  +  — , 


dh  —  Tds  +  -  dP  —  A  dq 
P 


(Add) 


Expressing  the  enthalpy  as  a  differential  quantity  in  terms  of  P.  s.  and  1/  yields 


dh  = 


dh 

ds 


ds  + 


P.t 


dh 

dP 


dP  + 


dh_ 

dq 


dq  . 


P,3 


Comparison  of  the  above  relation  with  Eq.  (A.  12)  reveals  that 


T=  dh 
ds 


P.O 


p  = 


dP 

dh 


and 


,  dh 
A  =  ~^q 


(A- 13) 

(A.l  1) 

(Add) 


P,3 
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As  with  the  Gibbs  equation  expressed  in  terms  of  the  internal  energy,  Eqs.  (A. 13) 
and  (A. 14)  may  be  evaluated  at  equilibrium,  that  is,  at  q  =  q'.  This  leads  to 


dh 

ds 


(A. 16) 


(A. 17) 


A-3  Definitions  of  Fundamental  Thermodynamic  Fluid  Properties 

In  this  section,  some  fundamental  thermodynamic  fluid  properties  are 
defined.  First  to  be  defined  are  the  frozen  forms  of  the  specific  heats  at  constant 
pressure  and  constant  volume,  the  isothermal  and  isentropic  bulk  moduli,  and 
the  coefficient  of  volume  expansion  (thermal  expansion  coefficient).  They  are, 
respectively, 

r°°  =  —  -T—  (AKI 

CP  *  O'T*  1  ( A .  K.  ) 


c“-ar|Pi  TdTPt 

Ccc  =  =T  — 

v '  dT  n ,, 


(A. 19) 


Kr  =  p 


(A.  30) 


.-CO  _  „  9P 

Ls=pTp 


(A. 21 


(A. 22) 


Q°°  =  _I  dp.  (A  ■>•>) 

The  definitions  of  the  equilibrium  forms  of  the  above  properties,  which  are  de¬ 
noted  by  the  superscript  u,  are  identical,  except  that  the  partial  derivatives  are 
evaluated  at  q  =  qm.  As  noted  in  the  introduction,  if  the  fluid  does  not  contain  a 
relaxation  mechanism,  the  need  for  the  superscript  is  gone,  and  it  is,  therefore, 
dropped. 

The  frozen  and  equilibrium  forms  of  the  ratio  of  specific  heats  are 


oo  s  _P_ 

'  rOO 

V 
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and 


7°  =  i 

7  c° 


(A. 21) 


It  can  be  shown  (see,  for  example,  Thompson  (1972,  p.  63))  that  q00  is  equal 
to  the  ratio  of  the  squares  of  the  frozen  isothermal  and  isentropic  sound  speeds, 
that  is,  that 


dP 

dP 


=  7°° 


dP_ 

dp 


(A.  25) 


T,q 


A  similar  relation  holds  for  the  equilibrium  values.  This  relation  is  used  in  the 
nondimensionalization  of  an  alternative  form  of  the  state  equation  that  involves 
the  temperature.  (See  the  derivation  that  precedes  Eq.  (B.23)  on  page  110.) 

A-4  Thermodynamic  Identities 

In  this  section,  relations  between  the  thermodynamic  fluid  properties  are 
developed.  The  first  relation  developed  is 


T(c°°a^°)2 


L-JLU 

c  <?) 


i 


(A. 26) 


Equation  (A. 26)  is  used  to  combine  the  entropy  and  state  equations  (see  Eq.  (2.95) 
on  page  28).  Our  starting  point  is  the  temperature  expressed  as  a  perfect  dilfer- 
ential, 

.  ,  arl 

ds  +  ~K~ 
dp 


Os 


P.9 


S,q 


,  dT 
dp  +  tt- 
dq 


dq 


P,» 


Dividing  through  by  ds  and  applying  constant  P  and  q  yields 


dr 

Ih 


P.9 


dT 

ds 


P.9 


dT 
+  dp 


S.q 


dp_ 

ds 


(A. 27) 


P.9 


By  applying  Eq.  (A. 6)  and  the  defining  relation  of  the  frozen  sound  speed, 
Eq.  (2.60),  to  the  following  identity, 


dp 

ds 

dP 

ds 

P,qdP 

dp 

P,q  ' 

».9 

hui  rearranging  and  inserting  the  results  into  Eq.  (A. 27),  we  obtain  the  following: 
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dT 

ds 


P.9 


dT 

ds 


P.9 


(£)' 


dT 

dp 


S.9J 
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Application  of  the  definitions  of  the  thermal  expansion  coefficient  and  the  specific 
heats  at  constant  pressure  and  volume,  Eqs.  (A. 22),  (A. 18),  and  (A. 19),  respec¬ 
tively,  leads  to  the  desired  result,  Eq.  (A. 26).  Use  of  an  identical  procedure  leads 
to  the  equilibrium  counterpart  of  Eq.  (A. 26),  namely, 


T(c' 


( i  u 
,0  0  \  2  I  1  1 

a’]  U  -  5, 


1 


(A. 28) 


For  nondimensionalizing  the  state  equation  [see  the  derivation  of 

is  required.  Our  starting  point 


Eq.  (B.20)  on  page  109],  an  expression  for 
is  the  expansion  of  frozen  specific  heat  at  constant  pressure, 


dP 

ds 


P,<7 


dh 

dh 

ds 

dp 

dT 

p,i 

ds 

P.q 

p.,  dT 

P.7 


Noting  the  definitions  of  the  temperature  and  the  thermal  expansion  coefficient, 
Eqs.  (A. 13)  and  (A. 22),  respectively,  we  see  that  the  above  relation  becomes 


(A. 29) 


P.7 


Applying  the  defining  relation  of  the  frozen  sound  speed,  Eq.  (2.60),  to  the  fol¬ 
lowing  identity, 

ds 

dp 


dp 

dP 

P,q  ^ 

ds 

3,q 

=  -1 


P>9 


and  then  rearranging  and  inserting  the  result  into  Eq.  (A. 29)  leads  to  the  desired 
result: 


dP_ 

ds 


p(c°°)2a™T 


Ip, 7  P 

The  equilibrium  counterpart  to  Eq.  (A. 30)  is 


dP_ 

ds 


p{c°)2a°T 


(A. 30) 


(A -31 ) 


P,7  =  7* 


For  nondimensionalizing  an  alternative  form  of  the  state  equation,  an 


expression  for 


dP  | 
dT 


P.7 


is  required.  (See  the  derivation  that  precedes  Eq.  (13.23)  on 


page  110.)  The  following  identity  is  our  starting  point: 

T.7 


dP 

dT 
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Using  the  definition  of  the  thermal  expansion  coefficient,  Eq.  (A. 22),  and  noting 
that  the  ratio  of  specific  heats  is  equal  to  the  ratio  of  the  isothermal  and  isentropic 
bulk  moduli,  Eq.  (A. 25),  we  readily  see  that 


OP 

Jrr 


P'Q 


^yOO 


(A.  32) 


The  equilibrium  counterpart  to  Eq.  (A. 32)  is 


i)P 

Jrr 


pa°(cy 


P,?=9 


(A. 33) 


A-5  Epilogue 

In  this  appendix,  the  fundamental  relations  from  thermodynamics  that 
are  commonly  used  in  acoustics  have  been  defined  (or  derived  as  need  be)  for 
fluids  with  a  single  relaxation  mechanism. 


APPENDIX  B 


ANALYSIS  OF  THE  NONDIMENSION AL  FORMS  OF  THE 
BASIC  EQUATIONS  FOR  A  HOMOGENEOUS, 
THERMOVISCOUS  FLUID 


B-l  Introduction 

Developed  in  this  appendix  are  nondimensional  forms  of  the  continu¬ 
ity,  momentum,  state,  and  entropy  equations  that  are  valid  for  a  homogeneous, 
thermoviscous  fluid  that  is  initially  both  quiet  and  irrotational.  When  nondimen- 
sionalizing,  it  is  assumed  that  the  signal  is  in  free  space,  far  from  any  boundaries. 
Also  developed  in  this  appendix  are  estimates  of  the  magnitude  of  the  nondimen¬ 
sional  coefficients  in  the  equations.  Moreover,  estimates  of  the  nondimensional 
signal  strength  are  obtained.  The  number  of  assumptions  required  for  small-signal 
acoustics  is  then  examined.  It  is  noted  that  only  one  fundamental  assumption  is 
required:  That  any  one  of  the  acoustic  field  variables  be  small.  It  then  follows 
that  the  others  are  also  small.  Next,  the  nondimensional  vorticity  equation  is  de¬ 
rived,  and  the  effect  of  the  irrotional  flow  assumption  is  considered.  The  notation 
used  in  this  appendix  differs  somewhat  from  that  in  the  main  text  because  the 
analysis  is  conducted  using  nondimensional  variables.  The  notation  is  introduced 
as  needed. 

B-2  Nondimensionalization  of  the  Equations  of  Motion 

Developed  in  this  section  are  the  nondimensional  forms  of  the  basic  equa¬ 
tions.  First,  the  dimensional  forms  of  the  equations  and  our  notation  are  intro¬ 
duced.  Then,  time  and  distance  are  nondimensionalized  such  that  their  deriva¬ 
tives  are  equal.  Next,  the  particle  velocity  and  the  excess  density  are  nondimen¬ 
sionalized  such  that  the  coefficients  of  the  nondimensional  continuity  equation 
are  unity.  The  procedure  is  repeated  for  the  excess  pressure  while  nondimen- 
sionalizing  the  momentum  equation.  The  procedure  is  repeated  again  for  the 
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entropy  anti  temperature  while  nondimensionalizing  two  different  forms  of  the 
state  equation.  Last,  the  nondimensional  entropy  equation  is  obtained. 

The  equations  of  motion 

The  dimensional  forms  of  the  continuity,  momentum,  entropy,  and  state 
equations  for  a  homogeneous,  thermoviscous  fluid  are1 


^  +  V.(pu)  =  0  , 


=  +  ^ u+(/iB  +  !/i)v(V.u)  , 


pi  ~ 
Dt 


i{pt  &)  —  i  o  +  a 


P~  Po 
Po 


d\Lj 

2  c  du  k 

di. 

3  ,}Wk 

4*  77T  D 

fp  -  Po\ 

2! 

{  po  ) 

(B.l) 

(B.2) 

(B.3) 

(B.4) 


where  p  is  the  density,  V  (=  iJL  -f  j  JL  +  kj|)  is  the  gradient  operator,  u  is  the 

A  A 

particle  velocity,  P  is  the  total  pressure,  T  is  the  absolute  temperature,  ,$  is  the 
entropy  (per  unit  mass),  p  is  the  shear  viscosity,  ps  is  the  bulk  viscosity,  n  is  the 
thermal  conductivity,  D / Dt  is  the  material  derivative, 


Dt  dt  v  ’ 


(B.5) 


and  the  constants  A  and  B  are  given  by 


*  ■  *  SIJ.  ■ 


n  — 


(B.G) 

(B.7) 


I  he  hat  symbol  “  indicates  the  dimensional  form  of  a  quantity.  The  static  value 
of  a  quantity,  denoted  by  the  subscript  0,  docs  not,  however,  need  a  hat  *  since 
it  is  inherently  dimensional. 

'As  in  (In-  main  text,  (lie  continuity,  momentum,  ami  entropy  equations  are  l\qs.  (12), 
(ia.ti),  ami  (19.5),  respectively,  in  the  hook  by  I.andau  and  Lifshitz  (1959). 
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Relationship  between  the  nondimensional  time  and  space  derivatives 

In  this  section,  time  and  distance  are  rendered  nondimensional  in  such  a 
way  that  the  nondimensional  time  and  space  derivatives  are  equal.  It  is  reasonable 
to  do  this  because  it  is  assumed  that  the  signal  is  far  from  any  boundaries.  Thus, 
the  only  dimensions  of  interest  are  those  associated  with  the  signal. 

Time  is  made  dimensionless  by  dividing  by  a  time  that  is  characteristic 
of  the  signal  tch  (=  4;  for  periodic  waves,  where  Cj  is  the  angular  frequency), 


t 

t  =  tTh 

Distance,  on  the  other  hand,  is  nondimensionalized  as  follows: 


(B.8) 


_ 

Xch 

where  z,  is  the  zth  cartesian  coordinate  and  xch  is  a  characteristic  distance.  The 
characteristic  distance  is  to  be  chosen  such  that  the  coefficients  in  the  classical 
wave  equation  are  unity.2  Expressing  the  classical  wave  equation  in  terms  of 
nondimensional  time  and  space  yields 


where  V  is  the  nondimensional  gradient  operator.  If  all  the  coefficients  in  the 
nondimensional  wave  equation  are  to  be  unity,  the  characteristic  distance  must  be 
related  to  the  characteristic  time  by  the  small-signal  sound  speed  Cq,  specifically, 


•Tch  —  Q)^c/i  • 


Thus,  the  definition  of  the  nondimensional  distance  becomes 

ii 

x,  =  ——  . 

Cfv'ch 


(119) 


Nondimensionalizing  the  continuity  equation 

1  he  continuity  equation,  Eq.  (B.l),  is  now  nondimensionalized.  As  part 
of  U  ic  procedure,  ; lie*  particle  velocity  and  excess  density,  which  appear  in  the 

2lt  w;»s  pointed  out  in  the  main  text  that  the  classical  wave  equation  is  a  good  model  of 
acoustic  propagation  for  the  important  case  of  small  signals  in  a  lossless  fluid. 
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continuity  equation,  are  also  rendered  nondimensional.  The  nondimensional  par¬ 
ticle  velocity  and  excess  density  are  defined  as  follows: 

u 

u  =  — 

M ch 


and 

5  =  Ll£l  , 

Pch 

where  uch  and  pch  are  the  characteristic  particle  velocity  and  density.  The  nondi¬ 
mensional  excess  density  is  referred  to  as  the  condensation.3  The  characteristic 
particle  velocity  and  density  are  to  be  chosen  in  such  a  way  that  all  the  coefficients 
of  the  nondimensional  continuity  equation  become  unity.  Expanding  the  conti¬ 
nuity  equation  and  expressing  it  in  terms  of  the  nondimensional  time,  distance, 
excess  density,  and  particle  velocity  yields 


|  +  ^5V-u  +  ^^(V.u)  +  -(u.?5)  =  0 

ot  Co  Co  pch  Co 


If  all  the  coefficients  in  the  above  relation  are  to  be  unity,  then  the  characteristic 
particle  velocity  and  density  must  be  defined  as  follows: 


^ ch  —  Q)  j 


Pch  ~  Po 


Accordingly,  the  definitions  of  the  nondimensional  particle  velocity  and  excess 
density  become 

u  =  -  (B.10) 


and 


5  = 


Co 

p-  Po 
Po 


(B.ll) 


Thus,  the  nondimensional  continuity  equation  becomes 


_  +  V((S+l)u|  =  0  .  (Ii.12) 

3Tho  symbol  5  has  been  chosen  to  denote  the  condensation  rather  than  the  more  commonly 
used  .s  because,  in  both  the  main  text  and  this  appendix,  s  is  used  to  denote  the  specific  entropy. 
Moreover,  the  condensation  is  not  used  outside  this  appendix. 
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Nondimensionalizing  the  momentum  equation 

The  momentum  equation,  Eq.  (B.2),  is  now  rendered  dimensionless. 
Only  one  new  variable,  the  pressure,  is  introduced  in  the  momentum  equation. 
The  pressure  is  nondimensionalized  such  that  the  coefficients  of  the  nondimen- 
sional  form  of  the  lossless  momentum  equation  become  unity.  The  viscous  loss 
terms  in  the  momentum  equation  are  ignored  in  this  procedure  because  losses  are 
assumed  to  be  of  secondary  importance  to  propagation. 

The  momentum  equation  is  first  rearranged  using  the  following  vector 
identities  (Gradshteyn  and  Ryzhik  1980,  Eq.  10.31.7): 

(u-V)u  =  (Vxu)xu  +  |Vu2  , 

V2u  =  V(V* u)  -  V x( V  Xu)  . 

Inserting  the  identities  into  the  momentum  equation,  Eq.  (B.2),  and  using  the 
definition  of  the  vorticity, 

J?  =  Vxu  ,  (B.13) 

leads  to 

p  +  (/?Xu)  +  |  Vti2j  +  VP  =  p  [VV(V*u)  —  (V x/?)]  .  (B.14) 

The  viscosity  number  V,  which  is  defined  in  Eq.  (2.82),  is  indicative  of  the  relative 
importance  of  the  bulk  viscosity  to  the  shear  viscosity.  In  the  next  section, 
estimates  of  the  magnitude  of  the  viscosity  number  and  the  other  nondimensional 
numbers  are  obtained. 

The  excess  pressure  is  nondimensionalized  using  a  characteristic  pressure 
Pck  that  is  to  be  determined, 


P-Po 

PcH 


Introducing  the  nondimensional  time,  distance,  particle  velocity,  excess  density, 
and  excess  pressure  into  Eq.  (B.14)  yields 


(5+i)  f^  +  (r?xu)  +  |vu2)  +  f±  Vp'  =  st[(Vxf?) 

\Ot  J  Poc o 
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(The  kinematic  vicosity  v  is  defined  in  Eq.  (2.81).)  The  magnitude  of  the  Stokes 
number  indicates  the  importance  of  viscosity.  If  the  Stokes  number  is  very  small 
(small  viscosity  or  large  characteristic  time  (low  frequency)),  then  the  momentum 
equation  may  be  approximated  as 


(5+1»^  +  ^rp'  =  ° 


If  the  coefficients  of  both  terms  in  the  foregoing  relation  are  to  be  unity,  then  the 
characteristic  pressure  must  be  given  by 


PcH  —  PoCq 


Thus,  the  definition  of  the  nondimensional  pressure  becomes 


,  P-Po 

P  =  - 2~ 

Po4 


(B.ic) 


Note  that  p0cl  is  the  static  value  of  the  isentropic  bulk  modulus. 

The  nondimensional  form  of  the  momentum  equation  is,  accordingly, 

(5+1)f^  +  («Xu)  +  iVu2^  +Vp/  =  St(vV(V.u)-Vxf?]  .  (B.17) 


Nondimensionalizing  the  state  equation 

The  state  equation,  Eq.  (B.4),  is  now  nondimensionalized.  The  only 
new  variable  introduced  in  the  state  equation  is  the  entropy.  The  entropy  is 
to  be  nondimensionalized  in  such  a  way  that  the  coefficients  of  the  linearized 
Taylor  series  expansion  of  the  state  equation  are  unity.  Our  starting  point  is  the 
definition  of  the  nondimensional  entropy, 


s  -  So 
$ch 


(BIS) 


where  sck  is  the  characteristic  value  of  the  entropy  that  is  to  be  determined.  Sub¬ 
stituting  the  above  into  the  linearized  state  equation  noting  the  defining  relation 
for  the  sound  speed,  c2  =  ,  yields 


p'  =  s  + 


$ck 

PoCq 
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However,  evaluation  of  Eq.  (A. 31)  from  Appendix  A  at  static  conditions  shows 
that 


dP_ 

ds 


P<%ap  T0 


p/  o 


-Po 


where  aPo  and  cPq  are,  respectively,  the  thermal  expansion  coefficient  and  the  spe¬ 
cific  heat  at  constant  pressure  evaluated  at  static  conditions.  Thus,  the  linearized 
state  equation  becomes 


ap0  To  , 


P  =  S  -f  sch 


The  coefficient  of  the  linear  entropy  term  is  unity  if  the  characteristic  entropy  is 
defined  as  follows: 

aPoi° 

Thus,  the  definition  of  the  nondimensional  entropy  and  the  nondimensional  equa¬ 
tion  of  state  become 

a„  To 

Po—  (s-s0)  (B.19) 


s  = 


"Po 


and 


p'  _s  +  s>  +  L!ls2  + 


(B.20) 


where  B/ A  is  a  nondimensional  number  referred  to  as  the  parameter  of  nonlin¬ 
earity  (Beyer  1974,  p.  99). 

Note  that  for  the  special  case  of  an  ideal  gas,  the  state  equation  is 
P  =  pRT,  where  R  is  the  gas  constant.  The  coefficient  of  thermal  expansion 
is,  accordingly, 

1 


a»  =  - 


(B.21) 


Thus,  for  the  special  case  of  an  ideal  gas,  it  is  seen  that  ctp<jTo  =  1  and  sck  =  cPq. 


Nondimensionalizing  an  alternative  state  equation 

Before  nondimensionalizing  the  entropy  equation,  an  alternative  form  of 
the  state  equation  is  nondimensionalized, 

P  =  P{p,  T)  . 
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The  reason  for  nondimensionalizing  this  equation  before  the  entropy  equation  is 
that  it  motivates  the  nondimensionalization  of  the  temperature.  The  nondimen¬ 
sional  temperature  is  defined  as  follows: 


r'  =  tzh 

~  Tch  ’ 

where  Tch  is  a  characteristic  temperature  that  is  to  be  determined  such  that  the 
coefficients  in  the  linearized  Taylor  series  expansion  of  the  alternative  state  equa¬ 
tion  are  unity.  Inserting  the  nondimensional  temperature  into  the  aforementioned 
equation  and  then  rearranging  yields 


P  =  ~ 2 


7V  0 


5+^4 

po4 


r 


But,  evaluating  Eqs.  (A. 32)  and  (A. 25)  from  Appendix  A 
and  inserting  them  into  the  above  leads  to 


at  static  conditions 


IP  =  S  +  apTchT' 


The  coefficient  of  T'  becomes  unity  if  the  characteristic  temperature  is  defined  as 
follows: 


Thus,  the  definition  of  the  nondimensional  temperature  and  the  nondimensional 
form  of  the  linearized  alternate  state  equation  become 


T'  =  aPo(T~T0) 

(B.22) 

7 P'  =  s+r  . 

(B.23) 

Nondimensionalizing  the  entropy  equation 


The  entropy  equation,  Eq.  (B.3),  is  now  nondimensionalized.  No  new 
variables  are  introduced  in  the  entropy  equation,  and  the  nondimensional  form 
of  the  entropy  equation  may,  therefore,  be  obtained  directly, 


(5+‘)<r  +  ,»w  = 

slV/“Uvi'r  + 


^f(V.u)2  +  i 


'  dn,  du, 
— — -  -f  - — - 
Kdxj  dx, 


(B-21) 


where  Pr  is  the  Prandtl  number,  which  is  defined  in  Eq.  (2.98).  The  Prandtl 
number  is  an  indicator  of  the  relative  importance  of  viscosity  to  heat  conduction. 
Note  that  the  coefficient  aPQcl/cPo  is  also  a  nondimensional  number. 

B-3  Estimates  of  the  Magnitude  of  the  Nondimensional  Coefficients 
and  Signal 

Estimates  of  the  magnitude  of  the  nondimensional  coefficients  in  the 
basic  equations  are  important  for  the  ranking  of  terms.  Estimates  of  the  following 
nondimensional  numbers  are,  therefore,  obtained:  (1)  the  Stokes  number,  (2)  the 
Prandtl  number,  (3)  the  viscosity  number,  and  (4)  the  nondimensional  coefficient 
ap0co/cp0-  Also  examined  is  the  magnitude  of  the  nondimensional  signal  by  way 
of  the  nondimensional  particle  velocity. 

The  nondimensional  number  that  indicates  the  effects  of  viscosity  is  the 
Stokes  number,  which  is  defined  in  Eq.  (B.15).  The  Stokes  number  is  proport  ional 
to  the  frequency  and  has  the  range  0  <  St  <  oo.  Its  value  is,  therefore,  estimated 
for  a  relatively  high  frequency,  1  MHz.  Estimates  of  Stokes  number  for  water 
and  air  at  20°C  and  1  atm  are  2.9  x  10~6  and  8.0  x  1 0— 4 ,  respectively.'  It  is 
apparent  that,  for  this  temperature,  pressure,  and  frequency,  air  is  more  viscous 
than  water.  Estimates  of  the  Stokes  number  at  1  MHz  for  a  variety  of  other  fluids 
are  listed  by  Truesdell  (1953,  Table  3.1). 

The  nondimensional  number  that  indicates  the  relative  importance  of 
viscosity  to  heat  conduction  is  the  Prandtl  number,  which  is  defined  in  Eq.  (2.98). 
From  its  definition  it  is  clear  that  the  Prandtl  number  has  the  range  0  <  Pr  <  oo. 
Estimates  of  the  Prandtl  number  cited  by  Truesdell  (1953,  Table  3.1)  indicate  it  is 
close  to  unity  for  many  gases,  but  generally  larger  for  liquids.  In  fact,  for  glycerin, 
Truesdell  cites  a  Prandtl  number  of  1000.  The  AIP  Handbook  (1977,  p.  2-203) 
cites  the  Prandtl  number  for  air  at  20°C  and  1  atm  as  0.71.  In  water  for  similar 
conditions,  the  Prandtl  number  is  about  7.O.5  Thus,  in  air  at  this  temperature 
and  pressure,  heat  conduction  is  a  slightly  more  important  loss  mechanism  than 

4The  estimate  for  water  is  based  on  density  and  sound  speed  data  taken  from  Kinsler  et 
al.  (1982,  p.  462)  and  on  shear  viscosity  data  from  the  CRC  Handbook  (1984,  p.  F-37):  pn  = 
998  kg/m3,  co  =  1481  m/s,  p  =  1.002  centipoises  (100  centipoises  =  1  poise  =  0.1  Pa  s). 
The  estimate  for  air  is  based  on  the  following  data  taken  from  Kinsler  et  al.  (1982,  p.  463): 
Po  =  1.21  kg/m3,  Co  =  343  m/s,  p  =  0.0000181  Pa  s. 

5This  value  is  based  on  shear  viscosity  and  specific  heat  at  constant  pressure  data  from 
the  CRC  Handbook  (1984,  p.  F-37  and  D-175)  and  on  coefficient  of  thermal  conduction  data 
from  the  AIP  Handbook  (1972,  p.  4—151):  p  =  1.002  centipoises,  cPo  =  4.2  x  103  J/kg°(l,  and 
k  =  0.6  YV/m°K. 
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viscosity,  whereas  in  water,  heat  conduction  apparently  plays  a  minor  role  in 
attenuation. 

Another  nondimensional  number  indicative  of  loss  mechanisms  is  the 
viscosity  number,  which  is  defined  in  Eq.  (2.82).  The  viscosity  number  indicates 
the  relative  importance  of  bulk  viscosity  to  shear  viscosity.  The  range  of  the 
viscosity  number  is  ~  <  V  <  oo.  Truesdell  cites  the  viscosity  number  for  air  at 
standard  temperature  and  pressure  as  1.9  and  for  water  at  15  °C  and  1  atm  as 
4.4.  Clearly,  bulk  viscosity  plays  a  more  significant  role  in  water  than  in  air. 

A  nondimensional  number  that  arose  in  the  entropy  equation  is  ap^j  cp^. 
Note  that  for  the  special  case  of  an  ideal  gas,  the  following  relations  hold: 
cPo  =  R~)/{ 7  —  1),  oPq  =  1  /T0,  and  c^  =  yRT0.  Thus,  we  see  that 


Since  the  value  of  7  for  ideal  gases  lies  in  the  range  1  <  7  <  5/3,  the  range  of 
oPoCq/ci>o  for  an  ideal  gas  is  0  <  qPqcI/cPq  <  2/3.  For  air  at  20°C  and  1  atm,  7 
is  about  1.4,  and  the  value  of  nPocl/cPo  in  air  is  therefore  about  0.4.  The  value 
of  OpoCo/cPo  in  water  at  20°C  and  1  atm  is  approximately  0.11. 6 

The  acoustic  Mach  number  is  a  nondimensional  number  that  indicates 
the  magnitude  of  the  signal.  The  acoustic  Mach  number  is  used  to  help  estimate 
the  magnitude  of  the  terms  in  the  basic  equations.  It  is  defined  as  the  ratio  of 
the  magnitude  of  the  peak  particle  velocity  to  the  small-signal  sound  speed, 


M  = 


|Upeak| 

CO 


(B.25) 


which,  in  terms  of  the  nondimensional  particle  velocity,  is  just  liip^ij.  The  value 
of  M  for  sound  levels  typical  of  a  normal  conversion  is  approximately  1.4  x  10-7, 
whereas  3  m  from  an  operating  jet  engine,  M  is  approximately  1.4  x  10~3.  By 
definition,  the  acoustic  Mach  number  is  independent  of  frequency.7 

"This  valid  is  based  on  the  previously  cited  data  for  the  sound  speed  and  specific  heat  at 
constant  pressure  as  well  as  thermal  expansion  coefficient  data  from  the  CRC  Handbook  (1984, 
P  F  4):  ttPo  =  206  x  10_6/°K. 

t  he  values  cited  are  based  on  the  values  of  60  dD  and  140  dB  re  20  /iPa,  respectively,  that 
are  cited  by  Pierce  (1981,  p.  62). 
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B-4  The  Number  of  Assumptions  Required  for  Small-Signal  Acous¬ 
tics 


From  the  signal  magnitude  estimates  above,  it  is  clear  that  for  ‘normal' 
acoustics,  the  acoustic  Mach  number  is  a  very  small  number  indeed.  In  fact, 
it  is  small  even  in  the  case  of  very  loud  sounds.  The  magnitude  of  the  other 
nondimensional  field  variables  in  the  basic  equations,  namely  <5,  p\  T\  and  .s', 
is  now  examined  for  the  case  of  small  Mach  number,  A/  < Cl.  It  turns  out 
if  the  acoustic  Mach  number  is  small,  then  the  acoustic  field  variables  are  of 
order  A/  or  of  higher  order.  Note  that  this  result  implies  that  the  number  of 
assumptions  required  for  small-signal  acoustics  is  one,  namely,  that  any  one  of 
the  field  variables  be  small. 

We  start  by  examining  the  nondimensional  continuity  equation, 
Eq.  (B.12),  with  the  expectation  of  obtaining  an  estimate  of  the  magnitude  of  the 
condensation.  The  three  possible  solutions  are  that  (1)  |5|  A/,  (2)  |«S|  3>  A/, 

and  (3)  I.S’I  ~  A/.  The  expanded  form  of  Eq.  (B.12)  is 

dS 

—  +  V*  u  -f  S  V*  u  +  u-V5  =  0 
at 

Recall  that  the  nondimensional  time  and  space  derivatives  are  equal.  An  exami¬ 
nation  of  the  above  equation  indicates  if  the  particle  velocity  is  of  order  A/,  then 
second  term  is  of  order  A/,  and  the  third  and  fourth  terms  are  of  order  |S|A/. 
The  first  term  is  of  order  |<S|.  Of  the  three  possible  solutions,  the  only  one  that 
balances  the  equation  is  the  third.  In  conclusion,  it  is  noted  that  if  A/  C  1 ,  then 
<S  1  also.  It  is  further  noted  that  if  M  1,  S  is  of  order  A/. 

A  similar  analysis  of  the  momentum  equation,  Eq.  (B.17),  is  now  con¬ 
ducted,  and  the  magnitude  of  the  nondimensional  pressure  is  determined.  First 
note  that  the  definition  of  the  vorticity,  Eq.  (B.13),  indicates  that  if  A/  <C  1 ,  then 
the  magnitude  of  the  vorticity  is  also  very  small,  \  fl\  <  1.  (It  turns  out,  as  is 
seen  later,  that  the  vorticity  is  actually  quite  a  bit  smaller  than  A/.)  Also  recall 
that  the  Stokes  number  is  generally  a  small  number.  Ignoring  the  terms  that  are 
products  of  small  quantities  simplifies  Eq.  (B.17), 


Ou 

Ik 


+  Vp'  =  o 


Recalling  that  the  nondimensional  time  and  space  derivatives  are  equal,  it  may 
be  seen  that  if  A/  1,  then  l//|  <C  1  also  Thus,  p'  is  of  order  A7,  if  A/  «  1. 

It  turns  out  that  initially  assuming  that  any  one  of  p',  u,  or  S  is  small 
(much  less  than  unity)  implies  that  the  other  two  are  also  small.  Thus,  the 
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number  of  assumptions  required  for  small-signal  acoustics  is  one,  namely,  that 
one  of  p',  u,  or  S  be  assumed  small.  That  the  other  two  are  small  follows  directly. 

Estimates  of  the  magnitude  of  the  nondimensional  temperature  and  en¬ 
tropy  fluctuations  are  now  obtained.  An  examination  of  the  alternate  state  equa¬ 
tion,  Eq.  (11.23),  indicates  the  nondimensional  temperature  fluctuation  T'  is  of 
order  (7  —  1)A/.  Simplifying  the  nondimensional  entropy  equation,  Eq.  ( 13.2-1 ), 
by  ignoring  products  of  small  quantities  leads  to 


ds' 

dt 


=  D°P,,r<wJr 


Thus,  s'  is  of  order  j^oPo7o(7  —  1)A/. 

In  this  section,  it  was  noted  that  if  M  <  1,  then  the  excess  pressure  and 
density  are  of  order  M .  Moreover,  it  was  noted  that  the  number  of  assumptions 
required  for  small-signal  acoustics  is  just  one.  It  was  further  noted  that  if  Af  <  1, 
then  T'  and  s'  are  of  orders  (7  —  1 ) A/  and  |j;crpo7o(7  —  1 ) A/,  respectively. 


B-5  The  Nondimensional  Vorticity  Equation  and  the  Irrotational 
Flow  Assumption 


The  nondimensional  form  of  the  vorticity  equation  is  of  interest  because 
it  will  assist  us  in  understanding  the  irrotational  flow  assumption.  The  nondi¬ 
mensional  vorticity  equation  is  formed  by  taking  the  curl  of  the  nondimensional 
momentum  equation,  Eq.  (B.17).  This  is  readily  done  using  index  notation  such 
as  that  described  in  the  book  by  Panton  (1984,  Chap.  3). 8  The  result  is 


DO 

UT 


=  (tt-V)u  -  !?( V*  u) 


+ 


1 

(5+1) 

St 

5+  1  . 


7  ( Vp' X VS)  -  St  V  [■ V5x  V( V.  u)] 

(^T)[V5x(Vxf?)]  . 


(13.26) 


1  his  is  an  equivalent  nondimensional  form  of  Eq.  (2.46)  in  the  book  by  Thompson 
( 1972).  I  lie  generation  and  transport  of  vorticity  is  now  briefly  examined.  For  a 
more  complete  discussion,  the  reader  is  referred  to  Thompson  and  the  references 
contained  therein. 

*  Pant  1  >11  (IDS),  i'11.  13. 3. 5)  uses  index  notation  to  develop  the  dimensional  form  of  I  lie 
vorticilv  ecpi.iiion  vaiid  for  incompressible  fluids. 
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Consider  a  lossless  fluid.  Note  that  the  Stokes  number  is  zero  for  an 
inviscid  fluid,  and  that,  in  a  uniform  lossless  fluid,  the  gradient  of  the  entropy  is 
zero  as  well.  Thus,  the  term  Vp'  X  VS,  which  is  proportional  to  Vp'x  Vs'  (re¬ 
call  that  |Vp'xVp'|  =  0),  is  zero.  Accordingly,  the  vorticity  equation  simplifies 

to  Dfi 

—  =  u-IJ(V.u)  . 

Note  if  the  vorticity  is  initially  zero,  the  vorticity  is  always  zero  because  rio 
vorticity  is  generated  in  the  absence  of  viscosity  or  entropy  gradients. 

In  the  case  of  a  heat-conducting  but  inviscid  fluid,  the  Stokes  number  is 
still  zero,  but  the  gradient  of  the  entropy  is  nonzero.  Equation  (B.26)  therefore 
simplifies,  but  not  quite  as  much  as  above, 

^  =  (O.V)u-fi(V-u)--j-i-j5(Vp'xVS)  . 

The  new  term  is  a  vorticity  generation  term  that  is  nonzero  at  time  zero  in  t  he  case 
of  an  initially  irrotational  fluid.  Thus,  even  in  the  absence  of  viscosity,  vorticity 
can  develop;  it  does  not,  however,  diffuse.  That  requires  viscosity  and  is  discussed 
below.  If  the  fluid  is  initially  irrotational,  the  new  term  generates  vorticity  at  a 
rate  proportional  to  Vp'xVs',  which,  according  to  our  ordering  system,  is  very 
small — of  the  order  f;  GpToil  ~  1)M2.  (It  may,  in  fact,  be  even  smaller  because 
a  cross-product  is  proportional  to  the  sine  of  the  angle  between  the  vectors.)  The 
other  vorticity  generation  (or  destruction)  terms,  such  as  (/?*V)u,  are  of  order 
M  higher  than  the  order  of  17.9  Thus,  the  vorticity  that  these  terms  generate 
(or  destroy)  is  also  very  small.  However,  our  nondimensionalization  scheme,  and 
hence  our  ordering  system,  rest  on  the  assumption  that  the  signal  is  far  from 
any  boundaries.  Accordingly,  all  that  may  be  stated  is  that  the  generation  of 
vorticity  in  the  middle  of  the  fluid  is  small.  Since  the  entropy  gradient  is  much 
sharper  at  boundaries  (due  to  the  presence  of  the  thermal  boundary  layer,  which 
has  not  been  discussed),  vorticity  generation  at  the  boundaries  is  larger. 

In  a  thermoviscous  fluid,  all  terms  in  the  vorticity  equation  must  be 
retained.  Now  in  the  case  of  an  initially  irrotational  fluid,  a  second  vorticity 
generation  term  is  nonzero  at  time  zero,  St  V(^r)  (VSxV(V.u)].  The  order 
of  this  term  is  StVA/2,  again,  a  very  small  number  in  the  middle  of  the  fluid. 
Once  more  it  is  important  to  note  that  near  boundaries  the  gradient  terms  are 
much  larger,  and  vorticity  generation  is  also  proportionally  larger. 

9Panton  (1984,  Sec.  13.5)  discusses  this  term  at  length  and  points  out  that  it  represents  the 
generation  or  destruction  of  vorticity  by  straining  vortex  lines. 
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Since  the  production  of  vorticity  in  the  middle  of  the  fluid  is  so  small  in 
comparison  to  that  near  a  boundary,  it  is  of  interest  to  understand  the  transport 
of  vorticity  from  the  boundaries  into  the  fluid.  To  that  end  the  leading-order 
terms  in  the  vorticity  equation  are  examined, 

^  =  stv2r?  . 


Before  solving  this  equation,  we  place  it  in  dimensional  form, 

dt 

Assuming  a  solution  of  the  form  fl  =  and  substituting  into  the  above 

leads  to 

k2=-J-  , 

v 

where  uj  is  the  angular  frequency  and  k  is  the  wave  number.  The  wave  number 
may  be  defined  as  follows: 

A  tlj 

k  - - ja  , 

Cph 

where  cph  is  the  phase  velocity  and  a  is  the  attenuation  coefficient  that  must 
be  positive  for  stability  reasons.  Solving  for  the  attenuation  coefficient  and  the 
phase  velocity  noting  that  yf]  =  ±(1  +  j)/y/2  yields 

1  1 

a  =  —  =  -  , 

Cph  v 

where  the  negative  root  was  chosen  so  that  a  is  positive  and  8  is  defined  as 


Thus,  the  assumed  solution  becomes 

Q  =  Ae-i/‘eJ“(i-i/c^h)  . 

Note  that  the  vorticity  decays  to  1/e  of  its  original  value  in  the  distance  8. 

Estimates  of  the  distance  8  for  water  and  air  may  be  obtained  using  the 
viscosity  data  in  footnote  4  on  page  111.  At  a  frequency  of  1  Mllz,  8  is  about 
5.7  x  It)-7  m  in  water  and  2.2  x  lO-6  m  in  air.  Clearly,  even  on  a  wavelength 
scale,  the  distance  the  vorticity  penetrates  is  very  small. 
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In  this  section,  it  was  noted  that  the  vorticity  generation  in  the  middle 
of  the  fluid  is  very  small — of  order  StVA/2  and  that  the  vorticity  generated  at 
the  boundaries,  for  which  no  estimates  were  provided,  is  restricted  to  a  very  thin 
layer  near  the  boundary.  It  is  therefore  concluded  that,  away  from  boundaries, 
vorticity  may  be  neglected  for  acoustic  problems  in  which  terms  of  order  StV.U2 
and  higher  are  neglected. 

B-6  Summary 

In  this  appendix,  the  continuity,  momentum,  state,  and  entropy  equa¬ 
tions  that  are  valid  for  a  homogenous,  thermoviscous  fluid  were  nondimcnsion- 
alized.  Also  obtained  were  estimates  of  the  magnitude  of  the  nondimensional 
coefficients  in  the  equations  and  of  the  nondimensional  field  variables.  The  num¬ 
ber  of  assumptions  required  for  small-signal  acoustics  was  then  analyzed.  It  was 
noted  that  the  number  is  one,  namely,  that  any  one  of  the  excess  pressure,  parti¬ 
cle  velocity,  or  excess  density  be  small.  Futhermore,  the  nondimensional  vorticity 
equation  was  obtained,  and  it  was  noted  that,  away  from  boundaries,  vorticity 
may  be  neglected  for  acoustic  problems  in  which  terms  of  order  StVM2  and  higher 
are  neglected. 


I 


BIBLIOGRAPHY 


Aanonsen,  S.  I.,  T.  Barkve,  J.  NazeTjptta,  and  S.  Tjptta  (1984).  “Distortion  and 
harmonic  generation  in  the  nearfield  of  a  finite  amplitude  sound  beam,” 
J.  Acoust.  Soc.  Am.  75,  749-768. 

AIP  Handbook  (1972).  American  Institute  of  Physics  Handbook ,  3rd  ed.,  edited 
by  D.  E.  Gray  (McGraw-Hill  Book  Co.,  Inc.,  New  York). 

Beyer,  R.  T.  (1974).  Nonlinear  Acoustics  (Naval  Ship  Systems  Command,  Dept, 
of  the  Navy,  Washington,  DC). 

Blackstock,  D.  T.  (1959).  “Normal  reflection  of  finite  amplitude  plane  waves  from 
a  rigid  wall,”  in  Proc.  Third  Intern.  Cong.  Acoust.,  Stuttgard,  1959,  edited 
by  L.  Cremer,  Vol.  I,  pp.  309-311  (Elsevier  Publ.  Co.,  Amsterdam,  1961). 

Blackstock,  D.  T.  (1972).  “Nonlinear  Acoustics  (Theoretical),”  in  American  In¬ 
stitute  of  Physics  Handbook,  3rd  ed.,  edited  by  D.  E.  Gray,  pp.  3-183  to 
3-205  (McGraw-Hill  Book  Co.,  Inc.,  New  York). 

Cotaras,  F.  D.  (1985).  “Nonlinear  Effects  in  Long  Range  Underwater  Acoustic 
Propagation,”  Master’s  Thesis,  The  University  of  Texas  at  Austin.  Ab¬ 
stract:  J.  Acoust.  Soc.  Am.  79,  2095  (1986).  Reprinted  under  the  same 
title  as  Tech.  Report  ARL-TR-85-32,  Applied  Research  Laboratories,  The 
University  of  Texas  at  Austin  (ADA  166  492). 

Cotaras,  F.  D.,  and  D.  T.  Blackstock  (1987).  “Refraction  of  finite  amplitude 
sound  at  a  plane  fluid-fluid  interface,”  J.  Acoust.  Soc.  Am.  82,  Si  1(A). 

Courant,  R.,  and  K.  O.  Friedrichs  (1976).  Supersonic  Flow  and  Shock  Waves, 
Appl.  Math.  Sci.  21  (Springer- Verlag,  New  York)  (Reprint  of  1948  edition 
published  as  Vol.  1  of  Pure  and  Appl.  Math,  by  Interscience  Publishers, 
New  York). 

CRC  Handbook  (1984).  Chemical  Rubber  Company  Handbook  of  Chemistry  and 
Physics,  65th  ed.,  edited  by  R.  C.  Weast  (CRC  Press,  Inc.,  Boca  Raton, 
FL). 


118 


119 


Gradshteyn,  I.  S.,  and  I.  M.  Ryzhik  (1980).  Table  of  Integrals,  Series,  and  Prod¬ 
ucts,  Corrected  and  Enlarged  Edition  (Academic  Press,  Inc.,  New  York). 

de  Groot,  S.  R.,  and  P.  Mazur  (1984).  Non-Equilibrium  Thermodynamics  (Dover 
Publications,  Inc.,  New  York)  (originally  published  in  1962  by  North- 
Holland  Publishing  Co.,  Amsterdam,  The  Netherlands). 

Feng  Shaosong  (1983).  “The  reflection  and  refraction  of  a  large  amplitude  plane 
sound  wave  in  two  dimensions,”  Chin.  J.  Acoust.  2,  293-302. 

Hamilton,  M.  F.  (1986).  “Fundamentals  and  Applications  of  Nonlinear 
Acoustics,”  in  Nonlinear  Wave  Propagation  in  Mechanics,  edited  by 
T.  W.  Wright,  Am.  Soc.  Mech.  Engin.  AMD-Vol.  77,  Book  No.  1100346. 

Hamilton,  M.  F.,  and  D.  T.  Blackstock  (1988).  “On  the  coefficient  of  nonlinearity 
/3  in  nonlinear  acoustics,”  J.  Acoust.  Soc.  Am.  83,  74-77. 

Henderson,  L.  F.  (1966).  “The  reflection  of  a  plane  shock  wave  at  a  gaseous 
interface,”  J.  Fluid  Mech.  26,  607-637. 

Jahn,  R.  G.  (1956).  “The  refraction  of  shock  waves  at  a  gaseous  interface,” 
J.  Fluid  Mech.  1,  457-489. 

Jahn,  R.  G.  (1957).  “Transition  processes  in  shock  wave  interactions,”  J.  Fluid 
Mech.  2,  33-48. 

Kinsler,  L.  E.,  A.  R.  Frey,  A.  B.  Coppens,  and  J.  V.  Sanders  (1982).  Fundamentals 
of  Acoustics,  3rd  ed.  (John  Wiley  and  Sons,  Inc.,  New  York). 

Kreider,  D.  L.,  R.  G.  Kuller,  D.  R.  Ostberg,  and  F.  W.  Perkins  (1966).  An  In¬ 
troduction  to  Linear  Analysis  (Addison- Wesley  Publishing  Co.,  Reading, 
MA). 

Lamb,  H.  (1932).  Hydrodynamics ,  6th  ed.  (Dover  Publications,  Inc.,  New  York, 
1945). 

Landau,  L.  D.,  and  E.  M.  Lifshitz  (1959).  Fluid  Mechanics  (Pergamon  Press, 
Oxford,  England). 

Lighthill,  M.  J.  (1956).  “Viscosity  effects  in  sound  waves  of  finite  amplitude,” 
in  Surveys  in  Mechanics,  edited  by  G.  K.  Batchelor  and  R.  M.  Davies 
(Cambridge  University  Press,  Cambridge,  England). 


120 


Myers,  M.  K.  (1980).  “On  the  acoustic  boundary  condition  in  the  presence  of 
flow,”  J.  Sound  Vib.  71,  429-434. 

Naze  Tjotta,  J.,  and  S.  Tjotta  (1987).  “Interaction  of  sound  waves.  Part  I:  Basic 
equat  ions  and  plane  waves,”  J.  Acoust.  Soc.  Am.  82,  1425-1428. 

Naze  Tjotta,  J.,  and  S.  Tj0tta  (1988).  Personal  communication. 

Panton,  R.  L.  (1984).  Incompressible  Flow  (John  Wiley  and  Sons,  Inc.,  New 
York). 

Pierce,  A.  D.  (1981).  Acoustics:  An  Introduction  to  Its  Physical  Principles  and 
Applications  (McGraw-Hill  Book  Co.,  Inc.,  New  York). 

Polachek,  H.,  and  R.  J.  Seeger  (1951).  “On  shock-wave  phenomena;  refraction  of 
shock  waves  at  a  gaseous  interface,”  Physical  Review  84,  922-929. 

Prigogine,  I.  (1961).  Introduction  to  the  Thermodynamics  of  Irreversible  Pro¬ 
cesses,  3rd  ed.  (Interscience  Publishers,  a  division  of  John  Wiley  and  Sons, 
Inc.,  New  York). 

Qian  Zuwen  ( 19S2).  “Reflection  of  finite  amplitude  sound  wave  on  a  plane  bound¬ 
ary  of  a  half  space,”  Scientia  Sinica  (Series  A)  XXV,  492-501. 

Rudenko,  O.  V.,  and  S.  I.  Soluyan  (1977).  Theoretical  Foundations  of  Nonlinear 
Acoustics  (Consultants  Bureau,  New  York). 

Shu,  Kun-Tien,  and  J.  H.  Ginsberg  (1988).  “Oblique  reflection  of  a  finite- 
amplitude  dilatational  wave  in  an  elastic  half-space,”  J.  Acoust.  Soc. 
Am.  84,  S5(A). 

Thomas,  G.  B.,  and  R.  L.  Finney  (1979).  Calculus  and  Analytic  Geometry,  5th 
ed.  (Addison-Wesley  Publishing  Co.,  Reading,  MA). 

Thompson,  P.  A.  (1972).  Compressible- Fluid  Dynamics  (McGraw-Hill  Book  Co., 
Inc.,  New  York). 

Trucsdell,  C.  (1953).  “Precise  theory  of  the  absorption  and  dispersion  of  forced 
plane  infinitesimal  waves  according  to  the  Navier-Stokes  equations,” 
./.  Rat.  Mcch.  Anal.  2,  643-741. 

Van  Buren,  A.  L.,  and  M.  A.  Breazeale  (1968).  “Reflection  of  finite-amplitude 
ultrasonic  waves.  I.  Phase  Shift,”  J.  Acoust.  Soc.  Am.  44,  1014-1020. 


121 


Van  Buren,  A.  L.,  and  M.  A.  Breazeale  (1968).  “Reflection  of  finite-amplitude 
ultrasonic  waves.  II.  Propagation,”  J.  Acoust.  Soc.  Am.  44,  1021-1027. 

Van  Wylen,  G.  J.,  and  R.  E.  Sonntag  (1976).  Fundamentals  of  Classical  Ther¬ 
modynamics,  2nd.  ed.,  SI  vers.  (John  Wiley  and  Sons,  Inc.,  New  York). 

Vincenti,  W.  G.,  and  C.  H.  Kruger  (1965).  Introduction  to  Physical  Gas  D  ynamics 
(John  Wiley  and  Sons,  Inc.,  New  York). 


3  January  1989 


Copy  No. 


1-3 


4-9 


10-21 


22 

23 


24 


25 


26 

27 

28 


DISTRIBUTION  LIST  FOR 
ARL-TR-89-1 

UNDER  ONR  GRANT  NO.  N00014-89-J-1109 


Office  of  Naval  Research 
Department  of  the  Navy 
800  North  Quincy  Street 
Arlington,  VA  22217-5000 
Attn:  L.  E.  Hargrove  (Code  1112) 

Director 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
Attn:  Code  2627 

Commanding  Office  and  Director 
Defense  Technical  Information  Center 
Building  5,  Cameron  Station 
5010  Duke  Street 
Alexandria,  VA  22314 

Naval  Postgraduate  School 
Monterey,  CA  93940 
Attn:  Technical  Library  (Code  0212) 
S.  L.  Garrett 

Naval  Research  Laboratory 
P.  O.  Box  8337 
Orlando,  FL  32856 
Attn:  D.  H.  Trivett 

Georgetown  University 
Physics  Department 
Washington,  D.C.  20057 
Attn:  W.  G.  Mayer 

Georgia  Institute  of  Technology 
School  of  Mechanical  Engineering 
Atlanta,  GA  30332 
Attn:  Y.  H.  Berthelot 
J.  H.  Ginsberg 
P.  H.  Rogers 


123 


Distribution  list  for  ARL-TR-89-1  under  ONR  Grant  No.  N00014-89-J-1109 
(cont’d) 


Copy  No. 


29 


30 

31 


32 


33 

34 


35 


3G 

37 

38-47 

4S 

49 


Johns  Hopkins  University 
Department  of  Mechanical  Engineering 
Baltimore,  MD  21218 
Attn:  A.  Prosperetti 

University  of  Mississippi 
Physics  Department 
University,  MS  39677 
Attn:  M.  A.  Breazeale 
L.  A.  Crum 

The  Pennsylvania  State  University 
Applied  Research  Laboratory 
University  Park,  PA  16802 
Attn:  A.  D.  Pierce 

University  of  Rochester 
Electrical  Engineering  Dept. 

Hopernan  Building 
Rochester,  NY  14627 
Attn:  E.  L.  Carstensen 
K.  J.  Parker 

Yale  University 
Department  of  Engineering 
Mason  Laboratory 
9  Mill  House  Avenue 
New  Haven,  CT  06520 
Attn:  R.  E.  Apfel 

Chief 

Defence  Research  Establishment  Atlantic 
P.  O.  Box  1012 

Dartmouth,  Nova  Scotia  B2Y  3Z7 
CANADA 

Attn:  Technical  Library 
D.  M.  F.  Chapman 
F.  D  Cotar  as 
I.  A.  Fraser 
H.  M.  Merklinger 


124 


Distribution  list  for  ARL-TR-89-1  under  ONR  Grant  No.  N00014-89-J-1109 
(cont’d) 

Copy  No. 

Vibron  Unlimited 
1720  Meyerside  Dr. 

Mississauga,  Ontario  L5T  1A3 
CANADA 

50  Attn:  Stephen  J.  Lind 

Technical  University  of  Denmark 
Industrial  Acoustics  Laboratory 
Building  352 
DK-2800  Lyngby 
DENMARK 

51  Attn:  L.  Bjprnp 

Osaka  University 
The  Institute  of  Scientific  and 
Industrial  Research 
8-1,  Mihogaoka,  Ibaraki 
Osaka  567 
JAPAN 

52  Attn:  A.  Nakamura 

The  University  of  Electro-Communications 

1-5-1,  Chofugaoka 

Chofu-shi 

Tokyo 

JAPAN 

53  Attn:  T.  Kamakura 

University  of  Bergen 
Department  of  Mathematics 
N-5014  Bergen-U 
NORWAY 

54  Attn:  Hans-Christen  Salvesen 

55  Kjell-Eivind  Frpysa 

Cambridge  University 
Department  of  Applied  Mathematics  and 
Theoretical  Physics 
Silver  Street 
Cambridge  CB3  9EW 
UNITED  KINGDOM 

56  Attn:  D.  G.  Crighton 


125 


Distribution  list  for  ARL-TR-89-1  under  ONR  Grant  No.  N00014-89-J-1109 
(ccnt’d) 


Copy  No. 

57 


58 

59 
CO 
G1 
62 

63 

64 

65 

66 
67 
GS 

69 

70 

71 

72 


Institute  of  Sound  and  Vibration  Research 

The  University 

Southampton  S09  5NH 

UNITED  KINGDOM 

Attn:  C.  L.  Morfey 

University  of  Bath 
School  of  Physics 
Claverton  Down 
Bath  BA2  7 AY 
UNITED  KINGDOM 
Attn:  H.  O.  Berktay 

Michalakis  A.  Averkiou  (ME:UT) 

Corinne  M.  L.  Darvennes  (ME:UT) 

Wang  S.  Lee  (ME:UT) 

Edel  Reiso  (ME:UT) 

David  T.  Blackstock  (ARL:UT) 

Charles  E.  Bradley  (ARL:UT) 

David  E.  Edmonds  (ARL:UT) 

Mark  F.  Hamilton  (ARL:UT) 

James  A.  Hawkins  (ARL:UT) 

John  D.  Sample  (ARL:UT) 

Deborah  A.  Sumina  (ARL:UT) 

James  A.  Ten  Cate  (ARL:UT) 

Jacqueline  Naze  Tjptta  (ARL:UT) 

Sigvc  Tjotta  (ARL:UT) 


126 


Distribution  list  for  ARL-TR-89-1  under  ONR  Grant  No.  N00014-89-J-1109 
(cont’d) 

Copy  No. 

73  Wayne  M.  Wright  (ARL:UT) 

74  Library  (ARL:UT) 

75-94  NAD  Reserve 


127 


